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Chapter 1 



Introduction 



1.1 Motivation 

This book is about nonlinear distortion in radio frequency (RF) power 
amplifiers (PAs). The purpose of the PA is to boost the radio signal to a 
sufficient power level for transmission through the air interface from the 
transmitter to the receiver. This may sound simple, but it involves solving 
several contradicting requirements, the most important of which are 
linearity and efficiency. Unfortunately, these requirements tend to be 
mutually exclusive, so that any improvement in linearity is usually 
achieved at the expense of efficiency, and vice versa. 

To avoid interfering with other transmissions, the transmission must 
stay within its own radio channel. If the modulated carrier has amplitude 
variations, any nonlinearity in the amplifier causes spreading of the 
transmitted spectrum (so-called spectral regrowth). This effect can be 
reduced by using constant-envelope modulation techniques that 
unfortunately have quite low data rate/bandwidth ratio. When using more 
efficient digital modulation techniques, the only solution is to design the 
amplifiers linear enough. 

The efficiency is defined as a ratio of the generated RF power and the 
drawn dc power. In modem radio telecommunication systems, the design of 
linear and efficient radio frequency power amplifier presents one of the 
most challenging design problems. In general, relatively high transmit 
power levels are needed, and the power consumption of the PA easily 
dominates over all other electronics and digital processing in a mobile 
terminal. Therefore, high efficiency is essential to extend the operation 
time of the terminals. In fixed-point wireless nodes (e.g., in base stations), 
efficiency is also important, because the transmitted power levels are 
essentially higher than in terminals. 
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1.2 Historical Perspective 

In first-generation systems, such as the Nordic Mobile Telephone (NMT) or 
Advanced Mobile Phone Service (AMPS), the RF signal was frequency 
modulated (FM). Highly efficient PAs are possible in FM systems because 
of the fact that no information is encoded in the amplitude component of 
the signal. Even so, the PA of a mobile phone consumed as much as 85% of 
the total system power at the maximum power level, thus limiting the on- 
time of the terminal. 

Unlike wired line communications, wireless systems must share a 
common transmission medium. The available spectrum is therefore limited, 
and so channel capacity (i.e., the amount of information that can be carried 
per unit bandwidth) is directly associated with profit. The demand for 
greater spectral efficiency was addressed by the development of second- 
generation systems, where digital transmission and time domain multiple 
access (TDMA) is used, where multiple users are time multiplexed on the 
same channel. For example, in the Global System for Mobile 
Communications (GSM), eight calls alternate on the same frequency 
channel, resulting in cost-effective base stations. The GSM modulation 
scheme retains constant envelope RF signals, but the need for smooth 
power ramp up and ramp down of the allocated time-slot transmissions 
imposes some moderate linearity requirements. This reduces the efficiency 
of the amplifier, but it is compensated by the fact that the PA in the mobile 
node is only active one-eighth of the time. This, together with the smart 
idling modes, allows GSM handsets to achieve very long operating times. 

The data transmission capacity of GSM is rather modest, so the 
obvious solution to increase the achievable bit rate was, as implemented in 
GSM-EDGE, to use several time slots for a single transmission and to 
replace the Gaussian minimum shift key (GMSK) modulation scheme with 
a spectrally more efficient 8-PSK that unfortunately has a varying 
envelope. So as wireless communication systems migrate towards higher 
channel capacity, more linear and, consequently, less efficient PAs have 
become the norm. 

Finally, the third generation wideband code-division multiple access 
(WCDMA) packs tens of calls on the same radio channel simultaneously, 
differentiated only by their unique, quasi-orthogonal spreading codes. This 
allows flexible allocation of data rates, while tolerance to fading is 
improved by increasing the signal bandwidth to nearly 4 MHz. The 
advantages offered by the WCDMA, however, come at the expense of more 
stringent requirements for the PA. The code-multiplexed transmission 
occupies a much larger bandwidth than in the previous systems, while 
exhibiting tremendous variations in amplitude. Furthermore, in WCDMA, 
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the mobile transmits on a continuous time basis. Designing an economical 
PA for these requirements is an enormous engineering challenge. 

The situation is not easier in the base stations, either, where the 
linearity requirements are tighter than in handsets. The trend is towards 
multicarrier transmitters where a single amplifier handles several carriers 
simultaneously, in which case the bandwidth, power level, and the peak 
power to average power ratio (crest factor) all increase. The efficiency of 
these kinds of power amplifiers is very low, and due to higher total 
transmitted power, this results in very high power dissipation and serious 
cooling problems. 



1.3 Linearization and Memory Effects 

The goal of this book is to improve the conceptual understanding needed in 
the development of PAs that offer sufficient linearity for wideband, 
spectrally efficient systems while still maintaining reasonably high 
efficiency. As already noted, efficiency and linearity are mutually exclusive 
specifications in traditional power amplifier design. Therefore, if the goal is 
to achieve good linearity with reasonable efficiency, some type of 
linearization technique has to be employed. The main goal of linearization 
is to apply external linearization to a reasonably efficient but nonlinear PA 
so that the combination of the linearizer and PA satisfy the linearity 
specification. In principle, this may seem simple enough, but several higher 
order effects seriously limit its effectiveness, in practice. 

Several linearization techniques exist, and they are reviewed in Chapter 
3; a much more detailed discussion can be found from [1-3]. Stated briefly, 
linearization can be thought of as a cancellation of distortion components, 
and especially as a cancellation of third-order intermodulation (IM3) 
distortion, and where the achieved performance is proportional to the 
accuracy of the canceling signals. Unfortunately, the IM3 components 
generated by the power amplifier are not constant but vary as a function of 
many input conditions, such as amplitude and signal bandwidth. Here, 
these bandwidth-dependent phenomena are called memory effects. 

Smooth, well-behaved memory effects are usually not detrimental to 
the linearity of the PA itself. If the phase of an IM3 component rotates 10° 
to 20°, or if its amplitude changes 0.5 dB with increased tone spacing in a 
two-tone test, it usually does not have a dramatic effect on the adjacent 
channel power ratio (ACPR, i.e., the power leaking to the neighboring 
channel) performance of a standalone amplifier, nor is it especially of 
concern if the lower ACPR is slightly different from the upper one. 
However, the situation may be quite different if certain linearization 
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techniques are used to cancel out the intermodulation sidebands; in fact, the 
reported performance of some simple techniques may actually be limited 
not by the linearization technique itself, but by the properties of the 
amplifier - and especially by memory effects. 

Different linearization techniques have different sensitivities to 
memory effects. Feedback and feedforward systems (see Section 3.2.2) are 
less sensitive to memory effects because they measure the actual output 
distortion, including the memory effects. However, predictive systems like 
predistortion and envelope elimination and restoration (EER) are 
vulnerable to any changes in the behavior of the amplifier, and memory 
effects may cause severe degradation in the performance of the linearizer. 

However, there is no fundamental reason why predictive linearization 
techniques should be poorer than feedback or feedforward systems since 
the behavior of spectral components, though quite difficult to predict under 
varying signal conditions, is certainly deterministic. Thus, in theory, real 
time adaptation or feedback/feedforward loops are not strictly necessary, 
provided that the behavior of distortion components is known or can be 
controlled. The primary motivation of this book is to develop a power 
amplifier design methodology which yields PA designs that are more easily 
linearized. The approach taken here proposes that, by negating the relevant 
memory effects, the performance of simple linearization techniques that 
otherwise do not give sufficient linearization performance, can be 
significantly improved. 

To achieve a significant linearity improvement by means of simple and 
low power linearization techniques requires detailed understanding of the 
behavior and origins of the relevant distortion components. This is a key 
theme that is carried on throughout this book. The actual linearization 
techniques themselves will not be discussed in detail, but instead, the 
fundamental aim of this book is to give the designer the crucial insights 
required to understand the origins of memory effects, as well as the tools to 
keep memory effects under control. 



1.4 Main Contents of the Book 

Obtaining meaningful data of signal bandwidth-dependent effects has been 
nearly impossible, as most commercially available RF power devices are 
supplied without simulation models, while those that are often fail even to 
fully reproduce the devices’ I-V and Q-V curves. Hence, the predicted 
distortion characteristics from computer simulations is generally regarded 
as unsatisfactory; the results may be accurate within 5 dB, but this is not 
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sufficient for analyzing canceling linearization systems, where subdecibel 
accuracy is a prerequisite. 

In laboratory measurements, the commonly used single-tone amplitude 
and phase distortion (AM-AM and AM-PM) characterization techniques 
actually have a zero bandwidth, and so they completely fail to capture 
bandwidth-dependent phenomena. Therefore, the accuracy of IM3 values 
resulting from AM-AM and AM-PM models suffers when attempting to 
model an amplifier that has memory effects. In addition, the AM-AM 
measurements also suffer from self-heating: The AM-AM measurements 
are performed using continuous wave (CW) signals, resulting in transistor 
junction temperatures quite different from those generated in practice, 
where modulated signals are applied to the PA. 

This book presents several techniques that help understand, simulate, 
measure, and cancel memory effects. The subsequent chapters will provide 
a detailed discussion of the following topics: 



1. A comparison between data available from AM-AM and AM-PM 
versus IM measurements. Normal single-tone AM-AM measurement 
has zero bandwidth, but it can be performed using a two-tone signal 
with variable tone spacing, as well. In this case, the same information 
about the nonlinearity of the device should be available in both the 
fundamental and IM3 tones, but the discussion will show that the large 
fundamental signal masks a considerable amount of fine variations in 
distortion in AM-AM measurements. 

2. To study the phase variations of the IM3 tones, a three-tone 
measurement system will be presented. 

3. Device modeling. Input-output behavioral models can be generated on 
the basis of a completed amplifier, but these do not yield any 
information to aid in design optimization. Instead, the analysis 
presented in this book models the transistor by replacing every 
nonlinear circuit element (input capacitance, g m , and so forth) by the 
parallel combination of a linear circuit element (small-signal 
capacitance, small-signal g m , and so forth) and a nonlinear current 
source. This leads to two important findings: 

a. There are several sources of distortion, and the distortion generated 
in any of these sources can undergo subsequent mixing processes, 
resulting in higher order distortion components than the degree of 
the nonlinearity suggests. 
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b. Distortion is originally generated in form of current, which is 
converted to a voltage by terminal impedances. Thus, the phase and 
amplitude of the distortion components can be strongly influenced 
by the terminal impedances, and especially by the impedances of the 
biasing networks. 

4. Based on the reasoning above, this book includes a review of a 
distortion analysis technique called Volterra analysis, which is based 
on placing polynomial distortion sources in parallel with linear circuit 
elements. The main benefits of this technique are: 

a. The dominant sources of distortion can be pinpointed; 

b. Phase relationships between distortion contributions can be easily 
visualized; 

c. A polynomial model can be accurately fitted to the measured data; 

d. The polynomial models can also be used in harmonic balance 
simulators. 

5. This book also introduces some circuit techniques for reducing 
memory effects in power amplifiers. The standard method of 
minimizing memory effects involves attempting to maintain 
impedances at a constant level over all frequency bands. 
Unfortunately, other design requirements often interfere with this aim 
and cause memory effects. To address this problem, an active 
impedance synthesis technique is introduced, which can be used to 
drive impedances to their optimum values. What is more, this 
technique can be used for electrical and thermal memory effects. 

6. Finally, the book presents a characterization technique for polynomial 
nonlinearities. Since many existing power transistor models are not 
sufficiently accurate in terms of distortion simulations, 
characterization measurements are the only way of obtaining this 
information. This is accomplished using pulsed ^-parameter 
measurements over a range of terminal voltages and temperatures. 



1.5 Outline of the Book 

The main emphasis of this book is on developing a detailed understanding 
of the physics underlying distortion mechanisms, while keeping the 
mathematical formulations in a tractable form. To lay the groundwork for 
the analysis of nonlinear effects in RF power amplifiers, Chapter 2 
discusses certain theoretical aspects related to amplifier circuits. Since RF 
power amplifiers are nonlinear, bandwidth-dependent circuits with 
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memory, it is important to define nonlinearity, bandwidth dependency, and 
memory, and to examine their associated effects. Chapter 2 also introduces 
a direct calculation method for deriving equations for the spectral 
components generated in such circuits. Due to its analytical nature, this 
method, based on the Volterra series, provides detailed information about 
distortion mechanisms in nonlinear systems. Later chapters of this book 
will describe the use of the method. 

Chapter 3 first discusses memory effects from the linearization point of 
view. Some of the most common linearization techniques are presented, 
and then the chapter highlights the harmful memory effects in more detail, 
with a particular focus on electrical and thermal memory effects. Electrical 
memory effects are those caused by varying node impedances within a 
frequency band, while thermal memory effects are caused by dynamic 
variations in chip temperature. Both kinds of memory effects are analyzed 
by comparing a memoryless polynomial model with measurements of real 
power amplifier devices. Memory effects tend to be considered merely in 
terms of modulation frequency, but Chapter 3 also introduces mechanisms 
that produce memory effects as a function of signal amplitude. These 
mechanisms are referred to as amplitude domain memory effects. 

Chapter 4 discusses transistor/amplifier models and introduces 
problems related to PA modeling. The amplifier models are classified as 
either behavioral or device-level models, which are based on some pre- 
defined, physically based functions or simply on empirical fitting functions. 
The Volterra model is an empirical model that is capable of providing 
component-level information that can be used for design optimization. The 
chapter also gives a derivation of the Volterra models for a common-emitter 
(CE) bipolar junction transistor (BJT) amplifier and a common-source (CS) 
metal-semiconductor field effect transistor (MESFET) amplifier. The 
models take into account the effects of modulation frequency, and 
temperature, and are therefore able to model memory effects. Moreover, IM 
products are presented as vector sums of each degree of nonlinearity, 
thereby providing insight into the composition of distortion, which is 
instrumental in design optimization. 

Chapter 5 discusses the characterization of the Volterra model. The dc 
characterization is briefly discussed for the sake of clarity, before shifting 
the focus on a new technique based on a set of small-signal ^-parameters 
measured over a range of bias voltages and temperatures. 

Chapter 6 presents a new simulation technique that offers insight into 
both amplitude and modulation frequency-dependent memory effects. A 
new measurement technique is introduced that allows both the amplitude 
and the phase of the IM3 components to be measured, which is an 
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important improvement over measurements based merely on the 
fundamental signal or amplitude. 

Chapter 7 introduces three techniques for canceling memory effects: 
impedance optimization, envelope filtering, and envelope injection. In 
addition, the chapter presents the source pull test setup for investigating the 
effects of out-of-band impedances. Then, a comparison is presented 
between envelope filtering and envelope injection techniques, and the 
superior compensation properties of the envelope injection technique are 
demonstrated. Finally, a detailed presentation of the envelope injection 
technique is given, and it is shown how both modulation frequency and 
amplitude domain effects can be compensated. A primary advantage of the 
memory effect cancellation approach is that the performance of a 
polynomial predistorter or other simple linearization technique can be 
significantly increased without a substantial increase in dc power 
consumption. Hence, good cancellation performance can be achieved by 
linearization techniques that consume little power, enabling the design of 
linear yet power-efficient PAs. 

Finally, additional supporting information is collected in the 
appendixes. Appendixes A and B discuss the background and limits of the 
Volterra analysis. Appendix C includes a full list of transfer functions, 
describing the path from all of the distortion sources to a given node 
voltage in a common-emitter type single-transistor amplifier. Appendix D 
includes a brief description of some practical aspects of the measurement 
setups and the RF predistorter linearizer used in the measurements 
presented in Chapter 7. 
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Chapter 2 

Some Circuit Theory and Terminology 



This chapter reviews the theoretical background needed for understanding 
nonlinear effects in RF power amplifiers. It begins comfortably by defining 
memory and linearity, and briefly reviewing phasor analysis and the most 
common ways to measure and define the amount of nonlinearity. It is also 
noted that nonlinear effects are more clearly and accurately seen as the 
structure of IM tones than as small AM -AM and AM-PM variations on top 
of the large fundamental signal. Sections 2.2 and 2.3 motivate the use of 
polynomial models, as the calculation of discrete tone spectrums in 
polynomial nonlinearities is easily done by convolving the original two- 
sided spectrums. 

Section 2.4 defines the memory effects as in-band variation of the 
distortion: the behavior of intermodulation distortion at the center of the 
channel is different from that at the edge of the channel. Nonlinear analysis 
methods are very briefly discussed in Section 2.5, and the rest of the 
chapter concentrates on presenting Volterra analysis using what is known 
as the direct method or nonlinear current method. The method is very 
similar to linear noise analysis: Distortion is modeled as excess signal 
sources parallel to linear components. The main advantages of the Volterra 
analysis are that we get per-component information about the structure of 
distortion as well as the phase of these components, so that we can clearly 
see which distortion mechanisms are canceling each other and how to 
change the impedances to improve the cancellation, for example. 

Finally, a simple example circuit is studied to see the analysis 
procedure, and the circuit-level presentation is briefly compared with a 
behavioral input-output model typically used in system simulations. The 
intention is to show that AM-PM can be modeled by an input-output 
polynomial with complex coefficients (or any complex function), but if the 
coefficients are fixed, it cannot predict bandwidth-dependent phenomena. 
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2.1 Classification of Electrical Systems 

Electrical systems can be classified into four main categories as listed in 
Table 2.1: linear and nonlinear systems with or without memory. An 
example of a linear memoryless system is a network consisting of linear 
resistors. Addition of an energy storage element such as a linear 
capacitance causes memory, as a result of which a linear system with 
memory is introduced. 

Nonlinear effects in electrical systems are caused by one or more 
nonlinear elements. A system comprising linear and nonlinear resistors is 
known as a memoryless nonlinear system. Nonlinear systems with memory, 
on the other hand, include at least one nonlinear element and one memory 
introducing element (or a single element introducing both). 



Table 2.1 

Classification of Electrical Systems 





Memoryless 


With Memory 


Linear 


Linear resistance 


Linear capacitance 
Nonlinear capacitance or 


Nonlinear 


Nonlinear resistance 


nonlinear resistance and 
linear capacitance 



2.1.1 Linear Systems and Memory 

Any energy-storing element like a capacitor or a mass with thermal or 
potential energy causes memory to the system. This is seen from the 
voltage equation of a linear capacitance, for example: 

t 

v c «j = £ • j an ■ df (2.i) 



Here, the voltage at time t is proportional to all prior current values, not just 
to the instantaneous value. This is the reason why capacitances and 
inductances are regarded as memory-introducing circuit elements. 

The well-known consequence of memory is that the time responses of 
the circuit are not instantaneous anymore, but will be convolved by the 




Some Circuit Theory and Terminology 



impulse response of the system; in a system with long memory, the 
responses will be spread over a long period of time. This is illustrated in 
Figure 2.1(b) where the time domain output of a linear system of Figure 
2.1(a) with and without memory is shown. Let the input signal be a ramp 
that settles to the normalized value of one. In a linear memoryless system, 
the output waveform is an exact, albeit attenuated (or amplified), copy of 
the input signal. If the system exhibits memory, the output waveform will 
be modified by the energy-storing elements. 




(a) 



(b) 



Figure 2.1 (a) Linear system and (b) its output in a time domain with and without 

memory. 



In the frequency domain, the consequence of memory is seen as a 
frequency-dependent gain and phase shift of the signal. To analyze 
frequency-dependent effects, phasor analysis is commonly used: sinusoidal 
signals are written according to Euler’s equation as a sum of two complex 
exponentials (phasors) 



"it! -Jtu, i 

x = zljCostWjf + <!>!) = — e + — - e , (2.2) 
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where the time-dependent part models the rotating phase that can be frozen 
to a certain point in time (like t= 0), and the complex-valued constant part 
contains both the amplitude A\ and phase information that fully describe 
a sinusoid with fixed frequency CO]. The reader should note that in linear 
systems no new frequencies are generated, and the system is usually 
analyzed using positive frequency +(Dj only. In nonlinear analysis, new 
frequency components are generated, and both positive and negative 
phasors are needed to be able to calculate all of them, as we will see. Also, 
the fact that the complex phasors contain the phase information will turn 
out to be very handy when the cancellation of different distortion 
components is calculated. 

The main advantage of phasor analysis (or using sinusoidal signals 
only, the derivatives and integrals of which are also sinusoids) is that the 
integrals and differentials involved in energy- storing elements reduce to 
multiplications or divisions with /CO, where the imaginary number j means 
in practice a phase shift of +90°. This way differential equations are 
reduced to algebraic equations again, and normal matrix algebra is used to 
quickly solve the circuit equations. Table 2.2 reviews the device equations 
for basic components to be used in phasor analysis. 



Table 2.2 

Impedances and Admittances of Basic Circuit Elements 





Impedance Z = V/I 


Admittance Y = I/V 


L 


j(OL 


1 / (j(OL) =-j 1 (CO L) 


C 


1 / (/'toQ =-j 1 (coC) 


y'coC 


R 


R 


1 / R 



We see that energy- storing elements cause phase shift, while memoryless 
resistive circuits do not. This is further illustrated in Figure 2.2 where the 
impedance Z of a series RC network is shown in a complex plane as a 
vector sum of the impedances of Z#=R and ^C= 1//C0C, calculated at a 
certain value of ft). As Z c is frequency-dependent, the magnitude and the 
phase of total impedance R+l/j(oC vary with frequency ft), which does not 
happen in a memoryless circuit. 

Here, the total impedance of a series circuit was drawn as a vector sum 
of two contributions. Later we will construct the phasors of distortion tones 
as similar vector sums of different contributions. 
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Figure 2.2 Impedance Z of a series connection of R and C shown as a vector sum of 
Z R and Zq. 



2.1.2 Nonlinear Systems 

Next, we discuss the nonlinear effects. A system is considered linear if the 
output quantity is linearly proportional to the input quantity, as shown by 
the dashed line in Figure 2.3. The ratio between the output and the input is 
called the gain of the system, and in accordance with the definition 
presented above, it is not affected by the applied signal amplitude. A 
nonlinear system, in contrast, is a system in which the output is a nonlinear 
function of the input (solid line) (i.e., the gain of the system depends on the 
value of the input signal). If the output quantity is a current, and the input 
quantity a voltage, Figure 2.3 represents a nonlinear conductance. If the 
output quantity is changed to a charge, nonlinear capacitance is presented. 



output 

quantity 

(y) 




Figure 2.3 Linear and nonlinear system. 
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The nonlinearity of a system can be modeled in a number of ways. One 
way that allows easy calculation of spectral components is polynomial 
modeling, used throughout in this book. The output of the system modeled 
with a third-degree polynomial is written as 

y = a x • x + a 2 ■ x 2 + a 3 • x , (2.3) 

where to a 3 are real valued nonlinearity coefficients at this stage of the 
analysis. The first term, a.\, describes the linear small-signal gain, whereas 
the a 2 and a 3 are the gain constants of quadratic (square -law) and cubic 
nonlinearities, introducing the curvature effects shown in Figure 2.3. In this 
chapter, the analysis is limited to third-degree, but up to fifth-degree effects 
will be discussed in Chapter 3 . 

The output of the nonlinear system can be calculated by substituting a 
single-tone sinewave (2.2), shown graphically in Figure 2.4(b), into (2.3). 
In the frequency domain, nonlinearity generates new spectral components 
shown in Figure 2.4(a) and Table 2.3. The output comprises not only the 
fundamental signal (cOj), but also the second harmonic (2a)!) and dc (0) 
generated by a 2 x 2 and the third harmonic (3a)!) generated by « 3 x 3 . This 
spectral regrowth, which will be discussed in more detail later, is not 
possible in linear systems. Figure 2.4(b) shows that, in nonlinear systems, 
the steady-state time domain output waveform is a distorted copy of the 
input waveform. Like spectral regrowth, this phenomenon is not possible in 
linear systems, in which the steady-state output signal is always identical in 
shape to the input (i.e., it can only be attenuated/amplified and/or phase- 
shifted). 



Table 2.3 

Amplitude of Spectral Components Generated by a Single-Tone Test and 
Nonlinearities Up to the Third Degree 



dc 


Fundamental 


2nd Harmonic 


3rd Harmonic 


(a 2 / 2)A 2 


ci ] A+(3<7 3 /4) A 3 


(a 2 / 2)A 2 


(a 3 /4)A 3 



If the nonlinearity coefficients in (2.3) have real values, the system is 
considered nonlinear and memoryless, because the fundamental output 
signal is in phase with the input over the whole frequency range. If the 
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Figure 2.4 Nonlinear effects in frequency and time domains, (a) Input and output 
spectrums and (b) waveforms. 



coefficients include a phase shift (which appears as a complex-valued 
coefficient), a constant, frequency-independent phase shift will exist 
between the input and output signals, thus modeling a nonlinear system 
with memory. Complex-valued coefficients are normally used in 
narrowband behavioral models, as will be shown later. Here it suffices to 
note that memory causes phase shift in nonlinear systems in much the same 
way as in linear systems. 

2.1.3 Common Measures of Nonlinearity 

We now look at the effects of nonlinearity as a function of signal 
amplitude. As noted earlier, new signal components occur at the dc, 
fundamental, second, and third harmonics. The fundamental signal consists 
of the linear term a \ A and the third-order term (3«3/4)A 3 , while the third 
harmonic only comprises the third-order term. The dc and second harmonic 
terms are equal in amplitude and are both caused by the second power term 
(«2/2)A 2 . Figure 2.5 presents the spectral components at the output as a 
function of input signal level, obtained from a polynomial system (2.3) for 
a single-tone sinusoidal input (2.2). As seen from Table 2.3, the second and 
third harmonics increase to the power of two and three of the input 
amplitude. The fundamental signal, however, increases to the power of one 
at low signal levels, but at higher values, the cubic nonlinearity (or any 
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odd-degree nonlinearity in general) starts to modify the linear behavior of 
the fundamental signal. This means that the nonlinearity of the system can 
be considered in two ways: either a generation of new spectral components 
and/or an amplitude-dependent gain of the fundamental signal gain. 

This gives two common measures for nonlinearity: 1-dB compression 
point P ldB where the large-signal gain has dropped 1 dB, and intercept 
points (Pnp 3 ), where the extrapolated linear and distortion products cross. 
By using a third-degree polynomial amplifier model (2.3) with negative a 3 
and single-tone test (2.2) for calculating the compression point and two- 
tone test (2.7) for IIP3, we get the common approximation stating that P ldB 
= P IIP3 - 10 dB and that the IM3 level at the compression point is as high as 
-20 dBc. 

Another widely used measure of nonlinearity is AM-AM and AM-PM 
conversions [1, 2], These figures model the amplitude and phase of the 
fundamental signal with increasing input amplitude. The linear and third- 
order spectral components of a fundamental signal are shown separately in 
Figure 2.6 at a certain amplitude value. Due to the third power dependency 
of the upper vectors, the fundamental signal is increasingly modified as the 
signal amplitude increases. Figure 2.6(a) presents the situation already 
depicted in Figure 2.5. The values of a 1 and a 3 are real and have opposite 
signs, producing amplitude compression at high amplitude values. The 
second plot, Figure 2.6(b), presents the opposite situation in which a\ and 
a 3 are both real and either positive or negative, resulting in AM-AM gain 




Figure 2.5 Illustration of nonlinear effects. The wanted (fundamental) output begins 
to change from its linear 1:1 slope at high amplitude levels and the 
generated spectral components increase as a function of signal 
amplitude. 
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(a) (b) 

Figure 2.6 Amplitude and phase conversions caused by third-order distortion, (a) 
AM-AM compression, (h) AM-AM expansion, (c) AM-PM, and (d) the 
situation shown next in Fig. 2.7. 

expansion. In the third plot, Figure 2.6(c), a { and a 3 display a phase 
difference that deviates from 0° or 180°, thereby producing an AM-PM 
conversion. Note that this combination of AM-AM and AM-PM cannot be 
predicted using a power series with real coefficients, but we need to have a 
complex value for a 3 in the phasor calculations. 

This reasoning can be extended to higher order distortion analysis, as 
well. If, for example the third-order term is in-phase and fifth-order term is 
in an opposite phase with the linear term, we have a response where the 
gain first expands due to cubic nonlinearity and then compresses due to 
fifth-degree nonlinearity, when the signal level is increased. 

We now consider the case shown in Figure 2.6(d), where the magnitude 
and phase of a 3 are 0.1 and 150°, respectively, while the corresponding 
values for aj are 1 and 0°. Figure 2.7 shows AM-PM as a function of 
fundamental gain compression (AM-AM), with a value of approximately 
3.5° at the 1-dB compression point. It must be emphasized here that a 
system operating at 1-dB compression is already heavily nonlinear. 
Linearity requirements are so demanding nowadays that amplifiers are 
backed-off well below the 1-dB compression point, and their AM-PM may 
be as low as 1° or 2° at full power and approach zero with decreasing power. 

The value of AM-PM is very small, so it is a difficult parameter to 
measure accurately. Phase changes in the fundamental signal introduced by 
AM-PM depend on signal amplitude, and very high values are needed to 
make a visible effect. The same observation holds for AM-AM. The 
problem with using amplitude conversions as a figure of merit for 
nonlinearity is that they measure nonlinearity on the basis of the 
fundamental signal, which comprises a strong linear term. Since nonlinear 
effects in the fundamental are small, the measurement of AM-AM and AM- 
PM is highly sensitive to measurement errors. 
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AM-AM [dB] 

Figure 2.7 AM-PM of a polynomial system as a function of AM-AM. From [3], 



Throughout this book, nonlinearity is considered by studying the 
behavior of generated new spectral components. Using the polynomial 
input-output model, the same information about nonlinearity (a 3 ) can be 
seen both from amplitude conversions and the third harmonic component 
(or third-order intermodulation term IM3 in the case of a two-tone test). 
Technically, it is easier and more robust to measure and analyze the 
behavior of distortion tones than AM-AM, in which the nonlinear effects 
appear only as small variations on top of a strong fundamental signal. 



2.2 Calculating Spectrums in Nonlinear Systems 

Integral transforms like Fourier or Laplace transform can be used to 
simplify the analysis of linear systems. With some care, their use can be 
extended to nonlinear or time-varying systems as well. 

It is well known that the time-domain response y(t) of a linear circuit is 
the convolution of the impulse response h{t) and the input signal x(t), as 
shown in (2.4). In the frequency domain this converts to a multiplication of 
the frequency response H(jO)) and the signal spectrum X(jt o). 

y(t) = h(t)*x(l)<* Y(j(o) = H(j( 0 ) ■ X(j( 0 h (2.4) 

where the convolution (*) is calculated with (2.5). A graphical 
interpretation of convolution (used later in Figure 2.8) is that for each value 
of t, we reverse the time axis of x(t), shift it by the amount of t, and then 
integrate the product of h(x) and time-reversed and shifted x(x) over all 
previous values of t, and store the result in place of y(t). 
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h(x)x(t - x)dx 



(2.5) 



For nonlinear systems the convolution operates the other way around: a 
time domain multiplication of two signals corresponds to frequency 
domain convolution of their spectrums. 

y(tl x(t) ■ x(t) O 7(yw) = X(j( 0 )*X(j( 0 ) (2.6) 

Similarly, the spectrum of y(t)=x(t) N is obtained simply by taking an A 1 - fold 
convolution of X(j(o) with itself. It may sound overly academic to calculate 
the spectrum of a nonlinear system as a multiple convolution of the linear 
signal spectrum, but in fact (2.6) is an extremely handy and effective way 
of calculating the line spectrum of a multitone signal numerically (see [4]), 
and either a symbolic or graphical convolution illustrated in Figure 2.8 is a 
rigorous way of obtaining all the possible mixing results falling to a given 
distortion tone. When performed with complex numbers, the convolution 
also preserves the phase information of the tones. 

As an example, the output spectrum of a two-tone test signal in 
quadratic nonlinearity can be calculated as 
given by 

x^: A j • costWj/ + ^j) + A 2 ■ 

V * 1 AV V '* 1 

= e + 

2 2 

A 2 e j j(o 2 t A 2 e ~ j * 2 

+ — - — • e + — - — • 



that is presented in Figure 2.8(a) using a two-sided spectrum. The right- 
hand side of the plot represents the positive frequency axis, and A ] and A 2 
are now complex numbers containing both the amplitudes (Aj/2) and phases 
of lower (cOj) and higher (co 2 ) tones, respectively. Due to odd phase 
response of real systems, the phasors A { and A 2 of the negative frequencies 
on the left are complex conjugates of A l and A 2 . Figure 2.8(b) is identical 
to Figure 2.8(a), whereas Figure 2.8(c) presents the original input spectrum 
with a reversed frequency axis: Positive frequencies are now on the left and 
negative frequencies on the right. Next, the reversed spectrum is slid from 
right to left and compared at all offsets to the original input in Figure 



follows. The two-tone signal is 



cos (to 2 t + (J) 2 ) 

-AV 

e 

-jw 2 t 



(2.7) 
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2.8(a). Figure 2.8(d) presents the situation at a single frequency offset, that 
corresponds to a single frequency in the output spectrum. Now we simply 
multiply all the aligning frequency pairs [shown with dashed line between 
Figure 2.8(a, d)] and place the sum of these products (A 1 A 2 +A 2 A 1 ) as the 
amplitude (actually a phasor) of the generated tone. The frequency offset 
between Figure 2.8(a), (d) corresponds to the envelope frequency f 2 -f\ 
(also called the beat, video, or modulation frequency), but the other tones 
are generated similarly. For example, a frequency offset 2a)! [i.e., the origin 
of the spectrum Figure 2.8(d) aligns with frequency 2aq in the original 
spectrum Figure 2.8(a)] causes the A j phasors in Figure 2.8(a), (d) to align, 
resulting in a second harmonic with amplitude A] 2 in spectrum (e). Finally, 
Figure 2.8(e) presents the complete spectrum generated by squaring the 
two-tone signal in Figure 2.8(a). The procedure demonstrated in Figure 2.8 
is known as spectral convolution. 

Note that it is necessary to use a two-sided spectrum to calculate the 
amplitudes of the distortion tones using the spectral convolution. Hence, all 
amplitudes except the dc term include the term 1/2. 




Figure 2.8 Spectral convolution, (a) The original and (b)-(c) flipped spectrum; (d) 
shows the flipped and shifted spectrum, and (e) is the final convolution 
result. Note that the phasors include the coefficient 1/2. 
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2.3 Memoryless Spectral Regrowth 

This section discusses the spectral regrowth in a memoryless nonlinearity. 
A block presentation of a nonlinear system modeled by an input-output 
polynomial (2.3) is given in Figure 2.9, where the output is the sum of the 
first, second, and third powers y±, y 2 , and y 3 of the input signal, weighted 
by the nonlinearity coefficients a\, a 2 , and a 3 , respectively. In phasor 
analysis, the coefficients may be complex to model the phase shift in the 
nonlinearities. The spectrums in the intermediate points A and B can be 
calculated as a two- and three-fold convolution of the two-sided input 
spectrum, respectively. As an example, the line spectrum of a squared two- 
tone signal in point A is shown in Figure 2.8(e). 

This polynomial system is usually analyzed by assuming that x(t) is a 
nondistorted two-tone signal. In this case, the linear term adjust amplifies 
the fundamental tones at C 0 j and (0 2 (co 2 >ft) i)- The quadratic nonlinearity 
a 2 x 2 rectifies the signal down to dc band to frequencies 0 Hz (dc) and 
co 2 -0J | . It also generates the second harmonic band consisting of tones at 
2to h 2w 2 and (0 1 +C0 2 , called the lower and higher second harmonic and the 
sum frequency, respectively. Similarly, the cubic nonlinearity a 3 x 3 
generates lower and higher IM3 at 2co 1 -co 2 , and 2 <b 2 -C0| and the 
compression/expansion terms (AM-AM) on top of the fundamental tones 
00; and 00 2 , all appearing in the fundamental signal band. It also generates 
the entire third harmonic band consisting of tones at 3o>i, 2(Q 1 +CQ 2 , 
CQj+2co 2 , and 3(0 2 , called the lower third harmonic, the lower and higher 
sum frequencies and the higher third harmonic, respectively. These tones 
are illustrated in the line spectrum shown in Figure 2.10. 

Distortion tones are classified as harmonic (HD) and intermodulation 
(IM) distortion, where the harmonic distortion is simply an integer multiple 
of one of the input tones and IM tones appear at frequencies 
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Figure 2.9 Block presentation of a memoryless system up to the third degree. 
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J ST0 1 +ZCO 2 , (2.8) 

where K and L are positive or negative integers. Another, more practical 
classification is based on the grouping of the tones: in RF applications, the 
dc, fundamental, second and third harmonic bands are far from each other 
and quite easily filtered separately, if needed. However, the IM3 distortion 
appearing in the fundamental band cannot usually be separated from the 
desired linear term. 

The third and the most important classification is based on the order of 
the distortion product, which in short means the number of fundamental 
tones that need to be multiplied to make a distortion product of a given 
order. In a two-tone excitation in Figure 2.10, the fundamental tones COj and 
a) 2 are first-order signals, while dc (0 Hz), envelope ( 02 - 0 )!, second 
harmonics 2a>i and 2co 2 , and the sum frequency aq+a^ are second-order 
signals. These build up the dc and second harmonic bands. Similarly, third- 
order signal components lay in the fundamental (2 cO|-oj 2 , co b co 2 , 2 co 2 - 0 J |) 
and third harmonic bands (3C0J, 2cO|+CD 2 , C0 i+2g> 2 , 3o) 2 ). The amplitudes of 
the Mh-order tones always are proportional to A N , where A is the amplitude 
of the fundamental tone(s). 

Using the notations of (2.8), the order N is sometimes written as N= 
|A)+|L|. However, this rule breaks down when higher order tones fall on top 
of the lower order ones. As an example, look at the fifth-order compression 
term (2.9) below that appears at frequency lCDj+Oct^ but still is of the fifth 
order. 



(A 5 / 32) • e 3 ' 



J<° 2 . P -J a i 



,-/'W 2 



(2.9) 



3 3 3 3 
3 3 



3 3 



3 3 
+ co 



freq. 



Figure 2.10 Spectral regrowth of a two-tone signal. AM- AM is shown as a dashed 
line next to fundamental tones. 
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Then what is the difference between the order of distortion and the 
degree of nonlinearity? So far the input signal has always consisted of first- 
order signals only, and the things have been simple: the first-degree term 
ayx in (2.3) generates first-order tones, the second-degree (quadratic) term 
a 2 x 2 second-order tones, and the third-degree (cubic) term third-order 
tones. However, the case is not so simple any more, if the input signal is 
already distorted, which is the typical case inside a real amplifier. A 
second-degree nonlinearity x 2 essentially makes a product x\x 2 , where the 
x 1 and x 2 are certain input tones. These need not be the same, and their 
order may already be higher than one. For example, multiplying the 
fundamental tone Wj with a second harmonic 2oo 2 inside a second-degree 
nonlinearity generates two third-order tones at 2co 2 -CO| and 2oo 2 +(0 1 . 
Hence, the order of the output tone is the sum of the orders of the input 
tones Xi and x 2 . In one extreme, a purely quadratic (second-degree) 
nonlinearity is capable of generating any order of distortion, if the distorted 
output is always fed back to the input. 

To summarize, the term order is a property of the final distortion 
product, and it is related to the amplitude dependency and frequency of the 
distortion tone. The term degree is a property of the nonlinear device, 
defining the shape of the nonlinearity. The order of the distortion caused by 
an Mh-degree nonlinearity depends both on the degree of the nonlinearity 
and the order of the input signals. In an Mh-degree nonlinearity, N tones 
are multiplied, and the total order is the sum of the orders of these N tones. 

This is illustrated in Table 2.4, where the amplitudes of all the tones 
generated by a third-degree polynomial are shown in a case where the input 
signal is a sum of the fundamental two-tone signal with phasors A i and A 2 
and the second-order distortion tones DC, E, H n , H n , and H 22 at 
frequencies 0, <d 2 -g>i, 2c0|, C 0 |+<b 2 , and 2 <b 2 , respectively.We see that in 
this case, also the second-degree (quadratic) nonlinearity a 2 x 2 can generate 
third-order distortion appearing at the fundamental and third harmonic 
bands. 

Note that contrary to most presentations in textbooks, Table 2.4 also 
contains the phase information and allows the calculation in a case of 
unequal tone amplitudes as well. The table gives the amplitudes and phases 
for a one-sided spectrum (i.e., they are directly the amplitudes of the 
sinusoids), and to make a two-sided spectrum, simply divide all but dc by 2 
and substitute the complex conjugates of the positive phasors to the 
negative frequencies. This table is already quite difficult to build 
analytically, but the encircled terms are easily found by drawing the 
spectrum of the second-order tones and convolving it graphically with a 
two-tone spectrum. 
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Table 2.4 

Spectral components generated in a third-degree polynomial nonlinearity 
y = a\X + a 2 x 2 + « 3 ^ 3 for a sum of two-tone signal phasors A \ and A 2 and 
second-order distortion phasors E, H\ \ , H 2 2 , and II\ 2 at (0 2 — i 2c0j , 2(0 2 , 
and (Oj + oo 2 , respectively. The terms inside the boxes present the third-order 
results generated from first and second-order signals in the input. 



Frequency 


Name 


a\ x 


(a 2 / 2) x 


(u 3 / 4) X 


0 


DC 


DC 


A \A 1 +A 2 A 2 




0 2 -Wi 


envelope 


E 


2A\A 2 




2a)j-co 2 


IM3L 




EAi+Hi 1 A 2 


3d 1 2 A 2 


Wl 


FUNDL 


A \ 


EA 2 +H\ \A] + 
2DCA ] +H l2 A 2 


6A 1 A 2 A 2 +3A 1 A 1 A 1 


oo 2 


FUNDH 


A 2 


EA | +H 22 A 2 + 
2DCA 2 +H n A l 


6 A \A iA 2 + 3 A 2 A 2 A 2 


2 (i> 2 -(i> l 


IM3L 




EA 2 +H 22 A\ 


3A,A 2 2 


2 m l 


2HL 


Hu 


V 




coj + a> 2 


2SUM 


H\2 


2A,A 2 




2 oo 2 


2HH 


h 22 


a 2 * 




3 C0j 


3HL 




H U Ai 


A l 3 


2a)j + co 2 


3SUML 




H \\A 2 +H l2 A l 


3 A, 2 A 2 


2co 2 + co 1 


3SUMH 




H 22 A x +H n A 2 


3 A,A 2 2 


3 co 2 


3HH 




h 22 a 2 
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2.4 Signal Bandwidth Dependent Nonlinear Effects 

Section 2.1 described the classification of electrical systems into linear and 
nonlinear systems with and without memory. This classification is 
presented graphically in Figure 2.11, in which the overlapping segment 
between two areas represents nonlinear systems with memory. This 
segment is further subdivided into two sections. The upper section 
represents a narrowband system, where the transfer function is dependent 
on the center frequency of the system only, while the lower section 
represents a system that is also affected by the bandwidth of the input 
signal. Since all practical systems are more or less affected by signal 
bandwidth, the upper section is referred to as a narrowband approximation 
of a real, bandwidth-dependent system. In this book, bandwidth-dependent 
effects are called memory effects. 

The narrowband single-tone signal used in Section 2.1 is insufficient 
for the characterization of memory effects. Instead, these effects can be 
investigated by applying a two-tone input signal with variable tone spacing. 
The alternative would be to use a real, digitally modulated signal, but it 
would yield less insight in the operation of the analyzed system, as will be 
seen later on. In addition, using a digitally modulated signal for the 
calculation of generated spectral components necessitates a time domain 
analysis tool with a Fourier transformation. The use of a sinusoidal input 
signal circumvents this problem, because spectral components can be 
calculated analytically. 

This book studies the effects of variable tone spacing in detail to 
characterize bandwidth-dependent effects. Applying a two-tone signal to a 
third-degree polynomial system (2.3) results in the following two 



nonlinear systems systems with memory 




Figure 2.11 Definition of memory effects used in this book. From [3], 




26 



Distortion in RF Power Amplifiers 



conclusions concerning IM3 signals at the output: first, they are not 
functions of tone spacing and, second, their amplitude increases exactly to 
the third power of the input amplitude. This is shown by the last column 
and third row in Table 2.4. The equation for the IM3L (lower IM3) 
component is proportional to the power of three while being independent of 
signal bandwidth. However, a comparison between the polynomially 
modeled and actual phases of the IM3L as a function of tone difference in a 
two-tone signal is sketched in Figure 2.12, where large differences can be 
observed between the two. The real phase (and amplitude) of the IM3 may 
deviate at low and high tone spacings (or modulation frequencies), 
indicating the existence of signal bandwidth-dependent nonlinear effects 
with memory, as marked by the lower overlapping area in Figure 2.11. This 
book refers to such effects as memory effects, and distinguishes between 
two distinct types: electrothermal memory effects, which typically appear 
at low modulation frequencies (below 100 kHz), and electrical memory 
effects appearing above MHz modulation frequencies. 




Figure 2.12 Phase of the IM3 component of a system with (solid line) and without 
(dashed line) memory effects. © IEEE 2001 [5], 



The fundamental output of a two-tone input is also modified by a third- 
degree nonlinearity, shown in Figure 2.10 and Table 2.4. As a result, the 
two-tone signals are also affected by the amplitude and phase conversions. 
It then follows that memory effects can be characterized as changes in these 
conversions produced by a varying two-tone input [6]. Unfortunately, a 
two-tone input is hampered by the same drawbacks as a one-tone input. 
Strong linear signals at the fundamental make nonlinear effects difficult to 
measure. This is particularly important in the characterization of memory 
effects, which are usually very weak compared to linear signals. Therefore, 
the analysis of intermodulation components is the most practical starting 
point for the exploration of memory effects. 
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2.5 Analysis of Nonlinear Systems 

Most nonlinear analysis/simulation methods operate either fully or partially 
in the time domain. Standard transient analysis based on numerical solving 
of nonlinear differential equations is an example of the former, and widely 
used harmonic balance method presents the latter. Here the passive 
components are modeled in the frequency domain, but still the responses of 
the nonlinear components are solved in time domain, and outputs and 
excitations are pumped back and forth between time and frequency domain 
using the discrete Fourier transform. Transient analysis can handle any 
form of input signal or even autonomous circuits (oscillators), but it suffers 
from ineffective modeling of distributed components and long-lasting 
initial transients that need to settle before the steady-state spectrum can be 
calculated. In the harmonic balance the signal is necessarily modeled by 
just a few sinusoids, but the initial transient is bypassed and more accurate 
frequency domain models can be used for passive components. An in-depth 
comparison of the basic simulation algorithms can be found in [7], 

The Volterra analysis technique used in this book is calculated entirely 
in the frequency domain, building higher order responses recursively using 
lower order results. Hence, no iteration is needed and it is a very quick and 
RF-oriented analysis method. What is even more important in studying the 
memory effects is that it can separate the sources of distortion exactly in 
the same way engineers are accustomed to doing in noise simulations: The 
dominant contributions can be listed, and the designer can attack them first. 
That kind of information is very valuable for design optimization, but 
usually impossible to derive from transient or harmonic balance 
simulations that usually display only the total amount of distortion. 

In the Volterra analysis, some simplifications and assumptions are 
made, though. The first simplification is that like in harmonic balance, only 
the sinusoidal steady-state response of a single or two-tone excitation is 
calculated. Second, the nonlinearities of the system are modeled 
polynomially (2.3). Using these assumptions, we may apply the Volterra 
method for calculating the output of a nonlinear system, which can give 
either numerical or analytical results for the distortion components. 

The Volterra analysis is reviewed in Section 2.5.1, while Section 2.5.2 
describes the direct or nonlinear current method for calculating nonlinear 
responses. Section 2.5.3, in turn, compares two Volterra modeling methods, 
the first of which provides merely input-output information, whereas the 
other one offers a true insight into the operation of the system. The latter 
method will be used throughout this book for its visualization and 
optimization benefits. More background information can be found in 
Appendix A. 
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2.5.1 Volterra Series Analysis 

Volterra analysis can be considered a nonlinear extension of linear ac 
analysis, and its main difference compared to the often-used power series 
analysis is that it contains also the phase information of the transfer 
functions. It is often calculated symbolically [8-11], in which case the 
transfer functions describing the amplitude and phase of the distortion 
tones as functions of input signals are derived. These transfer functions are 
illustrated in Figure 2.13(a), where H | is the linear (small-signal) transfer 
function, H 2 is the second-order transfer function (producing all the 
second-order tones in a two-tone test), and so forth; the total output y(t) is a 
sum of all these transfer functions applied to the input signal x(t) [8,12,13], 

The difference between linear and Volterra analysis is further 
illustrated in Figure 2.13(b). Linear small-signal ac analysis models the x-y 
input-output characteristic of a circuit element with its first derivative in the 
operating bias point. In Volterra analysis, the actual shape of the I-V or Q- 
V curve is modeled by a best fit, low-degree polynomial function of the 
controlling voltage, and the higher-degree coefficients of the polynomial 
are used to calculate the distortion components. 

The output of the second-order Volterra kernel for the one-tone 
sinewave (2.2) is derived in Appendix A for interested readers. However, 
since the spectral components at the output can be calculated using the 
direct calculation method explained in the next section, an in-depth 
understanding of Volterra kernels is beyond the scope of this book. In a 





Figure 2.13 (a) Schematic representation of a system characterized by a Volterra 

series and (b) comparison between small-signal and Volterra series 
analysis. From [3]. 
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fully numerical form, Volterra analysis has been implemented for example, 
in SPICE [14] and Voltaire XL [15] circuit simulators. 

A word of warning concerning the Volterra analysis is needed. First, 
the polynomial models are notorious for the fact that their response 
explodes outside the fitting range - hence, the model is only locally fitted 
around the desired bias point and applicable over a certain amplitude range. 
The applicable range depends on the fitting range of the polynomial (as 
anything can happen outside the fitting range), the nonlinearity of the 
device, and the degree of the modeling polynomial. As the degree of the 
model always needs to be limited at some rather low degree, some 
truncation error between the actual and the modeled response exists. The 
effects of this truncation error are discussed in some more detail in 
Appendix B. 

Figure 2.13 and the block diagram in Figure 2.9 model only an input- 
output nonlinearity. Due to intentional or nonintentional feedback, most of 
the controlling nodes in real amplifiers also contain distortion components 
and not just the linear contribution. This causes multiple mixing, as already 
pointed out in Table 2.4. For example, the second harmonic may mix with 
the linear term in a second-degree nonlinearity and generate IM3. These 
effects can be taken into account as well, by keeping track of the order of 
the calculated result. If vj contains all linear and v 2 contains all second- 
order voltage phasors and so forth, Vj+v 2 +v 3 +... can now be substituted into 
the polynomial as shown in (2.10). After expanding this we can collect the 
terms of a given order on separate rows, and the last complete row in (2.10) 
shows that the third-order output current i actually consists of three terms: 
input third-order distortion v 3 (if present) multiplied by the linear gain a j, 
input linear signal Vj distorted in the cubic nonlinearity a 3 x 3 , and finally, a 
mixing result of linear and second-order input signals vj and v 2 generated 
in the quadratic nonlinearity a 2 x 2 . Also, higher order terms like v 2 2 or v 3 v 2 
are generated, but they are ignored in this analysis. 



i = u 1 (v 1 + v 2 + ...) + a 2 (v 1 +v 2 + ...) 2 + a 3 ■ (vj + v 2 + ...) 3 
= a lVl 
+ a^v 2 + a 2 v\ 

+ ajV 3 + a 3 v\ + 2a 2 v l v 2 
+ ... 



( 2 . 10 ) 
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2.5.2 Direct Calculation of Nonlinear Responses 

A thorough study of the Vol terra series can be found in [1 1, 13], while this 
book focuses on the frequency domain analysis only. Furthermore, rather 
than determining all third-order products, researchers usually concentrate 
on IM3 responses, which makes it impractical to derive the general form 
A/th-order Volterra kernels. Instead, the direct method [11], also known as 
the nonlinear current method, can be employed to calculate only the desired 
signal components. 

The direct method is based on modeling the nonlinear I-V and Q-V 
characteristics given by the polynomial functions in Figure 2.14 with a 
parallel combination of a linear element and a nonlinear current source, the 
current of which depends on the polynomial coefficients and the controlling 
voltages. This is illustrated in Figure 2.14, where K lG to AT 3G model the I-V 
and K lc to V 3C the Q-V characteristic curves of a nonlinear conductance 
and capacitance, respectively. The linear terms V| G and K lc are modeled 
by a linear conductance and capacitance, and the higher degree 
nonlinearities are modeled by nonlinear voltage-controlled current sources. 
Furthermore, as the Q-V polynomial models the charge, it has to be 
differentiated with time to get the ac current. Note that the control voltage v 
may include also distortion voltages, which results in multiple mixing 
mechanisms, as illustrated in (2.10) and Table 2.4. 




Figure 2.14 Equivalent models for nonlinear (a) conductance and (b) capacitance. 



The procedure for the calculation of the response to a two-tone signal 
can be summarized as follows: 
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1. Evaluate the fundamental (first-order) node voltages using linear ac 
analysis for both tones. 

2. For each nonlinear component, evaluate the second-order distortion 
currents using the fundamental voltage amplitudes. These will appear 
at five sum and difference frequencies: dc, envelope ®i> second 
harmonics 2 a> 2 , 2(0 j, and sum frequency OO 2 +CO]. 

3. Use these distortion currents to calculate the second-order distortion 
voltages in each node using the ac analysis. Note that the distortion 
voltages are deterministic signals and they are summed as vectors, not 
as powers as in noise analysis. 

4. Using the first- and second-order voltages, calculate the third-order 
distortion currents in the nonlinear components. These will appear at 
eight frequencies, two of which are IM3 signals. 

5. Perform the ac analysis again at the frequencies of the third-order 
distortion currents to find the third-order node voltages. 

In short, the linear node voltages are solved first, using small-signal 
analysis. Then, nonlinear analysis is started by modeling the nonlinearities 
by current sources and short (open)-circuiting the linear voltage (current) 
sources. Using linear analysis again, the second-order voltage responses of 
the distortion currents are calculated, and the procedure can be repeated all 
the way to higher order responses. An example of the direct calculation 
method will be given later on in this chapter. 

Since the nonlinearities of the circuit elements are modeled by current 
sources, they will be explained in more detail. Each nonlinearity in the 
circuit is represented by a current source, which is placed in parallel with a 
linearized small signal element. The second-order current sources are 
calculated on the basis of the two-tone test signal, and the values of the 
second-order current source (one-sided) amplitudes at the envelope oq-a^ 
the second harmonic frequencies 2co \ are given in Table 2.5 and those for 
the IM3 results in Table 2.6. Note that Table 2.6 from [1 1] does not give the 
AM-AM term of the fundamental tones or third harmonics, but these can be 
derived using Table 2.4. 

In the tables, K 2G and K 2 c are the second-degree conductive and 
capacitive nonlinearity coefficients and V imn is the voltage of the node i at 
the frequency of m(i>\+n(i> 2 - For example, the third index is zero for 
responses at 2a>i, because 2(0j = 2(Q 1 +0(C>2. Similarly, intermodulation 
responses always include both frequencies, as indicated by the second 
column, and a negative sign corresponds to a negative frequency, which is 
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necessary to make the envelope frequency C0j-(0 2 . If a negative frequency is 
needed, the voltage phasor for it is the complex conjugate of the phasor at 
the positive frequency. 

The conductances are memoryless, but as i=dq/dt, the charge 
polynomial needs to be differentiated with respect to time, and this causes 
the y'co dependency in the nonlinear current of capacitances. Hence, 
capacitors cause very small distortion currents at low frequencies, but high 
currents at the harmonic bands. A two-dimensional conductance is 
controlled by voltages Vj and Vj (e.g., v be and v ce ). Here, both controlling 
voltages need to have one-dimensional polynomials of their own, but there 
are also terms consisting of the cross-products of voltages on both ports. 
These additional cross-terms are listed in the tables. 

Note again that the phase of negative frequency components is opposite 
to positive frequencies (i.e., tj 0 Normal rules of complex 

arithmetic apply, and for example, c 2 is still a complex number with twice 
the phase angle and frequency of c, while |c| 2 = cc is a scalar real number at 
dc. Some care is needed in calculating the responses of IM products 
consisting of both positive and negative frequencies, as terms v m 3 and 
v in 2v in = l v inl 2v in h ave different phase angles and frequencies, even though 
their amplitudes are exactly equal. 

As seen from Table 2.6, the third-order signal components are not just 
functions of cubic nonlinearities, but they are also affected by the second- 
order voltages and quadratic nonlinearities. For example, the distortion 
current i NL generated by a nonlinear conductance at the higher IM3 
frequency 2(j) 2 -Q>i has the amplitude and phase given by 



! nlg(!,-i, 2 ) 3 / 4 - 


^ 3 G 


■ v u 


1,0 


■n 


+ ^ 2 G 


■ V i, 


0, 1 ' 


v i. 


-1, 1 


+ * 2 G 


■ v~ u 


1, 0 ' 


v u 


0, 2 



where the first row shows the effect of cubic nonlinearity and the last two 
rows show the up- and downconversion of the envelope (fj.j j at CO 2 — CO 1 ) 
and second harmonic (fj 02 at 2co 2 ) tones, respectively, and Fj b 0 i s the 
same as V x q given in Table 2.6. This is the same result found in Table 2.4 
and illustrates again that if the controlling voltage of the nonlinear 
component is distorted, also a quadratic nonlinearity (X 2G v 2 in this 
example) can generate third-order distortion - or in a more general way, 
lower degree nonlinearity can also generate higher order distortion. 
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Table 2.5 

Second-Order Currents Caused by Second-Degree Nonlinearities in a Two- 
Tone Test [11] 



Type of 
Nonlinearity 


Nonlinear Current at Frequency 
C 0 1 +©2 or 


Nonlinear Current at 
Frequency 2c0j 


(Trans) 

Conductance 


^2G1 ‘ V i , 1.0 ‘ V i, 0,±1 


-(^io) 2 

2 2G1 ,,1,u 


Capacitor 


yX®! ± C0 2 ) • K 2C 
■r f , i.o -r.-o.ii 


/2oo. 

-r\ c"'-.',.) 2 


Two- 

dimensional 

conductance 

(cross-terms 

only) 


1 /2 ■ K 2G 1 G2 ' r it j o ' V j 0> ±1 
+ 1/2 • ^ 2G1G2 ' V t 0> ±i ' V j i_ o 


1/2 • if 2G1G2 
• ^.i.o- /y. i.o 




Table 2.6 




Third-Order Currents Caused by Second- and Third-Degree Nonlinearities 
in a Two-Tone Test [11] 


Type of 
Nonlinearity 


Nonlinear Current at Frequency 2 GO , ± 0J 2 


(Trans) 

Conductance 


^2G1 ' V i, 1. 0 ' V i, 1. ±1 + ^2G1 ' V U 0, ±1 ' V i, 2, 0 

+ 3/4-K 3G1 -Vl 10 -V it0M 


Capacitor 


(2jO) l ± y® 2 ) • [K 2C ■ Vj j 0 • V j j ±1 
+ ^ 2 C • Vi, o. ±1 • h 2, o + 3/4 -K 3C . Vl J Q • v. 0> ±1 ] 


Two- 

dimensional 

conductance 

(cross-terms 

only) 


1/2 • JST 2G ig2 • [ Vt. o. ±i • Vj , 2 , o + V t> lt o • V h ,. tl 

+ Vj, i, ±i ' V j' i o + V 1, 2 , o ' V j' o_ ±i ] 

+ 1/4- ^202 • [2-^0. ±1-^.1. o-Vj. 1.0 

+ ^1,0-^.0.±i] +1/4 ^3G1G22 

•[ 2 -^i,o-^,o.±i-^,i.o + ^o.±i-^.i,o] 
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2.5.3 Two Volterra Modeling Approaches 

A Volterra representation can be presented either as an input-output or a 
circuit-level description. Having already discussed the input-output model 
described by polynomial nonlinearity coefficients, we now look at the 
component-level model, which takes into account the nonlinearities of each 
circuit element. A comparison between the two models will be illustrated 
using the direct calculation method. Also, the existence of complex 
coefficients in the input-output Volterra model and amplitude conversions 
will be discussed. 

Figure 2.15(a) presents a simple nonlinear circuit with memory. The 
first-order lowpass filter includes a linear series resistor and a nonlinear 
shunt capacitor, introducing nonlinear effects with memory. 




(c) 



Figure 2.15 (a) Linearized first-order lowpass filter, and (b) circuit used for 

calculating nonlinear responses by a component-level Volterra approach, 
(c) The input-output Volterra modeling approach. From [3]. 



First, component-level calculations are applied using the direct 
method. The first-order, linearized transfer function is found to be 



f/lO'CO) = 



1 



1 + jtoCR ' 



( 2 . 12 ) 



The circuit describing the nonlinear responses is shown in Figure 
2.15(b). The linear voltage source is short-circuited and a nonlinear current 
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source is connected in parallel to the nonlinear element, in this case, the 
capacitance C. From Table 2.5, the value of the current source for the dc 
and second harmonic frequencies can be written as 



1/2 • y(cOj (Oj) • K 2 q v out(®i) ' v out( — ® i) 
= 0 



(2.13) 



and 



*nlc(2c°i) — 1^2 • y'2ft)j • K 2C v out( w i) > (2-14) 



where K 1C is the real-valued quadratic nonlinearity coefficient of the 
nonlinear charge. Note that the dc current caused by the nonlinear capacitor 
is zero, and no rectified dc voltage occur at the output. A second harmonic 
voltage exists, however, and is written as 

v out ( 2(0 1 ) = z N[ c (2d),) • 77 , ’(2c0|). i (2.15) 

where 77 7 (2(0 1 ) is the value of transimpedance transfer function from the 
current source to the output node at frequency 2o)i - in this case simply the 
parallel impedance of linear R and C. 

By combining (2.12), (2.14), and (2.15), the value of the second 
harmonic voltage in the output can be written as 



tfftc(2c°i) ' 2’F(2(0 1 ) 



= J ■ ®1 • K 2C 



(2.16) 



2C (I + y'co 1 CR ) 1 + jldt^RC 



It may be not be possible to minimize distortion by affecting the 
nonlinearity of the device (%lc)> but usually the designer has some control 
over the harmonic terminal impedances (as well as in the dc bias 
impedance) and hence 7’F(2co 1 ). This will be exploited later. 

This simple component-level example contains all the information 
available from the system. Let us now try to compress it to a plain input- 
output model with fixed complex coefficients, illustrated in Figure 2.15(c) 
and commonly used in system simulations. Here, the output is written 
directly as a complex function of the input as follows 
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v OUT a \ ' V IN + a 2 ' V ?N + a 3 ' V IN ■ (2-17) 

For a single-tone sine wave, the dc and second harmonic components at the 
output can be taken from Table 2.2 and are both equal to ( 2 . A 
calculation using the component-level Volterra model shows that the dc 
term is zero, while the second harmonic is given by (2.16). This means that 
the input-output Volterra model with fixed coefficients fails to 
simultaneously model both the dc and the second harmonic component. 
Thus, if the scope is restricted to the dc component, we simply define a 2 = 0. 
However, if we are interested in the second harmonic, (2.16) and Table 2.2 
allow us to write 



- 2 • ■ K *c ■ ■ (rr^cj • < 2 - 18 > 

The same result is obtained for the second harmonic using both Volterra 
methods. This permits us to draw two important conclusions regarding the 
differences between the component-level and input-output methods. First, 
the value of the polynomial coefficient a 2 depends on frequency. If the 
frequency changes, its value has to be recalculated. The input-output 
Volterra modeling is accurate only at one frequency, which is why it is 
known as a narrowband approximation of a real system. Second, the value 
of a 2 can be complex to model the phase shift, even though the Q-V 
curvature is modeled by a real coefficient K 2 q. Third, Volterra analysis 
includes phase information (ignored in the memoryless Taylor series), 
which makes it suitable for simulating high frequency effects such as AM- 
PM conversions as studied next using a cubic nonlinearity. 

Tables 2.5 and 2.6 do not show the response at the fundamental, but the 
response for a plain third-degree nonlinearity can be obtained by raising a 
single-tone signal (2.2) to the third power and picking up the fundamental 
terms, which gives a one-sided amplitude of (3/4 )A 2 A (the same we get 
from Table 2.4 with A 2 =0; or from Table 2.6 by setting (02 = ®i ar| d halving 
the amplitudes of the input tones). Here, the controlling voltage across the 
nonlinear C is Vqut ar| d hence 



3/4 • y'cOj • ^ 3C • Vou T ( < b 1 ) • v out( _<b i) • (2-19) 

The third-order compression/expansion of the fundamental output voltage 
can then be calculated as 




Some Circuit Theory and Terminology 



37 



v OUT3( <B l| -*NL3c(®l) ' ^(®l) 

_ , . ( V IN \ 2 V IN R 

r 3C \i +j(0 l CR) ■ 1-7'cOjCT? ' 1 +j(O l RC 



( 2 . 20 ) 



Here, v 0 ut is calculated using (2.12), and note that in the expression for 
the negative frequency v oux (-col), the complex conjugate of phasor v IN is 
used. The last term is again the transfer function from distortion current to 
output voltage, in this case again the parallel impedance of R and C, now 
calculated at ccq. Comparing (2.20) with the input-output modeled 
compression/expansion taken from Table 2.2 (equal to (3a 3 /4)A 3 ), the value 
of a 3 can be expressed as 



"3 * ' *3C • (t+^3c) ' (iTT^cs) 2 ■ (i^y • < 2 - 21 > 

Like a 2 , a 3 of the input-output model is a function of frequency and can be 
a complex value. The nonlinear blocks are often modeled by AM- AM and 
AM-PM conversions, which describe the gain and phase of the fundamental 
signal as a function of input amplitude. By taking the fundamental tone 
(including both the first- and third-order terms) from Table 2.2 and dividing 
the result by the input voltage phasor, the following equation can be written 
for amplitude conversions 

— = a.+3/4-a • |v IN | 2 . (2.22) 

V IN 

The absolute value of (2.22) represents the AM-AM, and its phase 
represents the AM-PM. This is a mathematical formulation of the situation 
presented graphically in Figure 2.6. As illustrated by the figure, the phase 
difference between a\ and a 3 determines the nature of the amplitude 
conversions. Finally, using (2.12) and (2.21) allows the conversions in this 
case study to be written as 



M3UT = 1 

v IN 1+ycOjCT? 
1 



3/4 • /CD, • K 



3C yi + ja x RC, 



l + j(o l CR) U-ycojCfl 



l 



(2.23) 
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Equation (2.23) shows that due to the 90° phase shift of the third-order 
term, AM-PM conversion necessarily appears at high amplitudes. Equation 
(2.23) also shows that, similar to input-output Volterra modeling, models 
based on AM-AM and AM-PM conversions depend on the center frequency 
of the system. Since the conversions provide a narrowband approximation 
of a real bandwidth-dependent system, they fail to take into account the 
memory effects appearing inside the signal band. 

This section has already given an example of the application of the 
direct method to distortion computations, illustrated the existence of 
complex coefficients in the input-output Volterra model, and provided some 
background information on AM-AM and AM-PM conversions. Finally, this 
section will also demonstrate the third-order distortion caused by quadratic 
distortion mechanisms that were neglected in (2.19) to simplify the 
analysis. 

Third-order distortion is produced not only by AT 3C but also by both the 
dc and second harmonic second-order voltages with the following 
mechanism. The controlling voltage of the quadratic nonlinearity contains 
both linear and second-order voltages V | and v 2 , respectively. The quadratic 
nonlinearity K 2 q acts now like a square law mixer, creating a third-order 
product 2vj v 2 as one term in the expansion of (vj+v 2 ) 2 . This is essentially a 
difference tone generated by a quadratic nonlinearity, and we can calculate 
its amplitude using the second row of Table 2.5 by replacing Wj with 2(0j, 
00 2 with -CO], and using the minus sign in the equation: 

^NLC^il = M ' K 2C ' v OUt( 2< ®i) ' v out( _(0 i) • (2.24) 

Since no nonlinear dc component exists in the nonlinear capacitance, (2.24) 
includes only the second harmonic contribution. This shows that the third- 
order distortion is generated also by the cascaded second-order distortion 
mechanisms. Using (2.12) and (2.16), the compression term produced by 
this mechanism is given by 




1 

1-y'cOj CR 



(2.25) 
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2.6 Summary 

This chapter has introduced some important theoretical aspects of electrical 
circuits. Before a system can be analyzed, it has to be classified. The 
classical circuit theory divides circuits into linear or nonlinear systems that 
either exhibit memory or not. If a system is linear, the output is directly 
proportional to the input, no new spectral components can be generated, 
and the steady-state output waveform is identical in shape to the input 
waveform. Nonlinearity, however, means that the output is a nonlinear 
function of the input, so the gain of the system depends on the amplitude of 
the applied signal. Nonlinearity also introduces spectral regrowth and 
modifies the steady-state signal waveform. 

In a memoryless system, the output is an instantaneous function of the 
input. Any change in the input signal occurs instantaneously at the output, 
therefore no phase difference exists between the input and output signals. 
Memory, however, makes the output also a function of previous input 
values. Thus, memory causes delays in transient signals, before the output 
settles to its steady-state value. This is caused by energy storing circuit 
elements. However, it is important to emphasize that memory itself does 
not modify the steady-state signal waveform, rather only introduces a phase 
shift between the input and the output. 

We can look at the nonlinearity of a system in two ways. It can be seen 
as a modification of system gain (and phase) as a function of the applied 
signal amplitude, or as generation of new spectral components. The first 
view describes the fundamental signal that is modified by nonlinear effects. 
The higher the signal amplitude, the more the fundamental signal (its 
amplitude and phase) is modified by nonlinearities. The drawback of this 
classical way of looking at nonlinear effects is that telecommunication 
systems are becoming increasingly linear, making the characterization of 
nonlinearity difficult by using just the fundamental tone. The alternative 
involves the analysis of generated distortion components. These are easier 
to measure, because strong fundamental signals do not cause disturbances 
and also because the distortion components provide more information 
about the analyzed system. 

RF power amplifiers are nonlinear circuits with memory. Nonlinear 
systems are often assumed to be unaffected by the input signal bandwidth, 
but this is not necessarily true for RF power amplifiers, as will be 
demonstrated in later parts of this book. Bandwidth-dependent effects, 
referred to as memory effects in this presentation, will be studied by 
looking at the amplitude and phase of IM3 components as a function of the 
tone spacing of the two-tone input signal. 
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Bandwidth-dependent nonlinear circuits with memory can be analyzed 
using the Volterra analysis. Here, the nonlinear I-V and Q-V characteristics 
are modeled as polynomial functions of the controlling voltages, where the 
first-degree terms are equivalent to the linear small-signal elements and the 
higher-degree terms are modeled by excess current sources parallel to the 
linear elements. The nonlinear response is calculated using the direct 
method, where the linear circuit is solved first to determine the 
fundamental node voltages. These voltages are inputs to the nonlinear 
current sources that model the nonlinearities of the circuit elements. 
Further, the currents transfer to node voltages, which are inputs to higher 
order responses, and so on. The key point here is the transformation from 
nonlinear current to nonlinear voltage, which is determined by the node 
(trans-) impedance at the frequency of the distortion current. Distortion can 
then be minimized not only by the fundamental impedance levels, but also 
by optimizing the out-of-band terminal impedances. 

A Volterra representation can be regarded either as an input-output or a 
component-level description. The first corresponds to a polynomial 
between the input and output quantities, and the second comprises the 
actual schematic of the circuit elements whose nonlinearities are 
characterized by real valued nonlinearity coefficients. These two can be 
compared by making the component-level model first and then extracting 
the input-output model from it. The input-output model is a narrowband 
approximation of a real bandwidth-dependent system, which is accurate 
only at one frequency at a time. If the frequency of interest is changed, the 
coefficients of the polynomial have to be recalculated. This makes the 
input-output model an insufficient tool for the simulation of memory 
effects. Similar to the polynomial input-output model, AM-AM and AM- 
PM curves, widely used as a figure of merit of nonlinearity, are narrowband 
approximations of real bandwidth-dependent systems. Nevertheless, being 
well suited to the study of memory effects, the component-level Volterra 
method is the method of choice in this book. 

One important property of the Volterra method is that the spectral 
components at the output can be expressed analytically, provided that the 
input signal is a sinusoid (or a sum of sinusoids). This is a unique property 
in polynomial modeling, which serves to reduce computational complexity, 
while providing an insight into the nonlinearity mechanisms of the system. 
This is one of the main reasons for the application of the polynomial 
Volterra approach throughout this book. 
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2.7 Key Points to Remember 

1 . Memory is caused by the storage of energy that has to be charged or 
discharged. 

2. The nonlinearity of a system is easier to measure on the basis of 
generated spectra than on variations of the fundamental signal. 

3 . Bandwidth-dependent nonlinear effects are known as memory effects 
and can be analyzed using the component-level Volterra method. 

4. The input-output Volterra method or AM-AM and AM-PM curves do 
not take memory effects into account. 

5. The direct method can be used in distortion computations for 
nonlinear systems characterized by the Volterra series. 

6. Analytical equations for distortion products can be calculated using 
the direct method that models the nonlinearities of the circuit elements 
using nonlinear current sources connected in parallel with a linearized, 
small-signal circuit elements. 

7. Since the amount of device nonlinearities cannot be affected much, 
distortion is most effectively minimized by optimizing the impedances 
seen by the distortion current sources. 
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Chapter 3 

Memory Effects in RF Power Amplifiers 



Memory effects, defined as bandwidth-dependent nonlinear effects, are the 
main topic in this chapter. RF power amplifiers play an important role in 
modern telecommunications, where opposite design goals make the 
performance optimization of amplifiers very difficult. Linearity was 
discussed in Chapter 2, while Section 3.1 investigates the quality of dc-to- 
RF conversion, which involves the primary problem with RF power 
amplifiers; namely, the trade-off between linearity and efficiency. To 
improve the trade-off, the amplifier can be designed to attain good 
efficiency at the expense of linearity. The linearity specification is then met 
by means of some external linearization technique. Unfortunately, memory 
effects cause a variation in intermodulation tones within the signal band. 
Although this may not dramatically decrease the linearity of the amplifier 
itself, it can considerably reduce the performance of the used linearization 
technique, thus deteriorating the trade-off between efficiency and linearity. 

Section 3.1 defines amplifier efficiency, and Section 3.2 reviews the 
most common linearization techniques and discusses the consequences of 
memory effects in them. Next, Sections 3.3, 3.4, and 3.5 concentrate on 
different types of memory effects in RF power amplifiers, such as making a 
difference between electrical memory effects caused by nonconstant 
impedances and thermal memory effects caused by dynamic self-heating. 
Section 3.5 introduces the topic of amplitude-dependent memory effects 
that arise at moderate signal amplitudes. 



3.1 Efficiency 

Efficiency in power amplifiers describes the part of dc power that is 
converted to RF power and can be expressed as follows: 
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where /’out ' s the output RF power and P dc is the power taken from the dc 
source. Power-added efficiency (PAE), however, takes the power of the 
input signal into account and can be expressed by 

pae - Pou C w - , 1 H> (3 - 2) 

where /*[ N is the power of the input signal and G is the gain of the amplifier 
stage. 

The maximum transmitting power level of mobile phones is usually in 
the region of 1W, and the power level of base stations is a great deal higher. 
However, modulators or upconversion mixers are only able to generate 
transmitter signals with a power of below 1 mW. As a result, a substantial 
power gain is needed in the transmitter chain, which uses cascaded stages 
to produce the desired output characteristic. The total efficiency of a two- 
stage cascade is calculated as 



^tot J j - ’ ( 3 - 3 ) 

rii • G 2 r | 2 

where rij and r| 2 are collector/drain efficiencies (not PAEs) of the first and 
second stages and G 2 is the gain of the last stage. We see from (3.3) that the 
total efficiency of the system is dominated by the efficiency of the last 
stage. Let us assume that the efficiency of the last stage is 50% and that it 
has a gain of 15 dB. If the efficiency of the first stage is now changed from 
20% to 30%, the total efficiency increases only by 1% (from 46% to 47%). 

The calculations above prove that significant improvements in total 
efficiency can be obtained by improving the efficiency of the last amplifier 
stage. Consequently, most effort should be put into the trade-off between 
efficiency and linearity at this stage. As the first stages can be designed 
without too much trouble, this book concentrates on the design of the last- 
stage amplifier. Multistage PA design can also entail some additional 
difficulties arising from interstage matching, for example, but since these 
effects are well covered in the literature [1-3], there is no reason to repeat 
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them in this book, even though the test setups and the analysis presented 
here can also be extended to multistage amplifiers. 



3.2 Linearization 

3.2.1 Linearization and Efficiency 

Back-off is the traditional way of meeting linearity requirements in class A 
amplifiers. Once the output power is reduced from its maximum value, both 
the amount of amplitude conversions and distortion products is reduced. 
Unfortunately, the back-off reduces efficiency [4], making it an unattractive 
linearization method for amplifiers. Figure 3.1 presents the level of IM3 as 
a function of efficiency at various back-off values for a class A amplifier. A 
third-degree polynomial input-output amplifier model is used in this 
simulation, and the results indicate that efficiency decreases rapidly when 
lower IM3 levels are required. 

Efficiency and linearity are opposite requirements in traditional power 
amplifier design, and if the goal is to achieve good linearity with 
reasonable efficiency, some linearization technique has to be employed. 
The main idea of linearization is that the power amplifier itself is designed 
to achieve good efficiency at the expense of linearity, after which the 
linearity requirements are fulfilled by external linearization. This is 
demonstrated in Figure 3.1. Let us assume the IM3 specification to be -45 
dBc. Without linearization, the amount of back-off that fulfils the IM3 
specification would result in an efficiency of approximately 10%. The 




Figure 3.1 Linearity of a PA as a function of efficiency in standalone and linearized 
configuration. 
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lower curve presents the linearized IM3 value, with a same linearity 
achieved at the power amplifier efficiency of better than 20%. In this 
example, the power consumption of the power amplifier is reduced by more 
than a factor of two. 

The calculation above considers only the power consumption of the 
amplifier, but in reality linearization also consumes a significant amount of 
power. Let us assume that the output power is compressed by 0.25 dB, 
which is a typical value for a modem telecommunications amplifier. Now 
some 10°- 25/1 ° - 1 — 6% additional power is needed in the output to restore 
the power of the fundamental output signal, and an additional 1% is enough 
for canceling the approximately -25-dBc IM3 components. The total 
additional power needed both to restore the fundamental and to cancel the 
IM3 signals is therefore close to 7% of the output power of the PA, which is 
not excessive. Unfortunately, it is large enough so that the efficiency and 
construction of the linearizer circuitry does matter. 

3.2.2 Linearization Techniques 

Several linearization techniques exist, and they are discussed in more detail 
in [5-7], Only the most common categories are briefly explained here, and a 
short comparison is present in Table 3.1. 

Feedback is commonly used, and it can suppress distortion as long as 
the feedback loop has sufficient incremental gain. To increase the loop 
gain, baseband error amplifiers in Cartesian [Figure 3.2(a)] or polar form 
are commonly used. The necessary up- and downconversions inside the 
loop increase noise sources and loop delays, limiting the stable bandwidth 
of the linearizer typically to below 100 kHz. As the amplifier is driven 
close to compression, also the loop gain and bandwidth vary with signal 
amplitude, complicating the analysis of the system. 

Feedforward [Figure 3.2(b)] is commonly used in wideband amplifiers. 
Here, the distortion generated in the main amplifier is extracted by 
subtracting the linear contribution from the output of the main amplifier. 
This distortion signal is amplified by an auxiliary amplifier and finally 
subtracted from the output. As this arrangement does not contain a 
feedback loop, it has no stability limitations, but still the bandwidth of the 
combiners and phase shifters limits the cancellation bandwidth. 
Unfortunately, the phase shifters and attenuators needed in the feedforward 
loop are quite tricky to tune automatically, and the use of the linear 
auxiliary amplifier noticeably reduces the overall efficiency. 

Predistortion [Figure 3.2(c)] is based on expanding the signal before 
the power amplifier, so that the predistorter-amplifier pair appears as a 
linear circuit. In principle, predistortion is a very power efficient and 
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wideband linearization method, although it typically needs a slow feedback 
to adapt the predistorting function. A simple RF predistorter may consist of 
just a couple of biased diodes, or the predistortion signal can be generated 
already in the digital baseband using adapted lookup tables. 

The last commonly used technique is called either envelope elimination 
and restoration (EER) or a Kahn transmitter. Here, the amplitude 
information is removed from the carrier by limiters and then returned by 
modulating the power supply of the power amplifier. Hence, rail-to-rail 
driving and high efficiency can be achieved. Unfortunately, EER is very 
sensitive to any time or phase difference between the carrier path 
(transmitting phase information only) and the supply modulating path, 
containing amplitude information only. 




Figure 3.2 Linearization principles: (a) Cartesian feedback, (b) RF feedforward, and 
(c) predistortion (PRED means a predistorter). 
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3.2.3 Linearization and Memory Effects 

The complexity of different linearization techniques vary, and so does their 
sensitivity to memory effects. A brief comparison of the presented 
linearization techniques is presented in Table 3.1. 



Table 3.1 

Comparison of Different Linearization Techniques 





Complexity Efficiency 


Band- 

width 


Cancell. 

perform. 


Main cause of 
memory effects 


Cartesian 

feedback 


Moderate 


High 


Narrow 


High 


Loop bandwidth 


Feedforward 


High 


Moderate 


High 


High 


Passive components 


EER 


Moderate 


High 


Moderate 


Low 


Time delays 


RF 

predistortion 


Low 


High 


High 


Low 


Power amplifier 


Digital 

predistortion 


High 


Moderate Moderate Moderate 


PA & BB and IF 
filters 



Feedback systems like Cartesian feedback are quite insensitive to 
memory effects in the power amplifier, as they sample the output distortion 
as it is and try to cancel it with sufficient loop gain. However, to stabilize 
the loop, the bandwidth must be limited, and this reduces the cancellation 
far away from the carrier. The feedforward technique also samples directly 
the output distortion, and then amplifies and subtracts it from the output. 
Here, the dominant memory effects come from the frequency response of 
the auxiliary path and passive components, both of which may reduce 
cancellation far from the carrier. The main concern in feedforward 
amplifiers is the complexity of tuning the two subtracting/summing loops. 

EER technique relies on the matching of two signal paths, one for 
phase and the other for amplitude information. The main concern tends to 
be the delays and linearities of these signal paths. 

Digital and analog RF predistortion are tempting alternatives in the 
sense that the distortion is corrected before the power amplifier; hence, the 
output power of the linearizer circuit is smaller and its efficiency is not 
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such an issue as, for example, in the auxiliary amplifier of a feedforward 
system. However, predistortion systems rely on exact inverse replication of 
the nonlinearity of the power amplifier, which means high sensitivity not 
only to memory effects but to drifting as well. Typically, some form of slow 
adaptation is needed for the predistorter. Digital predistorter is naturally 
more flexible, but it requires a high bandwidth and dynamic range from the 
digital baseband, and also all baseband and intermediate frequency (IF) 
filters between the predistorter and the power amplifier contribute to the 
memory effects (see [8]). 

Much of the discussion in this book is related to implementing analog 
RF predistortion systems, or reducing the memory effects of the amplifier 
to such a low level that a simple memoryless digital predistorter can be 
used. Still, the analysis techniques presented are applicable to most of the 
other linearization techniques as well. 

A simplistic way to look at any linearizer is to consider it a canceler: a 
certain amount of distortion is generated and it must be canceled with 
exactly the opposite phase replica of that distortion. Good cancellation 
performance places very tight requirements on the amplitude and phase 
match between the distortion components of the amplifier and the signal 
components generated in the linearizer. This cancellation is demonstrated 
in Figure 3.3(a). The power of the residual IM component can be calculated 
using the cosine rule, and the required matching for a given degree of 
cancellation is shown in (3.4), where A<|) and A4 are the phase and 
amplitude errors, respectively. Figure 3.3(b) shows the corresponding 
numerical values. To achieve a 25-dB reduction in the IM level, for 
example, the phase error cannot exceed 2° to 3° and a gain matching A 4/A 
(flatness) better than 0.25 dB (3% error in amplitude) is needed over the 
entire signal and IM band [9, 10]. 

CANC = 10 • log ( 1 - 2(1 + AA/A)cos(A<|)) + (1 + AA/A) 2 ) (3.4) 

Figure 3.3 clearly illustrates the importance of memory effects. In an 
amplifier with memory effects, the amplitude and phase of the distortion 
components vary as functions of modulation frequency (the tone spacing in 
a two- tone signal) and amplitude. Cancellation signals must track the actual 
IM signals of the amplifier over the modulation bandwidth, and tracking 
errors at any modulation amplitude or frequency will cause a deterioration 
in cancellation performance. This is illustrated in Figure 3.4. The relative 
phase of the IM signals varies with the modulation frequency, but in simple 
analog predistorters, for example, predistortion signals are usually 
memoryless with fixed phase. This problem is often avoided by using a 




50 



Distortion in RF Power Amplifiers 




Figure 3.3 (a) Principle of distortion cancellation and (b) the achieved cancellation 

as a function of phase and amplitude error. From [9], 



more complicated digital predistortion algorithm, a feedforward amplifier - 
or by a power amplifier with a low amount of memory effects. 

Distortion components are deterministic signals that vary with the 
instantaneous amplitude and modulation frequency of the signal. 
Nevertheless, they always behave similarly under similar conditions. The 
main contribution of this book is in finding out how distortion components 
behave under varying signal conditions. This can aid to improve the 
amplifiers so that good cancellation is achieved using simple RF 
predistorter type linearization techniques that normally do not provide 
enough cancellation. By carefully studying their distortion behavior, 
cancellation can be improved up to 20 to 30 dB that corresponds to the 
cancellation performance of more sophisticated linearization techniques. If 
this is achieved, more simple and low-power linearization techniques can 
be used. 



L_ 

\ 



two-tone 
signals IM3H 

/N /N <7^ 




Figure 3.4 Principle of distortion cancellation and its sensitivity to memory effects. 
©IEEE 2001 [11]. 
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3.3 Electrical Memory Effects 

To determine the mechanisms of memory effects, it is important to 
investigate why the real power amplifier device differs from the polynomial 
input-output model. This section first presents the impedance definitions of 
transistor amplifiers. Then, distortion composition is analyzed in more 
detail and compared with the single polynomial model. Finally, the effects 
of matching impedances are discussed in terms of memory effects. 

The notations for the common emitter/source BJT and MESFET 
amplifiers given in Figure 3.5 are used throughout this book. Z B (match) is 
the driving impedance of the stage, from which the base bias impedance 
Z B (bias) is excluded. These two correspond to the impedance that is 
measured by a network analyzer (NWA) upon disconnecting the transistor. 
Z B (int) is the bias-dependent internal base impedance. Similarly, the 
external collector impedance consists of a load impedance Z L and a 





Figure 3.5 Definition of impedances in (a) a CE BJT amplifier and (b) a CS 
MESFET amplifier. From [12]. 
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collector bias impedance Z c (bias), which are both measurable from the 
collector node. Z c (int) is the internal collector impedance. However, node 
impedance refers to the impedance level of the node, and the impedance of 
the base and collector nodes can be calculated by 

Z bb = Z g (match) II Z g (bias) II Z g (int) (3.5) 



and 



Z CC = z l 11 Z c (bias ) II z c (int) . (3.6) 

In the same way, changing the names of the terminals allows us to find 
equations for the node impedances of the MESFET presented in Figure 
3.5(b). Instead of source impedance, the term input impedance will be used 
for both the BJT and the MESFET to make a clearer distinction between it 
and the source terminal of the MESFET. Equations (3.5) and (3.6) describe 
the node impedance outside the transistor, although the impedance seen by 
the internal distortion generator is of primary interest here. These internal 
impedances will be calculated and used in the simulations and analyses 
presented in later chapters of this book. 

Real power amplifier devices contain more than one nonlinearity 
mechanism. As a result of their interaction, nonlinear responses are not just 
output signals as assumed in (2.3); rather, they act as inputs to other 
nonlinearities and are thus capable of generating new nonlinear responses. 
Thus, to improve our understanding of distortion mechanisms, the 
following simplified presentation regards the transistor amplifier as a 
cascade of two nonlinearities. Although this model lacks the feedback 
effects of real amplifiers, it is informative and provides an insight into the 
composition of distortion. 

A cascade of two connected Volterra kernels can be presented in the 
form of a block diagram given in Figure 3.6. Block H describes the base 
voltage as a function of the input signal, and block F the collector voltage 
as a function of the base voltage. Hj, H 2 , and H 3 correspond to the different 
order blocks represented by the coefficients ay, a 2 , and a 3 in (2.3), because 
polynomial input-output models reduce kernels to polynomial coefficients. 
The generation of IM3 by means of third-degree nonlinearities is 
straightforward [13], First, the third-order block H 3 of the base nonlinearity 
generates an IM3 signal at the base, which is linearly amplified in F | and, 
second, the linear signal at the base goes to the cubic nonlinearity of the 
transconductance F 3 , also producing IM3. The generation of IM3 by 
cascaded second-degree nonlinearities is somewhat more complicated. 
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Figure 3.6 Representation of the cascade connection with Volterra operators of the 
order of one to three. © IEEE 2000 [14], 



First, an envelope component is generated at the base in H 2 which, in turn, 
creates the IM3 component together with the linear signal at the base Hj in 
the quadratic nonlinearity of the transconductance F 2 . Similarly, mixing 
from the second harmonic adds to IM3 . 

A frequency domain combination of the different order spectrums is 
given in Figure 3.7. Figure 3.7(a) presents the output of the first block that 
includes the same frequency components as the third-degree polynomial 
model (2.3). The amplitude of the spectral components can be found from 
Table 2.4. This multitone signal is the input signal for the second block, 
and the output IM3 now combines with other frequency components. The 
envelope signal oo 2 — (Oj and the upper two-tone signal (0 2 , for example, will 
be mixed in the quadratic nonlinearity of the latter block, which results in 
the generation of the upper IM3 signal (and COj compression, as well). 
Similarly, the second harmonic of the upper input signal 2(0 2 and the lower 
input signal from the negative frequency side — will also mix to the IM3 
signal. As a result, IM3 sidebands are affected not only by the fundamental 
voltage waveforms, but also by the voltage waveforms of the different 
nodes at the envelope and second harmonic frequencies oo 2 — COj and 2oo 2 . 

The question is how to control the voltage waveforms of the different 
nodes and frequency components. Since the nonlinearities of the circuit 
components can be regarded as current sources, as explained in Chapter 2, 
their voltage waveforms can be affected by node impedances. The 
composition of IM3 in the real power amplifier device is sketched in Figure 







54 



Distortion in RF Power Amplifiers 



VA 



3 3 3 3 
3 3 



3 3 3 

r-4 + fS 



3 3 3 3 

cn + + co 



freq. 



(a) 




Figure 3.7 (a) Spectral components produced by nonlinearities up to the third order 

and (b) composition of IM3. © IEEE 2001 [11], 



3.7(b) with nonlinearities up to the third degree [15]. IM3 is largely 
generated by cubic nonlinearities, which are affected by fundamental 
impedances and signal levels. However, quadratic mechanisms that mix the 
envelope and second harmonic frequencies with the fundamental tones also 
have a significant effect on IM3 distortion, and these can be controlled by 
node impedances at these out-of-band frequencies. 

Electrical memory effects are caused by frequency-dependent 
envelope, fundamental or second harmonic node impedances. Figure 3.8 
gives the measured gate node impedances for the MESFET amplifier in the 
dc, fundamental, and second harmonic bands. The center and maximum 
modulation frequencies are 1.8 GHz and 20 MHz, respectively, which 
means that the dc band is important up to 20 MHz or beyond. The 
interesting fundamental band is between 1.77 GHz and 1.83 GHz, because 
the entire IM3 band of 60 MHz is relevant in terms of IM3 distortion. The 
second harmonic band lies between 3.58 GHz and 3.62 GHz. The 
fundamental impedance can easily be kept constant over the entire 
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modulation frequency range, because it is just 0.3% of the center frequency 
in our example. Also, the second harmonic band is quite narrow, and 
impedance matching is simple, provided that no harmonic traps are used. 
Such traps cause tremendous impedance variations and may cause 
significant memory effects. As the fundamental and second harmonic 
impedances play a minor role, memory effects are for the most part 
produced by envelope impedances. The envelope frequency varies from dc 
to 20 MHz, and the gate node impedance, for example, must be constant or 
very low over this region to eliminate memory effects. This is not the case 
in the practical implementation presented in Figure 3.8, however, where the 
gate impedance at the envelope frequency varies by approximately two 
decades. 

There is one important difference between the dc and the other 
frequency bands. If the center frequency of the system changes, both the 
fundamental and the second harmonic impedance change, while the 
envelope impedance remains the same. In other words, if a significant 
amount of memory is produced by the fundamental or second harmonic 
bands, the shape of these memory effects will change with the frequency 
channel. Despite that, it can be concluded that with careful design, memory 
effects introduced by various terminal impedances can be limited to those 
converted from the envelope frequency. A thorough analysis of distortion 
mechanisms will be presented in Chapter 4, where the distortion 
mechanisms and memory effects of the BJT and MESFET amplifiers are 
analyzed in detail. 




Figure 3.8 Measured magnitude of the Z GG of the MESFET amplifier. From [12], 
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3.4 Electrothermal Memory Effects 

Electrothermal memory effects are caused by electrothermal couplings, 
which affect low modulation frequencies up to the megahertz range. The 
dissipated power of the BJT can be expressed as 

^dissW = v ce(*) ‘ (3-7) 

where v CE is the collector-emitter voltage and i CE the collector-emitter 
current. Since two first-order fundamental signals are multiplied together, 
the spectrum of the dissipated power always includes second-order signal 
components (i.e., dc, envelope, sum, and second harmonics [16]). The 
temperature variations caused by the dissipated power are determined by 
the thermal impedance (Z TH ), which describes the ratio between 
temperature rise and heat flow from the device. Due to the nonzero mass of 
the component, thermal impedance in the active device is not purely 
resistive, but forms a distributed lowpass filter with a wide range of time 
constants. 

This means that the temperature changes caused by the dissipated 
power do not occur instantaneously, but due to the mass of the 
semiconductor and the package, a frequency-dependent phase shift always 
exists. Moreover, the surface of the silicon reacts surprisingly quickly, and 
thermal effects can be obtained in bandwidths up to 100 kHz to 1 MHz [16- 
20], Furthermore, since heat in the chip flows mostly vertically [21], it can 
be assumed that self-heating within the component produces more memory 
effects than the heat generated by surrounding heat sources. 




Figure 3.9 Heat flow from the device: (a) physical and (b) electrical lumped element 
models. ©IEEE 2001 [11], 
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A finite element model is employed here to simulate the thermal 
impedances of the package presented in Figure 3.9. For the sake of 
simplicity, the structure of the package presented in Figure 3.10 is modeled 
using brick elements. The silicon chip measures 600 pm by 600 pm by 300 
pm and the emitter of the transistor 400 pm by 400 pm. The thicknesses of 
the copper lead frame and the moulding compound are 100 pm and 1000 
pm, respectively, and the temperature of the bottom of the moulding 
compound is assumed to be constant. The structure includes 1859 nodes, 
that is, the grid is 50 pm in length and width and the thickness of the grid 
of the chip, the lead frame, and the molding compound are 50 pm, 50 pm, 
and 500 pm, respectively. 

Figure 3.11 illustrates three simulated thermal impedances for the 
entire structure. The first one is simulated at the center of the active area 
(marked with a square), the second at the corner of the chip (diamond), and 
the third a fitted lumped model with three time constants (cross). The figure 
indicates that the surface of silicon reacts quickly, and a thermal impedance 
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Figure 3.10 Simplified package structure. From [22], 





(b) 



Figure 3.11 Simulated thermal impedance at different locations of the integrated 
circuit (IC): (a) on a two-sided linear frequency axis, and (b) on 
logarithmic frequency axis. From [22], 
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of several ohms (approximately 10% of /? th ) can be seen up to frequencies 
of 1 MHz. To illustrate the odd phase response, Figure 3. 1 1(a) is plotted on 
two-sided frequency axis, while in Figure 3.11(b), typical log-log 
magnitude and log-lin phase plots are shown. A more detailed analysis of 
this thermal impedance can be found in [23-25], Thus it suffices to say here 
that the effects of the package and the heat sink are important from the 
point of view of thermal resistance, which determines the average 
temperature rise caused by self-heating. From the ac behavior point of 
view, layers near the active area (silicon and lead frame) are more 
dominant, because the responses of package and heatsink are usually far 
too slow compared to microsecond range variations in power dissipation. 
Thermal impedance simulations for the GaAs MESFET are presented in 
[26, 27], where similar time constants are reported for silicon to those that 
were found in the simulations presented here. 

Since power is dissipated at dc, fundamental, and second-order signal 
frequencies, but only dc and envelope components of the dissipated power 
fit into the passband of the thermal filter, the temperature of the chip takes 
the following simple form 

T = r AMB + ^TH ' ^DISS (0 ^ z) + Z TH (co l " ®2 ) ' ^DISS^l " ®2 ) ' 

The temperature of the chip consists of three components: one is 
simply the ambient temperature ^AMB and the other two consist of the 
thermal resistance multiplied by the dc power dissipation, and the envelope 
component multiplied by the thermal impedance at that frequency. It is 
interesting to note that the third term in (3.8) includes frequency, which 
means that the temperature variations at the surface of the chip also depend 
on the bandwidth of the signal. If any of the electrical parameters of the 
transistor are affected by temperature, thermal memory effects are 
unavoidable. This mechanism in which dynamic self-heating causes 
electrical distortion is known as thermal power feedback (TPF) [28] . 

A block diagram of TPF is shown in Figure 3.12, in which the basic 
amplifier is considered to be a polynomial input-output stage. Thermal 
impedance describes the relationship between dissipated power and 
temperature, and block K describes the relationship between temperature 
and the gain of the amplifier. Only the gain of the amplifier is considered to 
be temperature-dependent in this behavioral model. In practice, however, 
the output conductance [29] and the capacitances are also temperature- 
dependent at the transistor level, as will be seen in Chapter 4. Since some 
of the circuit parameters of the transistor are always functions of 
temperature, TPF cannot be avoided. TPF is a very difficult problem to 
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Figure 3.12 Block diagram of thermal power feedback (TPF). © IEEE 2001 [11], 



compensate for, because the exact chip temperature cannot usually be 
measured. For example, temperature-compensated external bias networks 
are incapable of detecting junction temperatures, and hence far too slow to 
compensate for changes therein and cannot offer an improvement in 
thermally induced distortion. 

We now demonstrate thermal power feedback. We choose 0.1 5-y'0.1 5 as 
the normalized cubic nonlinearity coefficient a value which 

corresponds to an IM3 level of -40 dBc at the signal level employed. The 
dG/dT is — 0.6%/K and the thermal impedance is taken from Figure 3.11. A 
two-tone input signal is modulated by thermally induced gain variations at 
the envelope frequency, thereby generating IM3 sidebands. Since the phase 
response of the thermal filter at the positive envelope (generating IM3H) is 
opposite to that at the negative envelope (IM3L), the IM3 sidebands 
produced by TPF turn in opposite directions as a function of modulation 
frequency, as shown in Figure 3.13. A linearity decrease of several decibels 
is observed at low modulation frequencies, and a phase imbalance exists 
between the sidebands at some relatively low values. 



3.5 Amplitude Domain Effects 

Memory effects that are introduced into power amplifiers affect their 
distortion performance related to both modulation frequency and 
amplitude. Modulation frequency-dependent effects were used in the 
previous sections to demonstrate the mechanisms of memory effects, while 
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Figure 3.13 IM3 caused by the basic amplifier and TPF. (a) Represents the magnitude 
of IM3 as a function of modulation frequency and (b) presents the IM3 in 
real-imaginary coordinates. From [22], 



this section examines the effects of amplitude on memory effects. It is 
important to note that, due to the definition, the memory effects considered 
here are actually both modulation amplitude and frequency dependent. The 
term “amplitude dependent” is justified, because effects higher than third- 
order are taken into account here, and the amount of fifth-order distortion 
in the IM3 tone depends on signal amplitude. These effects are slightly 
more complicated than frequency domain memory effects, wherefore third- 
and fifth-order distortion composition will be studied first without memory 
effects. They will be taken into account later by considering the amplifier 
as a cascade of two polynomials with a bandwidth-limited connection. This 
gives a useful insight into distortion mechanisms for the further simulations 
to be presented in Chapter 6. 

3.5.1 Fifth-Order Analysis Without Memory Effects 

Polynomial input-output relations up to the fifth degree can be written as 

, 2 3 4 , 5 

y = a { ■ x + a 2 ■ x + a 3 • x + a 4 ■ x + a 5 ■ x , 



(3.9) 
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where aj-a 5 are real valued coefficients. By applying a two-tone signal 
(2.7), the in-band intermodulation products shown in Figure 3.14(a) can be 
obtained. 

The relation between the degree of the nonlinearity (third, fifth,...) and 
the frequency of the tone (such as IM3, IM5,...) is demonstrated in Figure 
3.14(a). The IM5 tones are not affected by third-degree nonlinearities, but 
IM3 tones are functions of both third- and fifth-degree nonlinearities. This 
means that at low signal amplitudes, where the fifth-order distortion 
products can be neglected, the amplitudes of the IM3 tones are proportional 
to the third power of the input amplitude. With a fairly large signal 
amplitude, however, fifth-order products (which are dependent on a power 
of five) will start to affect the IM3 responses. As a result, the 3:1 amplitude 
estimate will no longer hold, as demonstrated in Figure 3.14(b). If the 
phases of the third- and fifth-degree coefficients are equal, the fifth-degree 
nonlinearity will expand the IM3 responses. However, if the phases are the 
opposite, the IM3 distortion will be locally reduced, as shown in the figure. 
This explains why notches in the IM3 sidebands have been reported at 
certain amplitudes [30, 31]. It is also interesting to note that the amplitude 
of IM3 (25/8) caused by the fifth-degree nonlinearity is five times greater 
than that of IM5 (5/8). This information is necessary for the identification 
of amplitude domain memory effects. 



a 3 A 3 

a 5 A5 

(a) 



3/4 



9/4 

4 ^ 



3/4 



IM3L FUND IM3H 
25/4 

5/8 2 1 8 ?? 2 ^ 5,8 
m/l IM'/L FUND IM3H «5H 




Figure 3.14 (a) Distortion components caused by third- and fifth-degree 

nonlinearities and (b) amplitude of IM3 and IM5 components as function 
of input amplitude. From [32], 
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3.5.2 Fifth-Order Analysis with Memory Effects 

The cascade representation of nonlinear systems introduced in Section 3.3 
will now be extended to fifth-degree nonlinearities. As noted in the 
previous section, IM3 is also generated by fifth-degree nonlinearities. The 
connection between the two blocks is bandwidth-limited to provide an 
insight into amplitude-dependent memory effects. A block diagram is given 
in Figure 3.15(a). 




7th-order result I M3^-^0=3ffl,-g I q 

”^|^_(2® 2 ±®iH2®2-®it 2 ^i 3/16*b 3 *a 3 *A 5 *a 1 2 

5th-order results j i^_(3m 2 H 2 ®2"®^)^®i+®2_ _^J2/16*b 3 *a 3 *A^*a ^ 

1 W (3 “2H 3 “2-2®II=2®i_ Vl/16*b3*a3*A 5 * ai l2 



Figure 3.15 (a) Cascade of two polynomials, and (b) mixing of distortion of block A 

in block B. From [32], 



The filter H(j(o) will be ignored at first, and since a cascade of two 
polynomials is still just a polynomial; no difference exists between it and 
the single polynomial presented in Section 3.5.1. The coefficients are 
different, but the 5:1 ratio between the fifth-order IM3 and IM5 terms 
remains constant. Four frequency bands have an impact on the IM 
responses: dc, fundamental, second harmonic, and third harmonic. Signals 
occur at the fourth and fifth harmonics as well, but since the distortion 
analysis is truncated at the fifth-order, these signal components do not 
affect the in-band intermodulation distortion. Each band consists of three 
or more individual tones. Figure 3.15(b) presents the significant signal 
components of the IM3 and IM5 distortion. 
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Next, the effects of the third harmonic to the in-band intermodulation 
distortion will be studied. The third harmonic band consists of frequency 
components at six frequencies, of which the outermost (4cO|-co 2 and 4oo 2 - 
Wj) are fifth or higher order distortion products. This is obvious, because 
these frequencies cannot be combined from three fundamental (Wj or oo 2 ) 
frequencies only. However, the middle four spectral components consist of 
both third- and fifth-order terms, of which the third-order ones are written 
in Table 2.5. We now neglect all fifth-order effects at the third harmonic, 
because mixing them to the fundamental band corresponds to at least a 
seventh-order effect. At this stage, we have to mix down to the 
intermodulation band the spectral components at 3c0j, 2oo 1 +co 2 , 2oo 2 +CQ], 
and 3(0 2 with amplitudes of a 3 /4, 3a 3 /4, 3u 3 /4, and a 3 /4. Furthermore, it is 
assumed that the role of the second harmonic in the middle of the blocks is 
negligible (i.e., the second and third harmonics do not mix back to the 
intermodulation band in the latter block B). The intermodulation distortion 
at the output from the third harmonic is now composed of the following 
phenomena: the third harmonic band generated by a 3 and the fundamental 
components (taken two times) produced by aj, mixed in the cubic 
nonlinearity of the latter block B. This is presented in Figure 3.16. 

This particular distortion mechanism produces the IM3 and IM5 
components of 5/16*a 1 2 *a 3 *h 3 *A 5 and l/16*aj 2 *a 3 *b 3 *A 5 that contribute 
the IM3 and IM5 in the same amplitude ratio of five. In other words, if the 
third harmonic is filtered out with filter H, the amount of fifth-order IM3 



block B 



< 

o 




a) (a 3 /4)A 3 

b) (3a 3 /4)A 3 




fundamental band 



c) 1/16*b 3 *a 3 *A 5 *a 1 2 

d) 5/16*b 3 *a 3 *A 5 *ai 2 
= 3/16*b 3 *a 3 *A 5 *a 1 2 
+2/16*b 3 *a 3 *A i *a 1 2 



Figure 3.16 The mechanism of IM3 and IM5 downconverting from the third 
harmonic band. 
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and IM5 changes, but the ratio of five obtained from the memoryless 
distortion composition remains constant. This observation can be 
generalized as follows: the same 5:1 ratio is also true for other distortion 
mechanisms that are converted from the dc and the second harmonic bands. 
Once the frequency band is filtered out, the ratio remains the same. 
Moreover, when the filtering inside the frequency band is flat (i.e., all 
spectral components are attenuated identically), the 5:1 ratio exists 
between the IM3 and IM5 responses caused by the fifth-degree 
nonlinearity. It is still important to emphasize that the ratio between the 
total IM3 and IM5 changes, because IM3 also comprises a third-order 
component that is not a function of third harmonic filtering. 

Let us now take a look at filtering inside the frequency band. Figure 
3.15 sketches the important signal components at the third harmonic band 
in terms of intermodulation distortion. The upper third harmonic 3co 2 mixes 
down to IM3 and IM5 at the amplitudes of 2/16 and 1/16, but the sum 
frequency 2 cd 2 +C0] only mixes to IM3 at the amplitude of 3/16. There is no 
mixing down to IM5, because the difference between the 2co 2 +co 1 and the 
3(0 2 -2(0i tones is 3co 1 -0 2 . Since the 2co 1 +(0 1 and 3cflj-co 2 are at least third- 
and fourth-order effects, IM5 converting down from the 2co 1 +(0 1 is at least 
a seventh-order effect As the spectral components at the third harmonic 
band are divided into its contributors, it becomes clear that if the sum 
frequencies 2 oo 1 +cq 2 and 2 cd 2 +co 1 exist, but the harmonics (3a>i and 3(0 2 ) 
are filtered, the value of the ratio between IM3 and IM5 caused by fifth- 
degree nonlinearities at the third harmonic band is infinity instead of five. 

Amplitude -dependent memory effects can be understood as deviations 
from the 5:1 amplitude ratio for fifth-order IM3 and IM5 responses or as 
phase differences between the two. This occurs whenever one of the 
frequency bands is filtered unevenly. To study these effects, the filter H(j(o) 
is considered to be tilted around the dc band, because it is the most 
important source of memory effects in practical applications. The 
composition of fifth-order distortion mechanisms to IM3 and IM5 is 
studied as a function of tone spacing. Since the frequency of the fourth- 
order envelope 2co 1 -2(0 2 is twice that of the second-order envelope (0j-(0 2 , 
the IM3 and IM5 generated by the fourth-order envelope will include a 
significant amount of memory that is a function of signal amplitude. 

This means that the 5:1 ratio fails to hold whenever the filter is tilted, 
and if the phase of the filter is not equal at the second and fourth-order 
envelopes, a phase imbalance will exist between the fifth-order 
contributors, as seen in Figure 3.17(a). A similar conclusion can be reached 
concerning the other frequency bands. In other words, the spectral 
components that are furthest away from the center of the frequency band 
are the most harmful in terms of amplitude-dependent memory effects. 
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Figure 3.17 (a) Composition of fifth-order distortion and (b) memory effects caused 

by it. From [32], 



The ratio between the fifth-order IM3 and IM5 contributors as a 
function of modulation frequency is presented in Figure 3.17(b). If the 
modulation frequency is close to zero, no memory effects occur and the 
amplitude ratio is 5:1. The ratio starts to increase with increasing 
modulation frequency, however, and acquires the value of 6 at 8 MHz. 
Also, a phase difference of 1 6° exists between the fifth-order IM3 and IM5 
contributors at high modulation frequencies. It is important to note that 
these values are closely dependent on the nonlinearity coefficients and on 
the filtering and that the result presented here cannot be generalized. It is 
evident, however, that deviations from the memoryless approximation are 
unavoidable, if the maximum modulation frequency is in the MHz range. 

One important feature of amplitude-dependent memory effects arises 
from the third harmonic band. Let us assume that all frequency bands up to 
the second harmonic are flat over the whole range of distortion bands, so 
that no memory effects exist at moderate amplitude levels, because the 
higher than third-degree nonlinearities are practically zero. If the signal 
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amplitude is now increased, fifth-order mechanisms will start to have an 
effect, and signal components from the third harmonic band will be 
converted down to IM3 . This means that if only the third harmonic band is 
tilted, memory effects will occur exclusively at high signal levels, where 
fifth-order distortion begins to play a significant role. 



3.6 Summary 

Once the actual transfer function of the amplifier is affected by the 
bandwidth of the applied signal, it will exhibit memory effects. There are 
some ways of measuring these effects, and in this book intermodulation 
signals as a function of tone difference of a two-tone signal are measured. 
If the amplitude and/or phase of the IM signals is affected by the tone 
difference, the amplifier exhibits memory effects. 

Smooth memory effects are not usually harmful to the linearity of the 
PA itself. A phase rotation of 10° to 20° or an amplitude change of less than 
0.5 dB as a function of modulation frequency has no dramatic effect on the 
linearity of the device, but the situation is completely different when 
predistortion type linearization is employed to cancel out the IM sidebands. 
This undertaking requires an extremely accurate amplitude and phase 
matching between the distortion components and the cancelling signals. If 
the IM3 components rotate as a function of modulation frequency, for 
example, but the canceling signals do not, the cancellation performance of 
the linearization method may be inadequate for wideband signals. Different 
linearization techniques have different sensitivities to memory effects, and 
feedback or feedforward type linearizers, for example, are not very 
sensitive to amplifier behavior. However, RF and digital baseband 
linearizers have potential for both high efficiency and wide bandwidth - 
provided that the deterioration in the cancellation performance due to 
memory effects can be cured. 

This chapter introduced two memory effects: electrical and thermal. 
Electrical memory effects are produced by nonconstant node impedances 
within frequency bands. Most of these effects are generated by a frequency- 
dependent envelope impedance, caused usually by the bias impedance. 
Thermal memory effects, for their part, are generated by the junction 
temperature, which is modulated by the applied signal. Since the chip 
temperature varies at the envelope frequency and some of the electrical 
parameters of the transistor are affected by the temperature, gain 
modulation and IM3 components are generated. These memory effects 
occur, because the temperature rise caused by dissipated power is highly 
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dependent on the modulation frequency, which also determines the 
behavior of the IM3 components. 

If the signal amplitude is reasonably large, fifth- and higher order 
distortion mechanisms will have an impact on IM3 performance. Fifth- 
order distortion will affect both the IM3 and IM5 sidebands, and in the 
polynomial case, the fifth-order IM3 and IM5 contributions will have 5 : 1 
amplitude ratio and a phase difference of zero. Amplitude -dependent 
memory effects that change this situation, however, arise from two 
phenomena: first, the frequency bands become wider, because fifth-order 
signal components are wider than third-order components alone and, 
second, also the third harmonic band converts down to IM at high signal 
levels. These effects are called amplitude-dependent memory effects, 
because the amount of fifth-degree nonlinearity in the IM3 components is 
dependent on the amplitude, and the total extent of amplitude-dependent 
memory effects is the sum of the memory effects converted from all 
frequency bands. 

Memory effects can be visualized by a vector presentation of IM3. 
Presenting the contributors to IM3, such as the envelope, fundamental, and 
second harmonic components, as individual vectors instead of the total 
result, makes it possible for us to identify the causes of the memory effects. 
This book will make extensive use of this presentation format. 



3.7 Key Points to Remember 

1. Efficiency and linearity are opposite design goals in traditional RF 
power amplifier design. 

2. The idea of linearization is that the power amplifier itself is designed 
to be not linear enough to achieve good efficiency, after which the 
linearity requirements are fulfilled by external linearization. 

3. The implementation and complexity of the linearizer affect the 
overall efficiency of the transmitter. 

4. The cancellation performance of especially the simple predistortion 
type linearization techniques is sensitive to memory effects generated 
in the RF power amplifier. 

5. There are two types of memory effects: electrical and electrothermal. 

6. Distortion appears as currents, and the terminal voltages caused by 
distortion can be shaped by the terminal impedances. 
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7. Lower order distortion may mix up to higher order distortion (i.e., 
A/th order distortion can be generated by lower degree nonlinearities). 

8. Electrical memory effects are caused by nonconstant terminal 
impedances at dc, fundamental, and harmonic bands, of which those 
within the dc band are the most harmful, because bias impedances are 
strongly frequency dependent. 

9. The thermal impedance describes the temperature increase caused by 
dissipated power, and it has a wide range of time -constants and a 
significant magnitude up to the megahertz frequency range in RF 
power transistors. 

10. Electrothermal memory effects are caused by dynamic temperature 
variations at the top of the chip that modifies the electrical properties 
of the transistor at the envelope frequency (tone difference 
frequency). As a result, IM3 signals that depend on thermal 
impedance are generated. 

1 1. Fifth-order distortion affects IM3 components, and the memory 
effects of fifth-order distortion can also be seen in IM3 components. 

12. Since fifth-order effects on IM3 are an amplitude-dependent 
phenomenon, memory effects also become amplitude dependent. 

1 3 . With higher-order distortion, the frequency bands get wider and new 
harmonic bands mix down to fundamental. This causes the memory 
effects to vary with signal level. 
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Chapter 4 

The Volterra Model 



In this chapter, we will first review the most common nonlinear modeling 
techniques. Section 4.1.2 discusses the properties of the polynomial 
Volterra modeling technique, and Section 4.2 explains in detail how the 
purely electrical and electrothermal nonlinearities are written and what 
terms are needed in the series expansion that models the nonlinearity. 
Section 4.3 illustrates how distortion of a common emitter amplifier is 
calculated using the Volterra analysis technique, and as a first case study, 
Section 4.4 presents the tear-down analysis of all the terms affecting the 
IM3 distortion in a BJT CE amplifier. Similar analysis is performed to a 
MESFET amplifier in Section 4.5. 



4.1 Nonlinear Modeling 

To be able to analyze the nonlinear behavior of an amplifier, we need 
accurate models for both the active, nonlinear circuit and the passive 
matching and biasing components. Both of these may be problematic to 
obtain. It is a well-known fact that the amount of Mh-order distortion is 
proportional up to Mh-order derivatives of the I-V and Q-V curves (see [1- 
3]). Hence, for accurate distortion simulation the I-V and Q-V curves of the 
active components must be modeled so that not only the dc value, but also 
the higher order derivatives are correct and continuous (for reference, in 
early simulation models, already the first derivatives could be 
discontinuous [4]). Moreover, capacitances are easy to model so that charge 
is not conserved, which may result in nonphysical rectification and self- 
biasing in purely capacitive nodes. Hence, especially if the capacitance 
values depend simultaneously on two terminal voltages, it is important to 
model the capacitances as charge equations instead of capacitances [5], 
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The passive components are also tricky to model at RF frequencies, due 
to their distributed nature at high frequencies. Lossy transmission lines are 
difficult to model for time-domain simulations, and, in general, the 
modeling of the passive components tends to be more accurate in the 
frequency domain. Still, the frequency-domain simulation models of some 
passive components like step changes in transmission line width may also 
be inaccurate at higher harmonics. As one extreme of modeling, it is 
possible to use the measured terminal impedance values of a test circuit in 
the circuit analysis. 

In short, for accurate distortion simulation we require from the 
simulation models that: 

1 . The Ath-order derivatives of the I-V and Q-V curves must be accurate 
enough for Ath-order distortion simulations. 

2. The frequency responses of the node impedances must be correct up to 
the highest relevant harmonic. At baseband frequencies, correct 
modeling of biasing impedances and the thermal impedance are also 
needed. 

3 . It would also be very helpful if componentwise information about the 
dominant distortion sources were available. 

4.1.1 Nonlinear Simulation Models 

Broadly speaking, there are two types of nonlinear models that PA and 
transmitter designers use: behavioral black-box models for system 
simulations, and device models for circuit simulations. These can be 
further divided based on the modeling technique: models may be analytical, 
based on some predefined and physics-based, parametrized model 
functions, or empirical, where the measured data is tabulated and 
interpolated or modeled with simple splines or polynomials with no clear 
physical meaning. These main groups are illustrated in Table 4.1. 

Behavioral baseband models are widely used for simulating and 
optimizing entire transmitters and transceivers, and new flavors have been 
added, for example to model memory effects, as in [6], However, 
behavioral models describe either an existing amplifier, or they are used to 
derive specifications for a so-far-nonexistent amplifier, but they are of 
limited use in the design of a new power amplifier. 

Just for reference, the properties of the most common behavioral 
models are described here briefly. Plain static AM-AM and AM-PM curves 
are not capable of modeling memory effects, but models where AM-AM 
and AM-PM curves depend on the modulation frequency have been 
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developed. In the Saleh model, linear filters are added both to the input and 
output side of the AM-AM and AM-PM nonlinear block, as illustrated in 
Figure 4.1. In the Blum & Jeruchim model, FFT and adequate over- 
sampling is used to find the instantaneous modulation frequency that is 
used to modify the AM-AM table (both models are described in [7]). 
Furthermore, Cadence has implemented its own behavioral K-model in its 
SpectreRF simulator [8] . An example of Volterra type behavioral modeling 
is called Volterra input output map (VIOMAP). It is conceptually a 
nonlinear extension of normal ^-parameters, including harmonic responses, 
and it has been successfully used in single-tone load-pull simulations [9, 
10 ]. 



freq. dependent nonlinear amp. 




freq. indep. 

nonlinear 

amplifier 




Figure 4.1 (a) AM-AM and AM-PM curves for a power amplifier, and (b) a 

frequency-dependent nonlinear model based on filters and memoryless 
nonlinearity. From [11], 
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Device models describe the operation of the semiconductor device, and 
together with proper models for passive and distributed components, 
models for amplifiers can be built and optimized. 

Early semiconductor models were analytical, using equations that were 
first derived in detail from semiconductor physics and then often simplified 
to reduce the simulation time. The basic problem with these equation based 
models is that the chosen functions and control parameters fix the possible 
shape of I-V and Q-V characteristics, and there may simply not be enough 
degrees of freedom to model things like / c - V CE curvature. As an example, 
in a slightly simplified form the collector current in the basic Gummel- 
Poon (GP) BJT SPICE model is given by 



In = 



(— )■ 


fl + Vc *- V ^ 


v v,) 


l VAF ) 



4 • /c 



( v BE 

’ ' eX P(^— 

IKF 



(4.1) 



where the basic exponential At^be dependency can be modified by only 
three control parameters: IS scales the current, VAF (so-called Early 
voltage [12], illustrated in Figure 4.2) makes an extremely simplified model 
of the output conductance, and IKF (so-called knee current) reduces the 
gain at high currents [13, 14], This simple equation covers the entire I-V 
plane and also fixes the derivatives d^/^/dK 11 , thus fixing the nonlinear 
behavior. The SPICE GP model can still be used reasonably well for 
simulating the fundamental signal [15], but especially due to the over- 
simplified and inherently linear output impedance model, it cannot be used 
for accurate distortion simulations, as will be seen later in this book. 

Better physical models have been developed, and for example Mextram 
and VBIC for BJTs [16-20] and Motorola MET model [21] for laterally 
diffused MOS (LDMOS) devices are quite sophisticated. Compared to 
early SPICE models, they are greatly enhanced. The I-V curve has a more 
realistic shape and continuous higher order derivatives, charge conserving 
and continuous capacitance models have been included, and the effects of 
self-heating have been added into the model. The latter is important for 
finding the correct dc bias, as the self-heating makes a big difference in the 
I-V curve. It can also be used to model the thermal memory effects, 
provided that the thermal model has enough time constants to model both 
the slowly warming package that affects mainly the dc bias and the 
microsecond range thermal memory of the chip’s surface. 
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Figure 4.2 Output impedance modeling in a BJT using the Early voltage VAF. 



Added degrees of freedom in the model increase the model complexity 
and the number of control parameters. In one extreme, MOS BSIM models 
have tens of parameters controlling the scaling properties alone. Hence, the 
complexity of the models tends to grow out of hand and their fitting gets 
increasingly complicated and sensitive to errors. 

Another approach in device modeling is to abandon the equations and 
use tabulated measured data or fully empirical fitting functions instead. 
Now any form of I-V and Q-V characteristics can be modeled, and this is 
the idea behind the Root models, sometimes called “the device knows best” 
models [5], There are some technical problems in interpolating tabulated 
data, as interpolating polynomials easily create oscillations between the 
data points and hence nonphysical fluctuation in the higher order 
derivatives. However, tabulated models are very flexible and quite easy to 
fit as no forcing to the predefined functions is needed. 

Volterra models are empirical models as they do not rely on 
semiconductor physics. The nonlinearities are described as polynomials, 
the coefficients of which may be obtained either by differentiating physics 
based I-V and Q-V functions, or by fitting polynomials directly to 
measured data. The latter approach is used here, and the properties of 
Volterra models will be studied more deeply in the following section. 

4.1.2 The Properties of the Volterra Models 

Polynomial models are not automatically quick to simulate; on the 
contrary, they may converge badly at signal levels higher than the original 
fitting range. However, the use of polynomial modeling allows the use of a 
noniterating and efficient Volterra analysis procedure. 
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However, here the main motivation for using the Volterra simulation 
technique is not the speed advantage but the fact that it provides an 
excellent tool for analysis. Dominant distortion mechanisms can be 
recognized in the same way as in normal ac noise analysis, and owing to 
the nonlinear analysis, multiple mixing mechanisms can also be 
recognized, which aids the design of harmonic terminal impedances, for 
example. Thus, Volterra analysis is one of the few ways of obtaining an 
understanding of memory effects and aiding design optimization. 

Still, polynomial modeling has some shortcomings that need to be 
recognized. First, polynomial models suffer from the fact that outside the 
fitting range their response explodes towards infinity. This is the opposite 
of typical nonlinear modeling functions, where smooth and limited 
behavior over the entire bias range is a desired property, as it aids 
convergence and the signal swings are not necessarily a priori known. 
Hence, the Volterra analysis is not a very general tool. Due to its speed it is 
used for quick distortion analysis and optimization for low-noise amplifier 
(LNA) type small-signal circuits in simulators like Voltaire XL [22] and 
early versions of SPICE, or even as standalone simulators [23], However, 
for successful power amplifier analysis, certain preliminary information is 
needed. 

Second, the actual large-signal dc bias voltages are needed in advance. 
The large signal operation often causes a shift in the dc operating point that 
affects both the gain and the amount of nonlinearity. This signal-induced dc 
shift slows the convergence in a harmonic balance simulation, and in the 
noniterating Volterra calculation procedure it can only be estimated but not 
completely predicted. To overcome this, we need either to check if the dc 
shift is significant or to fit the polynomial model at the actual large signal 
operating point. 

Third, during the fitting of the polynomial functions the extent of input 
and output voltage swings is needed. The real power of the polynomial 
modeling is that - besides the separation of distortion components - a local 
fit over just the required voltage range can provide more accurate high- 
order derivatives than the use of analytical models that have to cover a very 
broad range of operating regimes with a limited number of control 
parameters. The larger the fitting range, the less accurate a low-degree 
polynomial fit may be. Hence, it is desirable to fit along the desired 
maximum signal excursions, and not much beyond as it compromises the 
accuracy of polynomial modeling, and also not over a smaller range as the 
response of the polynomials may then be completely nonphysical outside 
the fitting range. In that sense, a good guess of the input and output 
trajectories is needed. Altogether, Volterra analysis is not necessarily an 
easy-to-use standalone simulation approach, but it provides much 
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debugging power when used in parallel with other simulation methods such 
as harmonic balance, for example. 

In this book, the examples studied are limited to single-transistor 
amplifier stages, and the Volterra analysis has been calculated semi- 
analytically, by deriving symbolically the transfer functions from each 
distortion source to all node voltages. Symbolic analysis is by no means 
necessary and limits the analysis to fixed CE or CS amplifier architectures 
and two-tone test signal. Instead, the Volterra analysis can - with almost 
the same resolution of distortion components achieved here - be calculated 
fully numerically by recursively using ordinary linear ac analysis on any 
circuit described by a standard modified nodal analysis (MNA) matrix 
representation and nonlinear current sources. To aid the study of more 
complicated amplifier topologies, the effects of multidimensionally 
controlled charges, and other complicated features, a fifth-order, multitone 
numerical simulator with a simple netlist interface has been designed [24], 



4.2 Nonlinear I-V and Q-V Characteristics 

Most transistor models are based on either a n or T model. Here, a n 
model is used, and this section presents the typical conductive (I-V) and 
capacitive (Q-V type) nonlinearities appearing in the n model of BJT, 
heterojunction BJT (HBT) and field-effect transistor (FET) devices. The 
BJT is used as an example, but the same models with different sets of 
polynomial coefficients can be used for FET transistors as well. 

As explained before, the Volterra model is based on polynomial 
modeling of I-V and Q-V curves. Measuring these may be somewhat 
complicated, as seen in more detail in Chapter 5. Charge as such cannot be 
measured directly, and we must rely on extracting the capacitances from ac 
measurements and integrating the charge equations from the obtained 
capacitance values. In a similar manner, the I-V curve can be mostly 
reconstructed from g m and g 0 values obtained from 5-parameter 
measurements, but the actual I-V curve is a safer starting point. 

The models presented are electrothermal models, which means that the 
junction temperature also appears as a free variable. However, the dc 
temperature rise is included in the bias point, and only the temperature 
variations caused by dynamic self-heating are considered. As the dissipated 
power is a product of voltage and current, the fundamental assumption in 
the following is that the ac component in the junction temperature is 
already a second-order phenomenon. Hence, a third-order model will 
contain only the first power of the temperature, which also means that the 
temperature dependencies of circuit elements are considered linear. 
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4.2.1 Ic- V BYT V CE Characteristic 

In most reported BJT/HBT Volterra series analyses the collector current is 
considered a function of base voltage only [25-27], which captures the 
dominant exponential input-output nonlinearity but assumes that the output 
conductance is constant. In MESFET Volterra series analysis the effects of 
drain voltage are usually implemented by a polynomial for g 0 (see [28]), 
but even that does not capture all nonlinearities. 

Already the simple expression in (4.1) is a three-dimensional function 
of V BE , V CE , and junction temperature T, as V t = kT/q. A polynomial model 
is derived simply by expanding the large-signal I-V function 

i c = nv BE>r C E>r), (4.2) 

to a three-input Taylor series around the dc operating point E BE q, F CE q, 
7q. Hence, an electrothermal third-degree series expansion of the ac 
current can be written as 

h “ &m v be + K 2GM v be + ^3GM v be 

+ &o v ce + K 2GO v L + K 3GO v ce ^ 

+ ^2GMGO v be v ce + ^3GM2GO v be v ce + ^3GMG02 v be v ce 
+ ^2GMT^J + ^3GMT^J v be + ^3GOT ? J v ce 



where v be - v BE - V BE q, v C 6 -v ce -F CE q and tj - T-Tq, and A jxxx is the zth 
degree nonlinearity coefficient of element xxx (alternatively, V MNP could 
be used to mark a term v be M v ce N Zj P ). The first row models the effect of v be 
alone, and the second row the effects of v ce alone (i.e., a nonlinear output 
conductance). These are not sufficient, though, as cross-products of v bc and 
v ce may also appear, and these are listed in row 3. Finally, the temperature 
change t } contributes to the current, and it, too, may mix with both terminal 
voltages, causing the last three terms on row 4. 

The effects of the electrical nonlinearities are demonstrated in Figure 
4.3, where the collector current at three base voltages is plotted as a 
function of collector voltage at three different base voltages. If all the 
coefficients except the g m are zero, we obtain the three equally spaced 
horizontal lines shown in Figure 4.3(a). Since the lines are exactly 
horizontal, the output conductance is zero and the collector voltage does 
not affect the amount of current. Furthermore, since the lines are equally 
spaced, the transconductance is linear. However, if K 2 qm or AT 3GM deviates 
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Figure 4.3 Demonstration of collector current nonlinearities. Vertical axis is the 
collector current and horizontal axis Vq E voltage, (a) Linear response, 
(b) non-zero g D , (c) non-zero V 3G0 , (d) non-zero ^GMGO ( e ) non-zero 
^3GM2GO ancl (f) non-zero V 3 G mg02- 



from zero, the lines in a I-V plane become unequally spaced, indicating a 
nonlinear transconductance. 

The effects of g 0 are demonstrated in Figure 4.3(b), in which just the 
linear terms g m and g 0 exist. Compared to Figure 4.3(a) the lines now have 
a nonzero slope that is proportional to g 0 and independent of vj, e . 

While Figure 4.3(b) is still fully linear, the nonlinearity of the output 
conductance is demonstrated in Figure 4.3(c), where the slope of the 
current varies with Vq%. In this case only K^qq has a nonzero value in 
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Figure 4.3(c), but both K 2 qq and K 2GO can be used to model curvature 
effects of the output conductance such as saturation and breakdown. 

Figure 4.3(d-f) illustrates the effects of the cross-terms ^2GMGO 
A 3 gm 2GO ar| d ^3GMG02> respectively, that model the interaction between 
base and collector nonlinearities. To aid comparison, the thin lines in 
Figure 4.3(d-f) are copied from (c) where only g m , g 0 and K^ G0 have 
nonzero values. ^2GMGO (corresponding to the v be v ce term) is nonzero in 
Figure 4.3(d), and as a result of it the slope of the lines varies not only as a 
function of collector voltage as in Figure 4.3(c), but also with the base 
voltage. This is essentially needed to model the Early effect illustrated in 
Figure 4.2. Similar reasoning applies to ^ 3 gmG 02 ar| d ^3 gm 2GO> that shape 
the output conductance as functions of v be v ce and v be 2 , respectively, as 
shown in Figure 4.3(e, f). 

The I-V curves of different modeling approaches are further compared 
in Figure 4.4. If the collector current is modeled as a one-dimensional 
function of base voltage and linear g 0 , just straight line I-V curves are 
generated, as seen from the thin dashed line in Figure 4.4. The I-V curves 
simulated using the SPICE Gummel-Poon model are also straight (and 
thin) lines, but their slope and hence the output conductance vary with the 
collector current, as suggested in Figure 4.2. In reality, however, BJT I-V 
curves are far from straight lines under large or semi-large signal 
conditions, due to quasi-saturation and breakdown effects. The curvature 
can be modeled by using one-dimensional polynomials for both v be and v ce , 
as illustrated by the thick solid line in Figure 4.4. However, elementary 
phenomena such as the Early effect cannot be modeled without the cross- 
terms that cause the Iq-Vq E slope to depend on the value of V BE . This is 
illustrated by the thick dashed line, corresponding to a full series expansion 
(4.3). The shapes of the saturation and breakdown also depend on the base 
voltage, and this makes the use of cross-terms mandatory to avoid large 
errors at the comers of the I-V plane, as seen in Figure 4.4. 

Finally, the electrothermal effects of the collector current are 
discussed, modeled by the last three terms in (4.3). A second-degree term 
V 2 gmt models a temperature-dependent shift in the current, as illustrated 
in Figure 4.5(a, b). It is worth noting that V 2GMT cannot be derived from 
the small-signal parameters g m and g 0 ; instead, actual current 
measurements are needed. V 3GOT is a third-degree term that includes a 
combined effect of temperature and collector voltage as indicated in Figure 
4.5(c). Essentially, it models the temperature dependency of the output 
conductance. Similarly, .fiT 3GMT models the combined effects of temperature 
and base voltage, visualized in Figure 4.5(d). Since the slope of that curve 
describes the transconductance, V 3GMT can be considered a change in the 
transconductance as a function of temperature. 
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4.2.2 g pi and r bb 



The I D -V DS characteristic is the only important conductive nonlinearity in 
FET-type transistors. In BJTs, two other conductive nonlinearities exist: 
The I B -V BE nonlinearity caused by exponential g pj and the nonlinear r bb . 
The effect of g pi conductance is usually more important, and it is also 
easier to model. In theory the I B -V BB equation should be roughly (4.1) 
divided by the current gain (3, but some simplifications can be made. Since 
the base current does not depend strongly on the collector voltage, we can 
use a two-dimensional model of v be and junction temperature t } dependence 
only: 



z 2> Spi ' v be + ^"2GPI ' v be + ^3GPI ' v be (4-4) 

+ ^2GPIT ' h + ^3GPIT ' h ' v be 

Here, the coefficients have similar meanings as before. The linear term is 
modeled by g pp and K 2 q?\ ar| d ^3GP\ model its exponential curvature. 
Again, A^gpit models the I B shift caused by self-heating, and ^qpit can 
be seen as a temperature dependence of the linear g pi term. 

The intrinsic base resistance r bb is a bit trickier to model. It is a series 
resistance between the intrinsic and extrinsic base points, but its value 
depends on the current crowding in the base area and also on the value of 
the intrinsic v be . Thus, it must be modeled as a three-dimensional 
conductance, being controlled by the voltage across the resistor (v bb =v bext - 
v bblt ), the intrinsic base voltage v bevnt , and the junction temperature t } . All 
the v be j nt * terms (b= 1,2,...) are zero, but the current crowding effect is 
modeled with the cross-terms between v bb and v bejnt , as shown in (4.5). 
However, r bb is usually small and it has been modeled simply as a linear 
conductance in the following examples. 

z rbb = &bb ' v bb + ^2 11 ' v bb ' v beint + - ( 4 - 5 ) 

4.2.3 Capacitance Models 

As explained before, the capacitances are modeled as polynomial charges 
that are then differentiated with respect to time to get the displacement 
current. The charge may be - and often is - controlled by more than one 
port voltage, in which case a multidimensional polynomial similar to (4.3) 
must be used. The charge may also be modeled as a transcapacitance, and 
in this case the charge does not appear between the controlling nodes but in 
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some other node. In the following examples, only one controlling voltage is 
assumed, and as an example, (4.6) represents the base-to-emitter charge as 
functions of base-to-emitter voltage and temperature. 

Qbe = ^pi ' v be + ^2CPI ' v be + ^3CPI ' v be (4-6) 

+ ^2CPIT ‘ h + ^3CPIT ' h ' v be 

From this equation, the corresponding measurable capacitance C ^ and the 
nonlinear current source z'nlcpi are obtained simply by differentiating the 
charge equation (4.6) with respect to v^ e and time, respectively. 

c pi ( < v be ) = C pi + 2 ^2CPI ' v be + 3 ^3CPI ' v be + ^3CPIT ' tj- ( 4 - 7 ) 

Z NLCPI = J 03 ' ( C pi ' v be + ^2CPI ' v be + ^3CPI ' v be ( 4 - 8 ) 

+ ^CPIT ' tj + ^3CPIT ' tj ' v be) 

In (4.8), CO is simply the frequency of the distortion tone; hence, 
capacitances do not cause dc distortion currents but distort most heavily at 
the harmonic frequencies. Equation (4.7) shows that the temperature- 
dependent charge term .K^CPIT cannot be derived from capacitance 
measurements; still, a time-varying junction temperature may cause a 
current proportional to it. Otherwise, the first term C pi in (4.6) describes the 
small-signal capacitance, and AT 2 cpi an d ^3 cpi define its electrical 
nonlinearity, ^ 3C p IT describes the charge being a function of both the 
controlling voltage and junction temperature and since C = dg/dv, the 
effects of ^ 3 cpit can be seen as the temperature dependency of the 
capacitance value. 

As seen from (4.6), a linear C-V trend /f 2 CPi causes quadratic charge 
nonlinearity. Similarly, a capacitance proportional to v 2 (A^cpi) causes 
cubic nonlinearity. Different types of capacitances have different 
characteristics, as illustrated in Figure 4.6(a). The base-emitter capacitance 
Cpi is exponential [see (4.27)] and therefore highly nonlinear. Reverse- 
biased p-n or Schottky junctions seen in BJTs and in FETs are just weakly 
nonlinear, and they can be linearized further by increasing the reverse bias. 
MOSFET-type transistors have peculiar gate capacitances and, for 
example, the C GS makes a clear dip around the threshold voltage. If the 
MOSFET is operated close to cutoff, this dip can cause a large amount of 
second-degree nonlinearity. 
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(a) V BE ! V GS / V CB IV] (b) V BE I V GS / V CB t V l 

Figure 4.6 Normalized (a) capacitances and (b) charges of BJT C BE , MOSFET C GS , 
and BJT C BC . 



4.3 Model of a Common-Emitter BJT/HBT Amplifier 

The direct method is now used to calculate the IM3 components in a 
common-emitter BJT/HBT amplifier, using the procedure outlined in 
Section 2.5.2. The analysis starts by building a model for the circuit, from 
which the fundamental amplitudes are found by a linear ac analysis. Then, 
the second-order currents and voltages and third-order currents and 
voltages are calculated using the procedure shown in Section 4.3.2. 

4.3.1 Linear Analysis 

The model of a common-emitter BJT amplifier, shown in Figure 4.7, 
includes the input impedance Z IN (lumping both the matching network and 
bias circuitry), base-emitter conductance g pi and capacitance C pj , feedback 
capacitance C BC , output capacitance C CE and output conductance g 0 , 
transconductance g m , load impedance Z L , and emitter impedance Z E . The 
input and load impedances include not only the impedances of the 
matching networks, but also the impedances of the bias networks and 
package parasitics, and Z IN further consists of the output impedance of the 
preceding stage and the intrinsic base resistance r bb , as shown in Figure 
4.7. 

To reduce the amount of equations, the input voltage source can be 
replaced by its Norton equivalent source 
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Z IN = 7 In( 5 )‘ V IN’ ( 4 - 9 ) 

and using the following shorthand notations 

= ^pi + 5 ^pi 

7 ce( s ) = §o + sC ce- ( 4 - 10 ) 

^BC^ 5 ) = 5( ^BC 

The matrix equation can be set up by changing all the impedances in 
Figure 4.7 to admittances and applying Kirchhoff’s current law at nodes 1 
to 3 . This results in 



(4.11) 



z 'lN 




^IN + ^BE + ^BC ^BE ^BC 




V B 


0 


.■§* 


SnT^E ^m + ^ E + ^E +7 CE ^CE 




V E 


0 




Sm-^BC -g m -^ C E ^L + ^CE +7 BC 




V C 
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and the responses to i IN can be found using Cramer’s rule. Thus, the 
voltages at the base, emitter, and collector are given by 

( Y BE Y CE + Sm Y \ J + V L + Y CE Y E + ^CE^V| i IN (j) 

v„(s) = • -P-— (4.12) 

l *t L Y E +Y BC (Y BE *T QE +Y E + gm ) J 

_ ^be^ce + ^ l + ^be^ l + ^bc(^be + ^ce + ^) „ 

V * (,) = d^T) Ws)(4-13) 



and 



v c (s) = 



^BE^CE Sm Y E + Y BC^ Y BE " 



■ y e + 8J 



det(s) 



-h M (4.14) 



where the determinant of the admittance matrix is written as 



det(s) = Y be Y ce (Y l + T e + Y m ) + Y m Y L (4.15) 

• (^BE + ^CE + Sm + Y e ) + We^IN + 

+ ^BC • [ V IN + Vl N + VlN + W E 

+ V L + ^CE Y L + Y E Y L + Y CE Y E + g m Y L 
+ &m 7 IN + ^m^E] 



From these, the base-emitter and collector-emitter voltages are simply 
v BE (s) = v B (s)-v E (s) and v CE (s) = v c (s)-v E (s), respectively. Finally, as both 
v BE and v^ E are frequently used in calculating the distortion generated by 
the g m element, for example, it is handy to derive their ratio 



TF(s) 



v CE ( ,y ) 

v Be(‘ s ') 



(4.16) 



The purpose of the linear analysis is to obtain the fundamental voltage 
amplitudes across all nonlinear components so that we can proceed in 
calculating the nonlinear currents generated in these components. Before 
doing that, a few observations about the signal swings are worth making. 

The exponential response of a BJT is extremely nonlinear and it does 
not tolerate higher than 10 to 30 mV signal amplitude in the BE junction 
without distorting excessively. That does not sound like very much for a 
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power amplifier, but a couple of things happen to help the situation. First, 
the device is not completely exponential, but when driven to high injection, 
a BJT linearizes considerably as modeled by the parameter IKF in (4.1). 
Second, the amplifier has several feedback mechanisms that reduce the 
signal level in the BE junction. The series emitter impedance causes a 
linearizing series feedback, and C BC causes a shunt feedback. The effect of 
C BC is very important, as the strong capacitive feedback considerably 
lowers the impedance at the base, and so reduces the BE voltage swing and 
the amount of generated distortion for a given driving power. 

4.3.2 Nonlinear Analysis 

In this section the nonlinear model of a CE BJT amplifier is presented and 
the equations for IM3 distortion are derived. The circuit has three two-input 
and one three-input I-V and Q-V nonlinearities, modeled by 27 first-, 
second-, and third-degree coefficients, of which 18 are purely electrical 
(Cpi* K 2CPb ^3CPI> Q)C> K 2CBC’ K 3CBC’ Spv K 2GPb K 3GPb K 2GM’ ^3GM> 
^2GO ^3 GO’ ^2GMGO ^3GM2GO’ an( l ^3GMG02) an d n i ne are related to 
dynamic temperature variations ( AT 2 C l J r I- A 3 (' 1 * 11- ^2 ( I H i' ^ 3 t i h i ■ -^2GPIT, 
^2GPIT’ ^2 GMT’ ^23 GMT’ an( i ^3 GOt)- 

Eventually, this analysis will present the IM3 tone as a vector sum of 
seven terms caused by cubic electrical nonlinearities, 21 terms caused by 
cascaded quadratic nonlinearities where the rectified envelope information 
is upconverted to IM3, 21 similar terms where the second harmonics are 
mixed down to IM3, and finally, five cubic and 24 cascaded second-degree 
electrothermal terms. This resolution may seem excessive, but it was 
chosen to illustrate the real multitude of different mechanisms that generate 
distortion; it also clearly illustrates that very much information is lost if 
only the effects of the cubic terms are analyzed. With larger circuits and 
higher order analysis it will be necessary to compress the data somehow, 
but the principle is still the same: we want to see how much of the total IM3 
is caused by mixing distortion from the dc or harmonic bands and to be 
able to minimize the total distortion (or its memory effects), we want to see 
what nonlinearities and impedances actually are causing the distortion 
voltages at these harmonic bands. 

Due to the multitude of terms, not all of them will be discussed 
separately. The full analysis is shown in Appendix C, and next, the 
calculation procedure is illustrated with some examples. 
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4.3. 2.1 Second-Order Distortion Currents 

The circuit for solving the second-order responses is shown in Figure 4.8, 
where the linear input voltage shown in Figure 4.7 is now short-circuited 
and the second-order distortion current sources %L2XX are added in 
parallel to all nonlinear circuit elements. Currents ending with T are 
electrothermal currents that are discussed in Section 4. 3.2. 5. As before, 
Z IN> Z E> and Z L lump the package parasitics, and biasing and matching 
impedances. 

For calculating the self-heating effects, the instantaneous power 
dissipation is calculated as P^ ss = v CE z' c , and the thermal impedance Z TH 
shown in Section 3.4 is used to calculate the instantaneous temperature 
fluctuation 7 , |(co 2 co | J at frequency co 2 — <» i - It is possible to use different 
temperatures for different circuit elements, but as they all are physically 
located close to the base area, a common temperature is used here. For 
large devices, however, it may be advantageous to split the device into 
smaller parallel devices that may see different temperature variations. 

The Tct^be'^ce nonlinearity is modeled as a three-dimensional 
function of the v BE , v CE , and temperature, including g m and g 0 nonlinearity 




Figure 4.8 Representation of a circuit containing current sources for second-order 
responses. 
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and all up to third-degree cross-terms. The g pi and C pi are nonlinear 
functions of the base-emitter voltage and temperature, and the weakly 
nonlinear C B c is controlled by the collector-base voltage and temperature. 

We start the analysis by calculating the second-order distortion currents 
; N l 2 XX- As an example, the second-order envelope current at co 2 -Cfl] caused 
by the g pj is, using Table 2.5, 



Z NL2GPl(®2 W l) ^2 GPI E ( ® 2 )K be (co,). (4.17) 

As another example, the second-order envelope current caused by the I c - 
^BE-^CE nonlinearity (g m ) is 



z NL2Gm( C 0 2 “ M l) = ^GM^BE^^^BE^l) > (4-18) 

+ ^2GO^Ce( (0 2)^Ce( G) i) 

+ 0.5 • ^2GMGC>( K Be(® 2) F CE( W l) + ^Ce(®2) 



which combines the effects of the second-degree input nonlinearity 
^2GM* V BE 2 ’ output nonlinearity 7f 2 GO* v CE 2 and the input-output cross- 
term ^2GMGO* v be* v ce> ail seen > n the I-V model (4.3). 

As noted from Table 2.5, the selection of the frequencies of the phasors 
and the values of possible constant terms depend on the tone frequency: a 
product F BE (col)*F BE (col) results in a tone at 2 col, for example. Above, the 
tones of the phasors are chosen so that they always result in distortion at 
the envelope frequency co 2 — cOj. The phasors F BE and V CE for the 
fundamental tones (Oj and (0 2 are calculated using (4.12)-(4.14). 



4.3. 2.2 Transimpedance Transfer Functions and Second-Order Voltages 

Next, we need to convert the distortion currents to distortion voltages in 
various nodes. Here, a semisymbolic analysis has been chosen, so that the 
transfer functions TF XY z f rom a nonlinear current source between nodes X 
and Y to a node voltage Z have been derived by hand. The general idea is 
that within each nonlinear element, the amplitudes of the distortion 
currents are calculated using the lower order voltage phasors, and using the 
transfer functions, the generated distortion currents are then converted to 
distortion voltages in the chosen nodes. The transfer functions can be 
derived from (4.11) by replacing z IN with a test current source between 
nodes X and Y. As an example, the transfer function 7Y BEB is of the form 
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V T3 

TF beb (s) = j* 

1 BE 

_ -IY e -Y l+ Y ce • Y e + 7 ce • 7 l + g m . r L + 7 bc • 7 e ] ( 4 - 19 ) 
det(s) 

where det(s) is given by (4.15). A complete set of transfer functions is 
shown in Appendix C. Using these notations, the complete second-order 
envelope voltage at the base node is now given by 

U B2 (C0 2 -0)i) (4.20) 

F F bbb ((£>2 ~ ®l) ' ( Z NL2GPI + Z NL2CPI + Z NL2GPIT + z NL2CPIt) 

+ TF CEb(®2 ~ “l) ' ( Z NL2GM + z NL2GMt) 

+ TF CBB (® 2 “ W l) ' ( Z NL2CBC + Z NL2 CBCt) 

where each distortion current source is multiplied by the proper 
transimpedance that converts the current to a voltage at the base. Note also 
that there are eight different current sources (four of them being 
electrothermal and marked with 7) and each one of them will have an 
equation resembling (4.17) or (4.18). 

To solve the second harmonic voltages, we need to rewrite the current 
equations of the ; NLX x sources for the desired harmonic, and to recalculate 
(4.20) at that frequency. After this procedure, we have the base envelope 
and second harmonic tones given as sums of eight different contributions. 
By storing them separately, we can see which one is dominant, and which 
ones may cancel each other. 

4.3. 2.3 Solving the Third-Order Voltages 

Next, we are ready to calculate the IM3 collector voltage at frequency 
2(0^ — co 2 . The third-order analysis uses exactly the same equivalent circuit 
as the second-order analysis, and the collector voltage is given by (4.21) 

F C3 (2co 1 - co 2 ) (4.21) 

^BEC^Wi - a> 2 ) • ( Z NL3GPI + Z NL3CPI + Z NL3GPIT + z NL3CPIt) 

+ ^CEC^®1 - ®2^ ' ( Z NL3GM + z NL3GMt) 

+ TF CBC (2(O l - (0 2 ) • (t NL3CBC + Z NL 3 CBCT ) 
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where the distortion currents are now calculated for tone 2 a)j-C 02 , and the 
transfer functions TF X Y c (also calculated at that frequency) are used to 
convert the nonlinear currents / NL from port X-Y to collector voltage. 
However, the equations of the nonlinear current sources will be messier 
than in the second-order analysis, because IM3 currents are not only caused 
by cubic nonlinearities, but by cascaded second-degree nonlinearities, as 
well. According to Table 2.6, both the envelope voltage and the second 
harmonic appear in the IM3 current caused by the nonlinear I B -V BE , for 
example, 

z NL 3 GPI^ W l - © 2 ) = 4 ^ 3 GPI Fg £-(©|) 2 ^ BE ^ 2 ) ( 4 . 22 ) 

+ ^2GPI^ BE^l - W l) + ^2GPI Fgj^fflj) K sjs .(C0 2 ) 



Similar responses are caused for example by the coefficient pairs 
^2GM> ^3GM> and K 2GO , K 3G q of the I C -V BE -V CE nonlinearity. In addition 
to these, the input-output cross-terms X 2GMG0 , X 3GM2GO , and ^ 3GM G02 
cause the following additional terms to Jnl3Gm(2®i-®2) 

1/ 2 • X 2GMGQ • [F BE (to 2 )F CE (2(o 1 ) + F BE (co 1 )F CE (co 2 -( 0 1 ) (4.23) 
+ F BE (w 2 - ®i)^ce( co i) + ^be(2Wi) F ce (w 2 )] 

+ 1 /4 • X 3G M2GO 

• [2f BE ((0 2 )F BE (0) 1 )f CE (w 1 ) + K be (w 1 ) 2 F ce (w 2 )] 

+ 1/4 • ^ 3GMG0 2 

• [2F BE (c0 1 )F CE (C0 2 )F C£ (C0i) + *W^T)F ce ( CO!) 2 ] 



where all such combinations of linear and second-order input and output 
voltages are shown that their products make a tone at 2co 1 -(0 2 . Note that all 
possible permutations are needed in the cross-terms. For example, the last 
term in (4.23) consists of one V BB and two V CB voltages at frequencies ccq, 
00;, and -0) 2 , and they can be arranged in the following three combinations: 
(Fbe(Wi), V G e(©i)> F CE (-© 2 )), (Fbe(©i)> F ce( _ ®2)’ F ge (©i)), and 
(Fbe( _ ®2)’ Fce(©i), F G e(©i)), the first two giving the same result. 
Numerical spectral convolution gives all these combinations automatically. 

Now, IM3 is described in terms of fundamental and second-order node 
voltages. Next, we take one step further and write the distortion as a 
function of fundamental input voltages only. This complicates the 
expressions quite a lot but has the major benefit that it maintains the 
information of the origin of second-order distortion voltages. 
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4.3. 2.4 IM3 Shown as a Function of Fundamental Input Tones Only 

To present the distortion as a function of input tones only, we need to do 
two things. First, F ce (/CQ) is written as TF(j(0)V BE (j(0), where TF is given 
in (4.16) and in Appendix C. Second, the second-order tones must also be 
written in terms of the fundamental inputs, which results in altogether 42 
cascaded second-degree terms, all listed in Appendix C. 

The purely cubic portion of (4.21) is written in (4.24) in terms of V BE 
only. Here, the nonlinear currents are written according to Table 2.6 as 
functions of V BE (j(0) and TF(j(0)V BE (j(0). For example, g 0 in the first term 
in (4.24) is controlled by the product F CE (col) 2 F CE (-co2). Input 
nonlinearities g m , g pi , and C pi are controlled by V BE values only, and the 
cross-terms are controlled both with input and output voltages. For C BC , 
the voltage V CB must be expressed as F CE -F BE =1 (TF-\)V be . 

^C3CUBIc( 2o) 1 - (4.24) 

= {TF cec (2co 1 - co 2 ) • (3/4 • K 3GO ■ TFiaJ 2 ■ TF«o 2 )) 

+ - co 2 ) • (3/4 • 7f 3GM ) 

+ 7 t F’ BE( -.(2(0j — (O 2 ) ' (3/4 • 7f 3Gpi ) 

+ TF BE( -.(2(0j — (O 2 ) • (y'(2c0j — 0 ) 2 ) ■ 3/ 4 • ^ 3GP j) 

+ TF G bc( 2 ®i — 0 ) 2 ) ■ (y(2co 1 — 0 ) 2 ) ■ 3/ 4 • ^ 3GBG 

• (//•(CD,) - l) 2 • (7 /'(C0 2 ) - 1 )) 

+ TFq^q{2(H^ — W 2 ) ' (1/4 • ^f 3G M2GO (277^(0) j) + TF( C0 2 ))) 

+ ^CEC^®1 _ W 2) ' (1/4 • ^ 3GMG 02 ' TF{ COj) 

• (mcoj) + 2TF(C0 2 )))} • ^be( w i) 2 • ^be( m 2) 

The equations for IM3 caused by cascaded quadratic electrical 
nonlinearities become quite complicated, and a full set of equations, 
including 21 mixing products from the envelope and 21 from the second 
harmonic, is given in Appendix C. As an example, a double mixing caused 
by ^2GM*^2GPi I s explained here. This contribution is written as 

*W 2 ®1 - ®2) = ^2GM ' K 2 GPI ' [ rF BEc( 2 «l " ^ 

■ T F CEBE ((d 2 ~ Wj) + TFq EG { 2 co j — W 2 ) ' 2^bebe(® 2 — ®j)] 

• 1^be( C0 i) 2 ' ^Be( M 2) 



(4.25) 
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This equation includes two low-frequency mixing mechanisms, which are 
sketched graphically in Figure 4.9. In the first mechanism, illustrated in 
Figure 4.9(a), a product of F BE (C0i) and F BE (-CQ 2 ) in ^ 2 GM generates the 
envelope current 0 ) 1-002 between the collector and the emitter, from where 
it is then transferred back to the base-emitter voltage by TF CEBE . These 
envelope and fundamental voltages F BE (oo 1 -co 2 ) and F be (G)|) are further 
mixed in generating IM3 current at 2co 1 -(D 2 between the base and 

the emitter. This is further converted to the collector voltage in the transfer 
function TF BEC . The second mechanism, shown by Figure 4.9(b), can be 
explained as follows: F be (G)|) and F BE (-CQ 2 ) generate the envelope current 
CO 2 — 0) 2 directly at the base in /f 2 GPi- This current between the base and the 
emitter is transferred to the base-emitter voltage by TF BEBE , and the 
resulting envelope voltage mixes with the fundamental base-emitter voltage 




Figure 4.9 IM3L caused by the cascaded second-degree distortion mechanisms 
AT 2 gpi and K 2 q M via the envelope frequency. From [11], 
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in ^2GM> as a result of which an IM3 current is generated between the 
collector and the emitter. This current is finally converted to the collector 
voltage by TF CEC . 

Such multiple mixing products are quite common. For example, one 
major cause of IM3 in deep class AB or class B amplifiers is the following: 
When clipping asymmetrically, the transistor causes a high second 
harmonic at the output (this is modeled by high K 2GM ). As the frequency of 
the second harmonic is high, it couples easily through C BC back to the 
input and mixes with the fundamental again in K 2 qm, causing IM3 current 
directly in the output. This mechanism (named as term V cm in Appendix 
C) can be minimized by attenuating the second harmonic either at the 
collector or at the base. 

4.3. 2. 5 Electrothermal Terms 

Finally, the third-order electrothermal distortion mechanisms are described. 
These consist of third-degree terms TF XYZ *(K iXXT *v*tj) and cascaded 
second-degree terms, where the thermally induced second-order distortion 
mixes with fundamental tones in electrical square-law nonlinearities. The 
third-order terms can be expressed by 

y C3T3^®1 — ®2^ (4-26) 

= TF CEC (2(d 1 - 0 ) 2 ) ' (^3 got ' T’-f (fflj) ■ — ® 2 )) 

+ 7 ^CF.C^®! _ W 2) ' ^3GMT ' ^Be(®i) ' ^j(®l - ®2) 

+ ^BEC^Wi ~ W 2 ) • ^ 3GPIT • ^BE^l) ' - W 2) 

+ j( 2®! ~ 0) 2 )^3CPIT ' ^BE C (2C0i - W 2 ) ' ' ^j(®l - ®2^ 

+ j( 20)! - w 2 )i: 3CBCT • TF CBC (2(o l - cd 2 ) 

• rmcoj)- !] • ■ Tj(co 1 -co 2 ) 

where Tj is the junction temperature, TF is the Fce/^be rat io, an d TF X YC is 
the transimpedance from the nonlinear source to collector. Note that the 
pattern is always the same: low-frequency temperature variations modulate 
the value of the circuit element, which causes mixing with the fundamental 
tone. The junction temperature is calculated based on the power dissipation 
Pdiss = V CE ; C’ which results in both first- and second-order tones. Only the 
envelope tone Pdiss(®2-«>l) = (F CE (® 2 ) / c(-®l3 + V C t(-(£>\)I c (<£> 2 ))a is 
picked and multiplied by the thermal impedance Z TH (a) 2 _ft >i) to obtain the 
variation in the junction temperature. The dissipating source is always the 
same, but if different temperatures are needed for different circuit 
elements, the thermal network can be modeled by a multiport Z-matrix. 
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The IM3 contributions of cascaded second-degree nonlinearities are 
again too complicated to be presented here, but these are shown in 
Appendix C, and their numerical values are plotted in the vector plots in 
Sections 4.4.2 and 4.5. Note that the terms including a temperature- 
dependent charge V 2CX xt cannot be calculated unless we find a way to 
measure or derive a value for it. 

As seen from the above and Appendix C, already the symbolic Volterra 
analysis of a simple CE amplifier gets quite involved, and more complex 
topologies are too difficult to handle analytically. Volterra calculations can 
be performed numerically as well, however. The Nlsim software [24, 29] is 
capable of calculating the nonlinear voltage components numerically, and 
more complex topologies can be analyzed. Compared to the fully analytical 
solution, Nlsim gives a slightly less detailed picture of distortion, as it does 
not separate the cascaded quadratic nonlinearities but simply displays the 
cubic term and up- and downconverted envelope and second harmonic 
terms for each nonlinearity. However, the second-order phasors can be 
plotted as vector sums to see the dominant second-order contributions. 



4.4 IM3 in a BJT CE Amplifier 

Here, IM3 distortion is studied in two different cases: first using a 
simplified cascade model of nonlinear input impedance and I-V curve, and 
then the full circuit, including also the feedback effects, employing the full 
analysis presented in the previous section. 

4.4.1 BJT as a Cascade of Two Nonlinear Blocks 

The distortion composition of a BJT amplifier is quite complicated, as 
noted in the previous section, and therefore some simplifications are first 
made here to find some of the basic characteristics of distortion 
composition. A complete analysis will be given later, but for the moment 
the feedback effects of the emitter impedance Z E and feedback capacitance 
C B c are ignored and the collector current is simply considered to be one- 
dimensional, a function of base voltage only. The simplified circuit is 
shown in Figure 4.10 and since all feedback effects have been ignored, the 
collector current of the BJT can now be regarded simply as a cascade of 
two nonlinear blocks, as shown in Figure 3.5. The first block contains the 
input nonlinearities C pj and g pi and the second the nonlinearity of the one- 
dimensional transconductance. 

It is commonly known that the linearity properties of the CE BJT 
amplifier are different when using either voltage or current excitation. In 
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Figure 4.10 The simplified BJT model, ignoring the feedback effects. © IEEE 2000 
[30], 

the idealized case the voltage excitation with zero Z IN makes the input 
nonlinearities negligible, because the input and base-emitter voltages are 
equal, and the voltage source shorts all distortion currents. Then the only 
effective nonlinear element in Figure 4.10 is the transconductance. The 
situation changes if a current source excitation is applied. Figure 4.11(a) 
shows the base-emitter voltage as a function of input current and Figure 
4.11(b) the collector current as a function of the base-emitter voltage, 
assuming purely exponential g pi and g m . Figure 4.11(a) can be recognized 
as a logarithmic function, while Figure 4.11(b) is an exponential one, and 
as a cascade of them, the collector current as a function of input current is 
perfectly linear, as shown in Figure 4.11(c). Thus, the strong exponential 
nonlinearity of a voltage-driven BJT disappears, when the transistor is 
driven by a current. Also, in practice the type of the excitation can be 
adjusted somewhat by the value of Z IN : the higher the value of Z ]N , the 
more the CE BJT appears to be current driven. 

In Figure 4.11 the two opposite nonlinearities cancel one another. The 
same phenomenon is also demonstrated in the Gummel plot in Figure 4.12, 
in which the collector and base currents as functions of F BE are presented. 
Since a logarithmic v-axis is used, these purely exponential nonlinearities 
appear as straight lines. It can be observed that the distance (i.e., the 
current gain) between the curves is constant and independent of the value 
of Z BE , indicating that the shapes of the nonlinearities are similar. These 
two nonlinearities canceling each other out are often referred to as tracking 
nonlinearities, resulting in linear operation with nonlinear elements. In 
practice, however, the nonlinearities do not track completely, as a result of 
which some amount of nonlinearity always exists. For example, the dashed 
parts of the lines in Figure 4.12 represent the practical situation in which 
the transconductance is not purely exponential at high F BE values anymore 
due to a high injection level at the base, thus reducing the current gain. 
Also, at very low Z BE values some leakage base current always exists. 
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Figure 4.11 Illustration of tracking nonlinearities: (a) F be as a function of 7 in , (b) 7 C 
as a function of *be, and (c) Ic as a function of 7 in . 



Next, the ac behavior of the circuit shown in Figure 4.10 is analyzed. 
The nonlinearity of the first block is caused by g pi and C^; g pi is considered 
to be purely exponential, and C pi is exponential, too, being calculated from 
the forward transit time (t f ) and transconductance g m as follows [31]: 

C pi = W m . (4-27) 

The latter block represents the nonlinear transconductance, and its 
nonlinearity can be calculated from (4.1). If we ignore the Early effect and 
use only one-dimensional collector current, the nonlinearity coefficients 
listed in Table 4.2 can be obtained. The second column corresponds to the 
strictly exponential low injection case, and the fourth column presents a 
high injection condition, in which the nonlinearity of the transconductance 
is reduced, corresponding to the dashed line in Figure 4.12. 




Figure 4.12 Base and collector current as a function of base-emitter voltage. Dashed 
parts of the curves present high injection and leakage current effects. 
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Table 4.2 

Nonlinearity Coefficients for Transconductance (From [30]) 



Degree of 
nonlinearity 


1 C v 1 KF 


i c = io - 1 ■ I KF 


I C = iO" 1 ' I KF 


1 


I C /V T 


0.981 • I C /V T 


0.845 • I C /V T 


2 


I c / (2 • V 2 t ) 


0.962 • I c /( 2 ■ V 2 r ) 0.724 • 7 c /(2 • V 2 T ) 


3 


I c / (6 ' V\) 


0.925 • I c /( 6 • v\) 0.535 • 7 c /(6 • v\) 



The distortion composition of cascaded nonlinearities was discussed in 
Chapter 3 and is shown graphically in Figure 3.6. The IM3L collector 
current of the cascade is found to be 

= ^Be(2®1 - ®2) • g m 

+ 3/4 . ^(©j) • V BE^) • ^ 3GM 28 

+ ' ^Be(2®i) • ^2GM 

+ ^BE(®l) ' ^Be(®1 _ ®2) ' ^2GM 

IM3L consists then of four components, the first two of which are 
generated directly inside the first and the second block, respectively. The 
last two are related to second-order interaction via the envelope and second 
harmonic frequencies between the blocks. In the first block the distortion is 
generated by the second-degree nonlinearity of the input impedance and the 
spectral components at the base for (4.28) can be calculated using Table 
2.6. 

The effects of the fundamental Z IN of the circuit in Figure 4.10 are 
studied first. The nonlinearity coefficients are calculated using equations 
(4.1) and (4.27), based on the GP model SPICE parameters taken from a 
BFG11 RF power BJT [32], The IM3 levels are calculated for a center 
frequency of 1.8 GHz and tone spacing of 1 MHz. If the value of the input 
impedance Z IN is high compared to the internal base impedance, the input 
will behave like a current source and the input-output nonlinearities will 
partially cancel each other out. 

Z 1N at the fundamental is swept over the range of reasonable values to 
check the effects of voltage/current excitation. The input of the BJT is 
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conjugately matched, which means that the input reactance was chosen so 
that the base capacitance is tuned out (/*imag(Z rN ) = +/0.6). z m values at 
the harmonic and envelope are set to infinity, which means that the base 
impedance at these frequencies consists of the internal base-emitter 
impedance only. IM3 distortion components at two values of the collector 
current are shown as functions of fundamental Z IN in Figure 4.13. The IM3 
asymmetry (i.e., the vector error between the lower and upper sidebands) is 
also given in Figure 4.13. 

Under low injection conditions, IM3 improves by increasing the value 
of Z IN up to 20 to 30 ohms. Increasing the value of Z IN any further does not 
reduce IM3 any more, because incomplete tracking between the input and 
output nonlinearities limits IM3 to a level of -57 dBc, as seen in Figure 
4.13. When approaching high injection conditions (7 C = 0. 1 IKF), the 
situation is different: although the linearity is quite independent of Z IN at 
very high values of the latter, a significant linearity improvement can be 
achieved at some relatively small values. The higher I c is compared to the 
knee current IKF, the better is the linearity improvement achievable at 
optimum Z IN . 

When the imaginary part of the fundamental Z IN is increased from 
conditions of conjugate match, the linearity starts to deteriorate, but even 
though it may be reduced slightly, the vector difference (i.e., the 




Figure 4.13 Amplitude of the IM3 sidebands of the collector current as a function of 
the fundamental input impedance. © IEEE 2000 [30], 
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asymmetry) between the IM3 sidebands may reduce. The lower curves in 
Figure 4.13 represent the asymmetry between the lower and upper IM3 
tones. By making the input match slightly inductive, a significant 
improvement of more than 10 dB may be achieved in the symmetry of the 
IM3 sidebands. Similar observations have been reported in [33], where the 
effects of input impedance on phase distortion were studied. The 
improvement in the symmetry can be a very important advance in 
applications involving linearization. 

Since a large amount of the output third-order distortion is generated 
via cascaded second-degree nonlinearities, careful optimization of out-of- 
band Z IN can improve the linearity. Third-order terms cannot be affected 
much by filtering, because their power overlaps the fundamental signal, but 
the power of the second-order signal lies well away from the fundamental, 
and filtering can be used to achieve the best possible linearity [27, 34-36], 
Also, there are quite a lot of requirements for the fundamental matching 
such as gain and efficiency, so that the fundamental impedance cannot be 
chosen based only on the linearity properties. However, the out-of-band 
matching has only secondary effects on gain and efficiency, and therefore it 
can sometimes be tuned mainly based on the linearity. 

We will now look at the effects of second-order signal components. 
The impedance around the fundamental is kept constant while Z IN at the 
envelope is swept. The upper eight-segment vector in Figure 4.14(a) 
represents the output IM3 product caused by input nonlinearities (C pi and 
gpi), while the two first parts of the lower three-segment vector represent 
the output third-order IM3 products caused by second-degree input-output 
nonlinearities. The first of these is generated via the envelope at the base 
and the other one via the second harmonic. The third part of the lower 
three-segment vector represents the cubic nonlinearity of the 
transconductance. It is interesting to note that the parts of the IM3 vector 
caused by the quadratic nonlinearities are opposite to the part caused by 
cubic nonlinearities. This means that the nonlinearities track each other in 
two different ways: First, through the input-output tracking explained 
earlier, which means that the distorted voltage waveform at the base 
cancels the distortion caused by nonlinear transconductance, and second, 
through the cancellation between the second- and third-degree 
nonlinearities. The latter may be difficult to note from the equations 
presented in this book due to extensive use of transfer functions. However, 
if we ignore C pi and write the expression for the third harmonic current in 
terms of circuit elements and coefficients in (4.29) [2, equation 8.76], the 
opposite signs of the second- and third-degree coefficients are clearly 
visible. 
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Figure 4.14 Representation of (a) partially tracking nonlinearities, and (b) the result 
of vectors. Both output IM3 signals consist of 1 1 contributors: K2GPIe, 
K2GPIh, K2CPIe, K2CPIh, K2K2GPICPIh, K2K2GPICPIe, K3GPI, 
K3CPI, K2GMe, K2GMh, and K3GM (e and h mean envelope and 
second harmonic, respectively). © IEEE 2000 [30], 



/out( 3o) ) = (“gpi ' ^3 GM + S m g pl K i GPI (4.29) 

+ 2g p i^2GM^2GPI - 2£m^2GPl) 

The result of the IM3 caused by these partially tracking nonlinearities 
is illustrated in Figure 4.14(b) for four values of the impedance at the 
envelope. If the magnitude of that impedance is high, the output IM3 signal 
will not be affected by its phase, but if it is relatively small, say at most one 
decade higher than the Z IN at the fundamental frequency, its phase will play 
an important role. The reason why Z IN at the envelope has an impact on 
IM3 only when it is small is simple: If Z IN is high, the total impedance of 
the base node is dominated by the internal impedance of the transistor, but 
when Z IN is small, it affects the total node impedance. 
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Minimizing the low frequency Z IN has often been suggested as a way to 
reduce overall distortion. As seen above, it is not always the case; instead, a 
small imaginary value for Z IN may reduce the total distortion, or improve 
the balance between the IM3 sidebands. Since distortion currents 
multiplied by Z [N (C 02 -OJ|) and Z IN (co 1 -( 02 ) mix to different sidebands and 
the impedances have opposite phases [as Z IN (co 1 -co 2 ) = Z IN (co 2 -(0 1 )], a 
small and reactive baseband bias impedance may be quite handy in 
correcting phase errors between the upper and lower IM3 sidebands. 

4.4.2 Detailed BJT Analysis 

Now the full Volterra model is employed to study the effects of individual 
distortion mechanisms. In this case the nonlinearity coefficients of the 
model are extracted by measurements, using a procedure that will be 
explained in detail in Chapter 5. In this section, the effects of the bias point 
and optimum terminal impedances at different frequencies on the linearity 
of the CE BJT amplifier are discussed. 

It is a well-known fact that the optimum impedances of a CE amplifier 
vary with the biasing conditions. It is impractical to examine 
experimentally all possible combinations of fundamental input and load 
impedances at different bias voltages, and if the envelope and harmonic 
impedances, which have a great impact on linearity, were to be taken into 
account, too, optimization of a single CE stage would become a very 
lengthy task. The Volterra model - provided it is accurate enough - can be 
used for optimizing the linearity of a CE amplifier by designing optimal 
input and load impedances under different sets of biasing conditions. 

An amplifier based on a BFG11 transistor [32] is constructed and the 
input impedance at a fundamental frequency of 1.8 GHz is tuned to the 
conjugate match for maximum power transfer. The load impedance is 
chosen so that the imaginary parts of the output reactance and matching 
network canceled each other out, and the real part of the Z L is determined 
by the desired I-V characteristics. Linearity can be improved by lowering 
the value of R L , because the voltage swing at the collector decreases, but 
unfortunately, the efficiency decreases at the same time. So R L was chosen 
to be 20 ohms as a trade-off between linearity and efficiency. The load 
impedance around zero frequency is small, and has to be so to supply the dc 
energy from the supply to the collector with minimum losses. The input 
and load impedances of the amplifier are 0.1+/76 and 31+/17 ohms, 
respectively, at the second harmonic and 13^/0.2 and 0.2-/1 .4 at the 
envelope frequency of 2 MHz. These figures are taken from measurements 
performed on the actual, implemented amplifier. 
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The base and collector bias voltages are first swept and the linearity 
monitored to find some basic characteristics of the biasing. To make the 
simulations comparable to each other, the input voltage swing and load 
resistance are tuned a little to keep the output voltage and current swings 
constant under varying biasing conditions. An output voltage swing of 
1.5 Vpp is chosen, which means that the IM3 values near the collector 
supply of as low as IV are not accurate. Calculated IM3L contours and 
measured IM3L points, in dBc, are shown in Figure 4.15(a), and a good 
correlation between these two is observed. 

Maximum linearity is achieved at a bias current of 75 mA, and linearity 
deteriorates with a decreasing bias current. The IM3 sideband increases by 
approximately 10 dB when the bias current is lowered to 40 mA. The 
collector bias voltage affects the linearity as well, so that 3.5V gives the 
maximum linearity. Collector voltages beyond 4V are not shown in Figure 
4.15(a), because higher collector voltages only reduce the efficiency 
without giving further linearity improvement. At low currents, the linearity 
is quite independent of collector voltage, but at high current values the 
collector voltage plays a more important role. This is obvious because at 




Figure 4.15 (a) Calculated and measured (dots) IM3L and (b) calculated asymmetry 

IM3L-IM3H as functions of the biasing condition for an output voltage 
swing of 1.5 Vpp. From [1 1], 
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lower collector current values the nonlinearity caused by g pj , C pi , and g m is 
much stronger than that caused by g 0 and cross-terms. Moreover, at high 
current values the nonlinearity caused by g p p C pi , and g m is weaker, 
because transconductance, for example, linearizes with an increasing 
current due to the effects of high injection. Also, since g 0 and its 
nonlinearity are quite strong at high collector current values, g 0 takes a 
more dominant role with respect to total distortion. 

The contour plot for the IM3H sideband is quite similar to that for the 
IM3L in Figure 4.15(a). The vector difference between the sidebands is 
plotted in dBc compared to the fundamental signal in Figure 4.15(b). It is 
important to emphasize that this figure is closely dependent on modulation 
frequency and not only on bias values and matching impedances. The 
asymmetry between the sidebands is an important figure of merit if 
predistortion linearization is employed, as discussed in Chapter 3 . 

The values for the fundamental input and load impedances in the 
previous example are chosen to achieve the desired power and gain 
characteristics with reasonable efficiency. For the best possible linearity, 
however, the out-of-band envelope and harmonic terminations have to be 
optimized. One commonly used approach for improving the linearity by 
tuning the out-of-band terminations is to minimize the load impedance at 
the second harmonics. Since the IM3 components are partially caused by 
the second harmonic voltages, the impedance at that frequency affects IM3 . 
Whether or not this improves the linearity depends on the nonlinearity 
coefficients and other impedances. To check the effects of the second 
harmonic matching, the real part of the load impedance at that frequency is 
reduced to 6 ohms. The result of the comparison is plotted in Figure 4.16, 
which shows improvements of 1 dB to 5 dB under all biasing conditions. 
These calculations show that optimization of out-of-band impedances is 
needed to achieve the best performance. Only the effects of the second 
harmonic load impedances were demonstrated, but other out-of-band 
effects are important, too. Carefully selected optimum out-of-band 
terminating impedances can improve the linearity by several decibels 
without reducing the power and efficiency performance at the fundamental. 

The asymmetry between the IM3 sidebands is greatly affected by the 
envelope impedances, as noted in Chapter 3 . In multicarrier transmitters the 
bandwidth of the signal may be very wide, and it is very hard to design 
constant impedances between dc and 20 MHz, for example. The input and 
load envelope impedances are usually not constant, and since the IM3 
sidebands are functions of these impedances, they vary as a function of the 
modulation frequency, causing memory effects. 

The memory effects are investigated by sweeping the tone spacing of a 
two-tone signal at collector and base bias voltages of 3V and 740 mV, 
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Figure 4.16 Calculated IM3L at different load impedances at the second harmonic 
(thin lines 35 ohms, thick lines 6 ohms). From [11], 



respectively. The phase of the Volterra-calculated IM3L is plotted in Figure 
4.17 as a function of tone spacing, and three different types of memory 
effects are met: thermal memory effects at low frequencies, resonance-type 
memory effect, and high-frequency memory effects, both caused by the 
biasing circuits. The resonance at 1 MHz is caused by a resonating 
collector impedance. It cannot be observed with harmonic balance (HB) 
simulations using the Gummel-Poon (GP) model, because its 
oversimplified output impedance masks the effect of the collector 
resonance. The smooth phase deviation at high modulation frequencies is 
caused by the input impedance at the envelope frequency, and it is 
nevertheless simulated correctly by the GP model, too. Since the dynamic 
self-heating effects are not implemented in the basic GP model, it naturally 
cannot predict the low frequency memory effects caused by thermal effects. 
The Volterra simulations agree reasonably well with the measured results, 
given also in Figure 4.17 and explained in more detail in Chapter 6. The 
Volterra model therefore seems to be a good tool for recognizing memory 
effects in a power amplifier. 

Next we will look at the fine structure of the IM3 phasors. This 
information can be used in design optimization, and IM3 is drawn here as a 
vector sum of tens of contributions to study why the phase of the IM3L 





106 



Distortion in RF Power Amplifiers 




Figure 4.17 Calculated, simulated, and measured phases of total IM3L as a function 
of tone spacing. From [11], 



depends on modulation frequency. Three frequencies are chosen for the 
vector plots, 100 kHz, 1 MHz, and 1.1 MHz, and the resonance at 1 MHz 
and the thermal memory effects are studied. 

We first look what happens just above the resonance at tone spacing of 
1.1 MHz. The IM3L phasor is plotted in Figure 4.18(a) as a vector sum on a 
real-imaginary axis, starting from the origin at the top center. The vector 
consists of 49 purely electrical terms, the first seven of which are caused by 
cubic nonlinearities and the following 42 are generated by cascaded 
second-degree nonlinearities via the envelope and second harmonic 
frequencies. The electrothermal terms are practically zero, because the 1.1 
MHz beat frequency already lies in the stopband of the thermal filter. 

The first seven points (1-7) are the cubic distortion mechanisms, of 
which Wjcpj (6) is the largest. The following 21 (8-28) are upconverted 
IM3 components from the envelope frequency, and finally, the last 21 (29- 
49) are downconverted from the second harmonic. One interesting finding 
is that there is not a single dominant contribution that we can attack; 
instead, the total IM3 is already smaller than any of the dominant 
contributions. This is due to several pairs that partially cancel each other: 





lmag(IM3L) [mV] _ ^ _ lmag(IM3L) [mV] 
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Figure 4.18 Vector representations of (a) IM3L@l.lMHz, and (b) the vector 
difference between IM3L@l.lMHz and IM3L@1.0MHz. From [11], 
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the third-degree nonlinearities AT 3GM (2) and K^q P] (6) get partially 
canceled, and so do the second harmonic mixing results ^ 2 GM*^ 2 GM (29, 
where the second harmonic caused by g m returns to the input and mixes 
again), and K 2 c P * K 2 q M (37, second harmonic caused by C pi mixing in 
^2 Gm)> or ^2GM*^2CBC (45, second harmonic caused by g m mixing in 
C BC ) and A' 2 cbc*^ 2CPI (48, second harmonic caused by C pi being 
amplified and mixed in C BC ). 

What could be done to improve linearity? Distortion proportional to 
Wjqpi is converted to voltage in fundamental Z IN> and its contribution (6) 
can be rotated towards the origin by adjusting the phase of Z IN . Further, the 
size of the entire 29-49 mesh can be reduced by lowering the second 
harmonic impedance. On the other hand, a 5% to 10% increase of the base 
impedance at the second harmonic would increase terms 37 and 48 and 
force the total sum closer to zero. 

From the memory effect point of view it is instructive to study what 
happens when the tone spacing is reduced to 1 MHz, hitting exactly the 
resonance in Z L . This is illustrated in Figure 4.18(b), which depicts the 
difference of IM3L@l.lMHz and IM3L@1.0MHz as a similar vector sum. 
First, the cubic nonlinearities (l)-(7) behave in exactly the same way, as the 
phases of the fundamental tones do not vary. The second harmonic mixing 
products (29-49) are also equal, indicating flat terminal impedances at the 
second harmonic. All of the dominant causes of the IM3 resonance arise 
from cascaded quadratic nonlinearities that are upconverting the envelope 
frequency - large terms include ^2GMGO*^2GM (14, output envelope 
mixing in the ^2GMGo)’ ^2GM*^2CBC (24, output envelope mixing in 
C BC ), and K 2 qpi*K 2 cbc (26, the input envelope being amplified and mixed 
in C BC ), that all involve the baseband frequency response of Z L . 

Memory effects generated by interaction between electrical and 
thermal behavior are studied next at a narrow tone spacing of 100 kHz. 
Purely electrical IM3 vectors are presented by the first segments of the 
vectors in Figure 4.19, and at such a low frequency they are almost 
identical. The electrothermal 30-segment vector is drawn in a termwise 
manner, and the electrothermal distortion is seen to be dominated by large 
^3 GMT’ ^2GMT’ ^3CPIT’ ar| d ^2CPIT terms, which after all cancellations at 
the end point 78 cause only a 15% difference to IM3 amplitudes caused by 
purely electrical mechanisms. Note, however, the big phase difference 
between IM3L and IM3H components. It is again due to the fact that 
baseband effects, including the thermal feedback, mix with opposite phases 
to lower and higher IM3 sidebands. Besides phase asymmetry, the thermal 
feedback also causes here amplitude asymmetry between IM3L and IM3H, 
as the distance between the origin and the final point 78 is different for 
IM3L and IM3H. 
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Figure 4.19 Vector representation of electrothermal IM3L (dashed) and IM3H (solid) 
distortion. Tone spacing is only 100 kHz to see the thermal effects. From 
[ 11 ]- 



4.5 MESFET Model and Analysis 

A Volterra model for a MESFET common source amplifier is derived in 
this section. Its small-signal equivalent circuit, given in Figure 4.20, 
includes an input impedance Z IN (consisting of the driver impedance, 
external input impedance, and internal series gate impedance), a gate-to- 
source capacitance (C GS ), a feedback capacitance (C GD ), a drain-to-source 
capacitance (C DS ) and a resistance (l/g 0 ), transconductance (g m ), a load 
impedance (Z L ), and a source impedance (Z s ). As in the BJT case, the input 
and load impedances include not only impedances of the matching 
networks, but also impedances of the bias networks. The drain current of 
the FET is modeled as a three-dimensional function of the gate and drain 
voltages and temperature similar to (4.3) for a BJT. C GS and C GD are also 
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Figure 4.20 Linearized first-order circuit for a common-source FET amplifier. 



regarded as nonlinear, being functions of the gate-to-source voltage and 
temperature and of the drain-to-gate voltage and temperature, respectively. 

The linearized circuit for a common-source FET amplifier presented in 
Figure 4.20 is pretty much the same as that for a BJT amplifier shown in 
Figure 4.7. By neglecting the g pi and changing the names of the circuit 
elements, the equations deduced in Section 4.3.2 can be used for the FET 
analysis as well. Similarly, by omitting g pi and its nonlinear current source, 
a circuit containing a distortion source for an FET can be obtained, as given 
in Figure 4.21. The nonlinearities of the device are characterized by 22 
first-, second-, and third-degree nonlinearity coefficients, of which 15 are 
purely electrical (C GS , A 2CGS , ^3CGS> QjD’ ^2CGD> ^3CGD> gm> ^2GM> 
^3GVb §0’ K 2GO’ ^3 GO’ ^2GMGO’ ^3GM2GO’ and ^3GMG02) and seven are 
related to temperature variations C^2CGST’ ^LtCGST' ^ 2 (( i DT* ^3CGDT’ 
^2GMT’ ^3 GMT’ and ^3 got)- Equations (4.3), (4.6), (4.24), (4.26), and 
(4.28) as well as Appendix C can be used for an FET simply by ignoring g pi 
and changing the voltages and subscripts BE to GS and BC to GD. 

Circuit elements and nonlinearity coefficients for an Infineon CLY2 
GaAs MESFET [37] are extracted using an ^-parameter characterization 
method discussed in Chapter 5. The input and load impedances measured 
from an implemented amplifier are listed in Table 4.3. 
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Figure 4.21 Circuit containing second-order distortion sources. The thermal ones can 
be obtained by connecting (nl2CGST ' n parallel with /’nl2CGS an d 

'NL2GMT with 'NL2GM- 



Table 4.3 

Input and Output Matching Impedances at Different Frequency Bands 



Frequency band 


Z IN 


Z L 


2 MHz (envelope) 


1.3k -/1.3k 


0.5 1 — jA.% 


1.8 GHz (fundamental) 


6.6+/11 


26 — y'6.6 


3.6 GHz (2nd harmonic) 


1 8 —747 


3.6— y'8.7 



The IM3 calculations over the range of bias current and voltage at the 
center and modulation frequencies of 1.8 GHz and 2 MHz and with an 
output voltage swing of 2 Vpp are depicted by a solid line in Figure 4.22. 
Large linearity variations of more than 30 dB are observed over the I-V 
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plane, but since the dc power consumption changes with bias point 
according to the linearity, no significant advance in terms of the efficiency- 
linearity trade-off can be achieved. The Volterra calculations agree 
reasonably well with the measured linearity values represented by the dots 
in Figure 4.22. 

Next, the effects of the out-of-band impedances are studied. The value 
of the input impedance at the envelope frequency is increased from 
1.3k-yl.3k to 10k - j \ Ok to see the effect of bias impedances. The new 
linearity contours drawn with dashed lines in Figure 4.22 show some 
significant changes. The linearity is reduced by a few decibels at low drain 
voltages, but very good linearity improvements of more than 10 dB are 
observed at high voltages of 4V and a drain current value of 150 mA. Since 
Z GG (env) is highly frequency-dependent, however, the improvements 
presented in Figure 4.22 are quite narrowband. 

The changes caused by the modified gate impedance at the envelope 
frequency are demonstrated in Figure 4.23, where IM3 vectors at 2-MHz 




Z<3e(2MHz) = lOk-jlOk 

Figure 4.22 Calculated (lines) and measured (dots) IM3L contours in dBc at two 
values of the input impedance at the envelope frequency. Center and 
modulation frequencies are 1.8 GHz and 2 MHz, and output voltage 
swing is 2 Vpp. © IEEE 2002 [38], 
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tone spacing for the two baseband input impedances are calculated at the 
bias point of V D =4.5V and I D =150 mA, shown in Figure 4.22. The total 
distortion is dominated by the cubic distortion mechanisms /f 3GM , 
^3GM2GO ^3GMG02> ar| d ^3CGS’ an( 3 some cancellation between the 
mechanisms can be seen in the figure. The effect of the cross-terms, 
especially ^3 qm 2GO’ i s ver y significant, and most of the large contributions 
are due to the I-V characteristic. However, the reason for improvement at 
the higher Z GG is caused by C GS (term A 2CGS 2 , shown for Z GG =1 Ok— /'I Ok 
only): The larger gate impedance amplifies the envelope current generated 
in ^2CGS> an( 3 after it mixes again in K 2 qgs to IM3, it mostly cancels the 
large K^CGS term, as shown in Figure 4.23. 

The effects of gate impedance at the envelope frequency have been 
investigated also in [39], where linearity changes caused by varying the 
gate impedance have been demonstrated by measurements. The same effect 
has been measured for the baseband drain impedance in [40, 41], Optimum 
gate and drain impedances were found to be low in all of these cases [39- 
41], but it is shown here that the optimum envelope impedance can be 
either high or low, depending on the amplitudes and phase angles of the 
nonlinearity mechanisms. An optimum nonzero envelope impedance was 




Figure 4.23 Vector representations of IM3 components for two gate bias 
impedances. The bias point is Vd=4V, Id=150mV and the tone spacing 
is 2 MHz. © IEEE 2002 [38], 
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also found by measurements in [35], In conclusion, the effects of out-of- 
band impedances can be used for optimizing the distortion behavior, and a 
certain amount of improvement in IM3 performance can be expected. In 
general the effects of out-of-band impedances are often undesirable, 
however, because the envelope impedance (mostly caused by the dc bias 
circuit) cannot be constant over a wide range of modulation frequencies. 
Since the IM3 components are affected by envelope impedances, memory 
effects will occur. 

To study the memory effects, a tone spacing sweep is again simulated 
using the Volterra model. The phases of the IM3L and IM3H tones, plotted 
in Figure 4.24, show a very good agreement with the measured results 
(presented in more detail in Chapter 6). The bump at 500-kHz tone spacing 
is also in this case caused by a resonance in the envelope drain impedance, 
and the phase drift at very high tone spacings is caused by a frequency- 
dependent baseband gate impedance. This 5° to 10° variation in the phase is 
certainly not serious in a standalone amplifier, but it is enough to reduce the 
achievable cancellation below 1 5 dB in a predistorted PA, as will be seen in 
Chapter 6. 




Figure 4.24 Calculated (smooth) and measured (rough) phases of the IM3 
components as functions of tone spacing. © IEEE 2002 [38], 
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4.6 Summary 

The requirements for the simulation models are challenging and since all 
the requirements cannot be fulfilled, the models are optimized for different 
purposes. Most of the models for transistors/amplifiers can be divided into 
two classes: behavioral and device models. Behavioral models just try to 
imitate the measured phenomena without any information on internal 
device operation, while device models more or less imitate the physical 
operation of the device. A lot of research is going on in developing the 
device models, but still they are not completely optimized for RF power 
amplifier design. Instead, most semiconductor models are aimed for analog 
IC design, where a wide range of operating conditions are used, and models 
need to be scalable for different sizes and geometries. From an RF power 
amplifier design point of view, accurate distortion simulations are essential, 
and the derivatives of the I-V and Q-V curves and self-heating effects 
should be as accurate as possible. Unfortunately, only a few manufacturers 
provide parameters for the most sophisticated power transistor models. 

The Volterra analysis is an extension to small-signal analysis, and the 
nonlinearities of the circuit elements are modeled by polynomial functions 
that are extracted around the desired bias voltages. In this way component- 
level information is achieved and the model is fitted locally just around the 
chosen operating point, thereby providing good accuracy in simulations. 
The main advantage of the model and using the Volterra analysis for 
distortion simulation is that elementwise information of IM3 can be 
obtained: IM3 can be drawn as a sum of vectors, each presenting the 
nonlinearity of one circuit element or mixing mechanism. In this way the 
Volterra analysis provides insight into distortion mechanisms, and gives 
information about the dominant ones and possible cancellation schemes. 
This is very helpful when optimizing the matching networks and selecting 
proper bias voltages. 

The third-order Volterra model presented in this chapter can be used for 
both BJT/HBT and MESFET common emitter/source amplifiers, with 
slight modifications. The model includes a nonlinear three-dimensional 
collector/drain current which is a function of the base-emitter/gate-source 
and collector-emitter/drain-source voltages and of temperature. The input 
capacitance C pj /C GS is a nonlinear function of the base-emitter/gate-source 
voltage and temperature. The feedback capacitance C bc /C gd is also a 
nonlinear function of the collector-base/drain-gate voltage and 
temperature. In a BJT model, g pi is a function of the base-emitter voltage 
and temperature, although its effect is small whenever the transistor is used 
at high frequencies, where C pi dominates the input distortion. Since the 
impact of the other circuit elements on total distortion is small, these are 
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considered to be linear components. The model presented in this chapter 
also includes the feedback capacitance C bc /C gd and the emitter/source 
impedance Z E /Z S , which have been neglected in most previously published 
Volterra models. The circuit is solved analytically, and complete equations 
for the IM3 components are presented here and in Appendix C. 

The collector/drain current used in this analysis is a complete Taylor 
expansion of its variables, and compared to most other Volterra models, the 
input-output cross-products if 2G MGO ^3GM2GO ar| d ^3GMG02> which have 
a significant impact on distortion performance, are taken into account. The 
temperature on the surface of the chip is proportional to the instantaneous 
power dissipation v CE z G . Thus the junction temperature is already a second- 
order phenomenon, and it is modeled as an independent, low frequency 
variable. The thermally induced IM3 can be calculated using the Volterra 
model, and it affects IM3 at low modulation frequencies of up to hundreds 
of kilohertz. 

The IM3 at the output is not only the sum of the effects of cubic 
nonlinearities, but also the cascaded quadratic nonlinearities have an 
impact on the total amount of distortion. Therefore, the out-of-band 
impedances at the envelope C0 2 -tflj and second harmonic 2 c 0 | can be used to 
optimize the distortion behavior, and linearity improvements of some 
decibels can be achieved by using optimum out-of-band impedances. Since 
the impedances, and especially the envelope impedance, cannot be constant 
over a wide range of modulation frequencies, the amplitude and/or phase of 
IM3 becomes dependent on the modulation frequency, which is very 
harmful with many linearization techniques. These memory effects can be 
simulated with the Volterra model, which is capable of predicting the 
measured memory effects with sufficient accuracy. 

The impact of out-of-band impedances and electrical memory effects 
seem to be stronger in BJTs than MESFETs. In BJTs, there are two strong, 
almost exponential nonlinearity mechanisms that partially cancel out each 
other. This phenomenon is very sensitive not only to fundamental but also 
to out-of-band impedances, and IM3 can be significantly affected by the 
latter. In both BJT and MESFET, collector/drain impedance at the envelope 
frequency is difficult to design, because large LC time constants are needed 
for energy storage, causing frequency-dependent envelope impedance and 
hence memory effects. These effects are more serious in BJTs, but special 
attention also has to be paid to designing the baseband Z GG in MESFETs. 
Due to C GS , this impedance is high and markedly tilted, and some amount 
of memory effect is generated at high modulation frequencies. 

Dynamic thermal effects and TPF are more important in BJTs than in 
MESFETs. The thermal impedances of the chips and packages are quite 
similar, but due to the fact that the electrical circuit elements in a BJT are 
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more sensitive to temperature, more TPF is generated in BJTs. In most 
cases, TPF has to be taken care of in a BJT, whereas in a MESFET the 
phase of the IM3 starts to be affected by dynamic temperature variations on 
the surface of the chip only when the electrical IM3 value is already very 
small. 

The total IM3 consists of a number of distortion mechanisms, and 
many of them partially cancel out each other. This tracking phenomenon is 
dependent on matching impedances and nonlinearity coefficients, which in 
turn are dependent on the bias voltages. It is often observed in practice that 
linearity suddenly improves at some value of the bias voltages, for 
example. This is caused by canceling nonlinearities, and in some situations 
the tracking is very good and good linearity is achieved. Unfortunately, 
tracking is very sensitive to changes in impedances, voltages, and 
temperature, which makes it difficult to exploit experimentally. The 
Volterra model presented in this chapter nevertheless provides a systematic 
way for studying these effects, so that it is easy to see whether or not 
tracking is possible and to which parameters the tracking is sensitive. A 
careful investigation into distortion mechanisms is the key to successful 
exploitation of the internal cancellation mechanisms. 



4.7 Key Points to Remember 

1. Accurate distortion simulations are needed in RF power amplifier 
design. 

2. In general, simulation models can be divided into two classes: 
behavioral and device models. Device models may be based on 
predefined functions, or just on tabulated measured data. 

3. Accurate derivatives have traditionally been only a secondary 
requirement in developing new device models. 

4. The derivatives up to the order of N of the I-V and Q-V curves must be 
accurate enough for Mh-order distortion simulations. 

5. Using the polynomial Volterra model, the derivatives can be locally 
fitted into the actual behavior. 

6. The Volterra analysis can give - either analytically or numerically - 
the response of each distortion mechanism separately. This makes it 
possible to look at the total distortion in a termwise manner that gives 
a lot of information for design optimization. 
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7. There are many partially canceling distortion mechanisms in RF 
power amplifiers. For example, the g pi and g m nonlinearities cancel 
each other in a current-driven BJT. 

8. PA designers have very limited control to the intrinsic nonlinearities 
of the device, but they can affect the transfer functions that convert the 
distortion currents generated by the nonlinearities to node voltages. 

9. The IM3 caused by cascaded quadratic nonlinearities is visible in both 
BJT and MESFET amplifiers, especially in BJTs. 

10. TPF is more important in BJT/HBT than MESFET amplifiers, because 
the circuit elements of the BJT are more temperature-dependent. 

1 1 .The Volterra model is a powerful tool for recognizing the memory 
effects of the amplifier. 

12. Low- frequency bias impedances cause most of the electrical memory 
effects. 



References 

[1] Maas, S., “How to model intermodulation distortion,” 1991 IEEE MTT-S 
International Microwave Symposium Digest, Vol. l,pp. 149-151. 

[2] Wambacq, P., and W. Sansen, Distortion Analysis of Analog Integrated Circuits, 
Norwell, MA: Kluwer, 1998. 

[3] Webster, D., et al., “Effect of model derivative discontinuities on cold FET 
distortion simulations,” 1997 Workshop on High Performance Electron Devices 
for Microwave and Optoelectronic Applications, pp. 97-102. 

[4] Tsividis, Y., and K. Suyama, “MOSFET modeling for analog circuit CAD: 
problems and prospects,” IEEE Journal of Solid-State Circuits, Vol. 29, No. 3, 
1994, pp. 210-216. 

[5] Root, D., “Principles and procedures for successful large-signal measurement- 
based FET modeling for power amplifier design,” Gain Without Pain seminar 
material, Agilent Technologies, 2000. 

[6] Le Gallou, N., et al., “An improvement behavioral modeling technique for high 
power amplifiers with memory,” 2001 IEEE International Microwave 
Symposium, Phoenix, AZ. 

[7] Kenington, P.B., High Linearity RF Amplifier Design, Norwood, MA: Artech 
House, 2000. 




The Volterra Model 



119 



[8] Chen, J., “What the nonlinear K-model does and how to make it do more,” 
Cadence Design Systems, 2000, at http://www.cadence.com/datasheets/ 
rf_notes.html 

[9] Verbeyst, F., and V. Bossche, “VIOMAP, the S-parameter equivalent for weakly 
nonlinear RF microwave devices,” IEEE Trans, on Microwave Theory and 
Techniques, Vol. 42, No. 12, 1994, pp. 2531-2535. 

[10] Wang, T., and T. Brazil, “A Volterra mapping-based S-parameter behavioral 
model for nonlinear RF and microwave circuits and systems,” 1999 IEEE MTT-S 
International Microwave Symposium Digest, Vol. 2, pp. 783-786. 

[11] Vuolevi, J., “Analysis, measurement and cancellation of the bandwidth and 
amplitude dependence of intermodulation distortion in RF power amplifiers,” 
Doctoral thesis, University of Oulu, Oulu, Finland, 2001. 

[12] McAndrew, C., and L. Nagel, “Early effect modelling in SPICE,” IEEE Journal 
of Solid-State Circuits, Vol. 31, No. 1, 1996, pp. 136-138. 

[13] HSPICE User’s Manual Release 96.1, Meta-Software Inc., 1996. 

[14] Laker, K., and W. Sansen, Design of Analog Integrated Circuits and Systems, 
New York: McGraw-Hill, 1994. 

[15] Versleijen, M., and A. Bauvin, “Accuracy of bipolar compact models under RF 
power operating conditions,” Proc. 1994 IEEE MTT-S International Microwave 
Symposium Digest, Vol. 3, pp. 1583-1586. 

[16] De Graaff, H., et al., “Experience with the new compact MEXTRAM model for 
bipolar transistors,” Proc. 1989 Bipolar Circuits and Technology Meeting, pp. 
246-249. 

[17] De Vreede, N., et al., “Advanced modeling of distortion effects in bipolar 
transistors using the Mextram model,” IEEE Journal of Solid-State Circuits, Vol. 
31, No. 1, 1996, pp. 114-121. 

[18] Van Rijs, F., et al., “RF power large signal modeling with MEXTRAM,” Proc. 
1996 Bipolar/BiCMOS Circuits and Technology Meeting, pp. 57-60. 

[19] Kloosterman, W., J. Geelen, and D. Klaassen, “Efficient parameter extraction 
for the MEXTRAM model,” Proc. 1995 Bipolar/BiCMOS Circuits and 
Technology Meeting, pp. 70-73. 

[20] McAndrew, C., et al., “VBIC95, the vertical bipolar inter-company model,” 
IEEE Journal of Solid-State Circuits, Vol. 31, No. 10, 1996, pp. 1476-1483. 

[21] Curtice, W., et al., “A new dynamic electro-thermal nonlinear model for silicon 
RF LDMOS FETs,” 1999 IEEE MTT-S International Microwave Symposium 
Digest, Vol. 2, pp. 419-422. 

[22] Microwave Office™ User’s Manual II, Applied Wave Research, Inc., 2000. 

[23] Kolding, T., and T. Larsen, “High order Volterra series analysis using parallel 
computing,” International Journal of Circuit Theory and Applications, Vol. 25, 
No. 2, 1997, pp. 107-114. 




120 



Distortion in RF Power Amplifiers 



[24] Heiskanen, A., and T. Rahkonen, “5th order multi-tone Yolterra simulator with 
component-level output,” Proc. 2002 IEEE International Symposium on Circuits 
and Systems, Phoenix, AZ, 2002, pp. 591-594. 

[25] Bin, L., and S. Prasad, “Intermodulation analysis of the collector-up 
InGaAs/InAlAs/InP HBT using Volterra series,” IEEE Trans, on Microwave 
Theory and Techniques, Vol. 46, No. 9, 1998, pp. 1321-1323. 

[26] Lee, J., et al., “Intermodulation mechanism and linearization of AlGaAs/GaAs 
HBT’s,” IEEE Trans, on Microwave Theory and Techniques, Vol. 45, No. 12, 
1997, pp. 2065-2072. 

[27] Fong, K., and R. Meyer, “Fligh-frequency nonlinearity analysis of common- 
emitter and differential-pair transconductance stages,” IEEE Journal of Solid- 
State Circuits, Vol. 33, No. 4, 1998, pp. 548-555. 

[28] Crosmun, A., and S. Maas, “Minimization of intermodulation distortion in GaAs 
MESFET small-signal amplifiers,” IEEE Trans, on Microwave Theory and 
Techniques, Vol. 37, No. 9, 1989, pp. 1411-1417. 

[29] Rahkonen, T., and J. Vuolevi, “Term-wise Volterra analysis of nonlinear analog 
circuits for recognising dominant cause of distortion,” Proc. Norchip 2000 
Conference, Turku, Finland, November 6-7, 2000, pp. 198-203. 

[30] Vuolevi, J., and T. Rahkonen, “The effects of source impedance on the linearity 
of BJT common-emitter amplifiers,” Proc. 2000 IEEE International Symposium 
on Circuits and Systems, Geneva, Switzerland, 2000, pp. IV-197-IV-200. 

[31] Quarles, T. et al., SPICE3 Version 3f3 User’s Manual, University of California, 
Berkeley, CA, 1993. 

[32] BGF II /X NPN 2 GHz RF Power Transistor Datasheet, Philips Semiconductors 
1995 

[33] Yamada, H., et al., “Self-linearizing technique for L-band F1BT power amplifier: 
effect of source impedance on phase distortion,” IEEE Trans, on Microwave 
Theory and Techniques, Vol. 44, No. 12,2 1996, pp. 2398-2402. 

[34] Aparin, V., and C. Persico, “Effect of out-of-band terminations on 
intermodulation distortion in common-emitter circuits,” 1999 IEEE MTT-S 
International Microwave Symposium Digest, Vol. 3, pp. 977-980. 

[35] Sevic, J., K. Burger, and M. Steer, “A novel envelope-termination load-pull 
method for ACPR optimization of RF/microwave power amplifiers,” 1998 IEEE 
MTT-S International Microwave Symposium Digest, Vol. 2, pp. 723-726. 

[36] Staudinger, J., “The importance of sub-harmonic frequency terminations in 
modelling spectral regrowth from CW AM-AM & AM-PM derived non- 
linearities,” Proc. 1997 IEEE Wireless Communications Conference, pp. 121- 
125. 

[37] CLY 2 GaAs Power MESFET Datasheet, Infineon Technologies, 1996. 

[38] Vuolevi, J., and T. Rahkonen, “Extraction of nonlinear AC FET model using 
small-signal S parameters,” IEEE Trans, on Microwave Theory and 
Measurements, Vol. 50, No. 5, May 2002, pp. 1311-1315. 




The Volterra Model 



121 



[39] Tarui, Y., et al., “An improvement of IM and power of high power amplifiers 
using RC-paralleled circuits with frequency selectivity,” 1998 IEEE MTT-S 
International Microwave Symposium Digest, Vol. 3, pp. 1655-1658. 

[40] Takenaka, I., et ah, “Improved IMD characteristics in L/S-band GaAs FET 
power amplifiers by lowering drain bias circuit impedance,” IEICE Transactions 
on Electron, Vol. 82, No. 5, 1999, pp. 730-736. 

[41] Le Gallou, N., et al., “Analysis of low frequency memory and influence on solid 
state HPA intermodulation characteristics,” Proc. 2001 IEEE International 
Microwave Symposium, Phoenix, AZ. 




Chapter 5 

Characterization of Volterra Models 



In this chapter we will see how the parameters of Volterra models for RF 
power transistors can be extracted from measured data. Various indirect 
techniques have been used to build Volterra models, and for example the 
nonlinearities of g m and g 0 [1], and even the input-output cross-terms [2] 
have been estimated from the measured level of the harmonics. Some time- 
domain characterization methods have also been published [3,4], However, 
figures like harmonic distortion lump the effects of several nonlinearities, 
and we would like to build separate electrothermal models for each 
nonlinear I-V and Q-V source. Hence, the methods used here are quite 
similar to the techniques used in the extraction of any function based 
nonlinear device models. The main differences compared to conventional 
small-signal device characterization are: 

• Power devices suffer from serious self-heating, and changing the drain 
bias voltages affects the power dissipation and hence the junction tem- 
perature. As the bias and the temperature vary simultaneously, these 
effects are difficult to separate. To overcome this problem, the self- 
heating can be minimized by arranging pulsed (isothermal) measure- 
ments. 

• Volterra models are fitted locally around the desired bias point, and 
there is no need to model the entire range of bias voltages. Instead, the 
fitting range can be chosen according to the expected signal swing. 
Note also that the numerical sensitivity to small measurement errors 
can be affected by the location of the measurement points. 

• Semiconductor manufacturers usually measure unpackaged dies. If we 
are characterizing packaged devices, de-embedding techniques are 
needed to remove the effect of the package from the measurements. 
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This chapter starts with a review of polynomial fitting techniques in 
Section 5.1. Then, the effects of self-heating and pulsed measuring 
techniques are presented in Section 5.2. I-V nonlinearities can be 
characterized by dc current measurements, and the required measurement 
setups and fitting techniques are presented in Section 5.3. 

Q-V nonlinearities must be fitted using measured capacitance values, 
as charge cannot be measured using the normal measuring instruments. The 
ac measurements are usually made at RF frequencies and include several 
important topics that are discussed from Section 5.4 onwards. In more 
detail, pulsed S-parameter measurements and fixture calibration techniques 
are discussed in Section 5.5, removing the effects of a package from the 
measured data in Section 5.6, extracting the circuit element values from the 
measured T-parameters in Section 5.7, and finally, fitting Volterra models 
based on the dQ/dV and dl/dV data obtained from the ac measurement in 
Section 5.8. 

As examples, the measured and fitted results of three power transistors 
are presented. The extracted models for a 1W BJT and MESFET and a 30- 
W LDMOS device are discussed in Sections 5.9, 5.10, and 5.11, 
respectively. 



5.1 Fitting Polynomial Models 

5.1.1 Exact and LMSE Fitting 

In the Volterra analysis, a nonlinearity J(x) is presented as a series 
expansion (5.1) around the desired operating point x 0 : 

f(x ) = f(x Q ) + a 1 ■ (x - x Q ) + a 2 • (x - x Q ) 2 + ... (5.1) 

where 



(5.2) 

X = X Q 

However, it is by no means necessary to find the coefficients «o" a N by 
calculating the higher derivatives. Instead, normal polynomial fitting can be 
performed to directly find the coefficients. If we have exactly K equations, 
the fitting results will be exact at each measured point. However, the degree 



k\ dx* 



fix) 
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of the polynomial should not be too high compared to the actual 
nonlinearity to be modeled. This is especially true if the original data is 
noisy, because errors in data points can cause slight oscillations in the 
fitting function between the data points. To avoid this, it is safe to use as 
low a degree polynomial as possible to model the nonlinearity accurately 
enough and always check the fitted result visually. If we have more data 
points, a least-mean-square error (LMSE) fitting can be performed. This 
result may not fit exactly to any of the measured points, but it is usually less 
sensitive to small errors in the measured data. 

Polynomial fitting can be easily described using matrix operations. To 
fit a Nth -degree model using N results measured at points X|,x 2 ,...,x N , we 
group the measured results /(xj) into a A/xl matrix Y and describe the model 
as a Nx(K+l) matrix M containing the different powers of xyi, 1=1... N and 
j=Q ...K, calculated at each measurement point i, and a (K+\ )xl coefficient 
vector A containing the coefficients a 0 to a K that need to be solved. Now 
the matrix equation (5.3) describes the system to be solved: 

M A = Y (5.3) 

where 





1 x, .. xf 




a 0 




/(* i) 


M = 


1 x 2 .. xf 


A = 


a i 


Y 


f(x 2 ) 




1 X N .. x§ 




a K 




/(%) 



If now N=K+ 1, an exact solution for a 0 -a K is obtained simply by 

A = M~ l ■ Y . (5.5) 

If the number of data points is larger than K+ 1, a LMSE fit is achieved by 

A = ( M T ■ M)~ l • ( M T ■ Y) (5.6) 

The contents of the matrix M is not limited to a one-dimensional 
polynomial. Instead, it may contain the selected powers and cross-products 
of F gs , E ds and temperature, for example. However, the measurement 
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points Xj must be chosen so that the matrix M will not become ill- 
conditioned, because the relative error of solution A is the relative error of 
measurements Y multiplied by a term proportional to the condition number 
of the model matrix M. 

Note that the Volterra analysis considers ac responses only. Hence, all 
voltages used in the M matrix are ac values, that is, v - V-Vq, where Vq is 
the bias point where the fit is performed. The fitted function fix) may 
consist of large signal values like currents or charges (if measurable), or be 
obtained from small-signal measurements, like capacitances or 
conductances. 

5.1.2 Effects of Fitting Range 

As the correct bias point must be obtained before the polynomial fitting can 
be applied, also the fitting range has to be correct. The limit values for the 
fitting range can be summarized as follows: 

1 . Large enough to encompass the nonlinearities; 

2. Small enough to avoid effects lying outside the signal swing; 

3. The ratio between the ranges must be correct in the case of 
multidimensional fitting. 

The first requirement is obvious, because nonlinearities are difficult to 
detect if too narrow a range is selected, and a larger range will reduce the 
sensitivity to numerical errors. But, if the range is too large, nonlinear 
effects outside the signal swing will start to affect the extracted 
coefficients, even if they have no impact on the electrical distortion. This is 
illustrated in Figure 5.1(a), where the thick line represents the load line set 
by optimum i? L , the ellipse a dynamic load line, and the rectangular box the 
fitting range. The nonlinearities of the transistor are important only within 
this range, but if the nonlinearities are extracted over a larger range, 
nonlinear effects such as saturation will cause errors in the extracted 
coefficients and lead to erroneous or less accurate simulations. 

The third requirement, illustrated in Figure 5.1(b), arises from two or 
more dimensional nonlinearities. The load resistance describes the slope of 
the load line corresponding to the ratio between the voltage and current 
swings. However, points cannot be chosen only from the load line, as this 
results in an insolvable group of equations. If all drain current 
measurements are chosen from a single load line so that = — -4 V V GS’ the 
model functions like vqs 3 and v D § 3 become linearly dependent (as v D g 3 = 
— ^v 3y GS 3 ) an d the matrix equation (5.3) will be insolvable. 
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too wide 




Figure 5.1 (a) A fitting range that is too wide in relation to the load line, and (b) 

good and poor fitting ranges in multidimensional fitting. From [5], 



Hence, it is necessary to scatter the measurements to a wider area than 
just one line. However, if the drain voltage fitting range is too wide and the 
gate voltage too small, as drafted in Figure 5.1(b), some nonlinear effects 
of the output conductance outside the actual signal swing will be taken into 
account, but all nonlinear effects of the transconductance inside the actual 
signal swing will not be seen. The fitting range in both directions must 
correspond to the extent of the signal swing, as indicated in Figure 5.1(b). 

As a conclusion, the fitting range should equal the actual signal swing 
in all dimensions. From a practical point of view, these requirements for a 
Iq- Kfip- Fce curve can be fulfilled by two parameters, output power and 
load resistance. Signal swings in different directions can be calculated from 
these two parameters, and throughout this book, a box limited by the input 
and output voltage swings is used as a fitting area. This is assumed to 
provide a good approximation whenever the nonlinearities are so weak that 
they can be modeled by the third- or fifth-degree polynomials used here. 
Alternatively, a tilted area following the estimated load line [drawn with 
dashed line in Figure 5.1(a)] can also be used. 



5.2 Self-Heating Effects 

Nonlinear extraction is impeded by variations in chip temperature. Because 
the junction temperature responds to changes in biasing conditions, the 
circuit elements of the transistor model at different bias points will be 
extracted at different chip temperatures. The chip temperature of the biased 
and zero-input device can be calculated by 
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'/'.I - T a + R tu .V d -I d , (5.7) 

where T A is the ambient temperature and R JH is the thermal resistance 
(thermal impedance at dc) describing the temperature rise caused by 
dissipated power. In the characterization measurements, the power of the 
measured signal is small so it can be assumed that the dissipated power 
equals V U I D . As a result of (5.7), the constant chip temperature contours 
presented in Figure 5.2 are obtained. The temperature variations as a 
function of the bias values cause problems to the extraction of electrical 
and thermal nonlinearity coefficients. Let us illustrate this by rewriting 
(4.3) using just the two terms below: 

*D = So ■ V D + ^3GOT ' V D ' h • (5-8) 




Figure 5.2 Constant temperature contours and extracted I-V curve with a CW (solid) 
or pulsed (dashed) signal source. From [5], 



For simplicity, let us further consider g 0 to be zero and A^got to be 
negative. Now isothermal measurements yield the dashed line shown in 
Figure 5.2. This is a horizontal line with respect to v D , corresponding to a 
zero g Q , as expected. However, if CW measurements are used, the solid line 
with a negative slope is obtained. This is due to the negative K^qqj, which 
apparently seems to decrease the value of g 0 . Actually this effect is due to 
changes in the chip temperature (proportional to v D z D ), and if the 
temperature effects are not taken into account during the extraction, g 0 
wrongly derives a negative value. In other words, the value of the nonlinear 
output conductance changes with respect to both v D and temperature, and 
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in CW measurements, these two effects are difficult, although not 
completely impossible, to separate from each other. Similar problems are 
encountered also with the other model elements, not just with the output 
conductance g 0 . 

Separation of electrical and thermal effects is required, because in 
power amplifier applications the drain voltage varies at the RF frequency, 
but the temperature changes only a little at the modulation frequency. 
Therefore, the terminal voltages and the junction temperature have to be 
independent variables, as discussed in Chapter 4. It was noticed that it 
would be possible to derive the terms for (5.8) using CW measurements 
also, if R xh is accurately known. In some works the problems of self- 
heating are handled by measuring the transistor in a steady-state condition 
and then mathematically separating the effects of self-heating from the 
purely electrical behavior [6, 7], Unfortunately, this is not possible in most 
power transistors because the dc power consumption can cause thermal 
breakdown at high bias values, destroying the device. This is the other 
reason why pulsed measurements are commonly used. 

5.2.1 Pulsed Measurements 

To avoid self-heating, the device must either be actively cooled, or, more 
simply, the dc bias must be pulsed with a low duty cycle to keep the 
average temperature constant [8, 9], The effects of the width and duty cycle 
of the biasing pulse in pulsed measurements will be discussed next. 

First of all, the measuring equipment must be fast enough to be able to 
measure during the pulsed bias. The pulse must be long enough to produce 
an electrical steady state, while at the same time it must be as short as 
possible to keep self-heating insignificant. 

The range of the thermal time constants is illustrated with two 
examples. First, the thermal time constants of a packaged 1-W CLY2 
MESFET transistor is measured by triggering the gate voltage to zero (on) 
and then monitoring the drain current as a function of time at ambient 
temperature of 20°C with an oscilloscope. Any drift in the drain current 
indicates changes in the chip temperature, and thermal settling times of up 
to 3 seconds are seen in Figure 5.3. If the package is immersed in water 
(modeling a perfect heat sink), the largest visible settling time is 100 ms. In 
both cases, 1 ms can be considered a good measurement time, because 80% 
to 95% of the change in / DS occur after this time. A similar order of 
magnitude of results is obtained in [10] by simulations. The optimum pulse 
length is dependent on the transistor type and package, so that the result 
cannot be generalized. To be sure that the self-heating can be neglected 
without considering the type, size, or transition frequency of the transistor, 
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the pulse should be shorter than the 1 ms used here. However, very fast 
measurements may cause problems with electrical settling. Thus, the 
measurement time should be as short as possible, but without causing 
electrical settling problems. 




Figure 5.3 Thermal step response of a CLY2 chip and package after triggering on 
the gate bias. © IEEE 2002 [11], 




Figure 5.4 Small-signal S 21 of an LDMOS transistor as a function of time after 
triggering on the bias voltage. © IEEE 2002 [12], 






